나는 함수를 통합하려고한다. 그래서 이런 짓을 확인하기 위해, 내가 무엇을 기대 내가지고있어 결과가 없습니다Scipy 적분이 잘못 나온다
def integrand(a, b):
return scipy.integrate.quad(function, a, b)
: 나는 그것을 통합하고있어
이제def function(x):
something = ...
something_else = ...
return exp(something)/sqrt(something_else)
:이 기능은 비 부정적 보장
for x in range(0, 10000):
if integrand(0,x+1) < integrand(0,x):
raise ValueError("Weird!")
물론, '이상한'예외가 있습니다. 어떻게 그렇게 될수 있니?
문제를 나타내는 실제 예제를 표시 할 수 있습니까? – BrenBarn
실제 예제는 10,000 포인트 C14 데이트 보정 곡선이 포함되어 있기 때문에 어렵습니다. – zmbq
작은 테스트 케이스 없이는 문제를 디버그하기가 어렵습니다. 문제를 작은 테스트 케이스로 줄이기 위해 노력하십시오. – BrenBarn