贝尔斯托法求根

先上总结，紧跟代码，后上讲解

总结

之前总结不完整，完整的以后敬上

 

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 -