在下面的代码中,我在Python中实现了对分方法。正如一般概述,我的代码执行以下操作:我的函数能够找到任意连续标量值函数f的根,作为lambda函数的一个给定的函数容差。在
我的例程在区间[0,1]上查找函数f(x)=cos x-sin x的根,容差为10^{14}
记录达到此公差所需的迭代次数。在
但是现在我想在同一区间绘制收敛图。这将是迭代次数的绝对误差。
为了做到这一点,我必须在一个列表中收集一系列错误数字,并将其与一系列整数1到iter的最终值进行比较。在
我正在寻找一些帮助,因为我被困了。我已经用不同的迭代方法做了另外两个代码,所以一旦我能看到它如何在这一个我应该能够实现它在其他人也!非常感谢所有的帮助import math
def root(x):
return(math.cos(x)-math.sin(x))
def bisection_method(f, a, b, tol):
if f(a)*f(b) > 0:
#end function, no root.
print("No root found.")
else:
iter = 0
while (b - a)/2.0 > tol:
midpoint = (a + b)/2.0
if f(a)*f(midpoint) < 0: # Increasing but below 0 case
b = midpoint
else:
a = midpoint
iter += 1
return(midpoint, iter)
answer, iterations = bisection_method(root, 0, 1, 10**(-14))
print("Answer:", answer, "\nfound in", iterations, "iterations")