我有以下简单的 python 程序来使用二分法求方程根:
from numpy import exp
#
def fun(x):
return 5.0+4.0*x-exp(x)
#
a=3
b=10.0
eps=1.0e-15
#
fa=fun(a)
fb=fun(b)
#
if fa*fb>0:
print("wrong interval!!!",fa,fb)
exit()
#
iter=1
while (b-a)>eps:
c=(a+b)/2.0
fc=fun(c)
if fc==0:
print("x = ",c)
exit()
if fc*fa>0:
a=c
fa=fc
else:
b=c
fb=fc
iter+=1
#
print("x = ",c)
print("accuracy = ",'{:.2e}'.format(b-a))
print("f(",c,") =",fun(c))
print(iter," iterations needed")
如果我将 a 放入错误的间隔(例如 a=3),它会说这是错误的间隔,但无论如何它都会继续给出(显然)错误的结果和四行
错误:root:别名无效:名称 less 不能使用别名,因为它是另一个魔术命令。
而且,内核死掉了(我使用的是 jupyter)。你能帮助我吗?
撒科打诨
交互式爱情
相关分类