Open mapengfei-nwpu opened 3 years ago
# 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))