无条件稳定显格式

[mapengfei-nwpu]

# newton ecant method
# xian ge fa

# default interval is [0,1]
def secant(f, x0=0, x1=1):
    if abs(x0 - x1) < 1e-7:
        return x1
    else :
        f0 = f(x0)
        f1 = f(x1)
        k = (f1-f0)/(x1-x0)
        # print("k",k,"x0",x0)
        if abs(k) < 1e-7:
            return secant(f, x1, 0.5*(x1+x0))
        x = x0-f0/k
        return secant(f, x1, x)

if __name__ == "__main__":
    def f(x):
        return x*x

    print(secant(f,99,100))
    print(secant(f))

[mapengfei-nwpu]

IBFE的能量守恒.pdf