先上总结,紧跟代码,后上讲解
总结
之前总结不完整,完整的以后敬上
Python代码(重修订)
注:代码中的下标正好与上面推论中的相反
import cmath
import math
def f(x):
return x ** 5 - 3.5 * x ** 4 + 2.75 * x ** 3 + 2.125 * x ** 2 - 3.875 * x + 1.25
rate = 0.01
r = s = -1
a = [1, -3.5, 2.75, 2.125, -3.875, 1.25]
b = []
c = []
def deal() :
for i in range(len(b)) :
a[i] = b[i]
del a[len(a) - 2 : len(a)]
flag= 0
for k in range(9999):
size = len(a)
del b[0: len(b)]
b.append(a[0])
b.append(a[1] + r * b[0])
for i in range(2, size):
b.append(a[i] + r * b[i - 1] + s * b[i - 2])
del c[0: len(c)]
c.append(b[0])
c.append(b[1] + r * c[0])
for i in range(2, size - 1):
c.append(b[i] + r * c[i - 1] + s * c[i - 2])
R = (b[size -