import numpy as np
def romberg(f,l):
val=np.zeros(l)
for i in range(0,l):
h=2**(1-i)
N=2/h
val[i]=trapez(f,h,N)
return val
def trapez(f,h,N):
result=0
stuetz=np.zeros(N+1);
for j in range(0,N+1):
stuetz[j]=0+j/N
sum=0
for k in range(1,N): #bilde summe von 1-N-1
sum+=np.polyval(f,stuetz[k])
a=np.polyval(f,0)*0.5*h
b=(np.polyval(f,stuetz[N]))*0.5*h
result=a+(h*sum)+b
return result
main:
p=[1,0,0]
romberg(p,5)
error:
----> 3 romberg(p,5)
----> 7 val[i]=trapez(f,h,N)
类型错误:“float”对象无法解释为整数
我该如何解决这个问题请帮助我?
我很感激!
change N=2/h
to N=2//h
执行整数除法。
def romberg(f,l):
val=np.zeros(l)
for i in range(0,l):
h=2**(1-i)
N=2//h
val[i]=trapez(f,h,N)
return val
edit:
您还可以使用int(2/h)
def romberg(f,l):
val=np.zeros(l)
for i in range(0,l):
h=2**(1-i)
N=int(2/h)
val[i]=trapez(f,h,N)
return val
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)