我编写了以下代码来使用辛普森规则来集成任意函数,在我的例子中是 sin(x):
import math
delta_x=(x2-x1)/N
def simpson(x1,x2,f,N):
sum=0
i=1
for i in range(1,N+1):
sum+=f(x1+i*delta_x)
sum1=(3*delta_x)*sum
return(sum1)
print(simpson(0,math.pi,math.sin(x),100))
但是我在 sum+=f(x1+i*delta_x) 行上收到错误“float object not callable”。有谁知道可能出了什么问题?
谢谢 :)