欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

牛顿迭代法python求x^3-x-1=0的根

程序员文章站 2024-02-02 21:40:22
...
def f(xi):
    return xi*xi*xi-xi-1
def f1(xi):
    return 3*xi*xi-1
x=[]
x.append(0.5)
eps=1e-14     #误差限制
error=abs(f(x[-1]))    #最新加的x在最后
number_iteration=0
while error>eps:
    x.append(x[-1]-f(x[-1])/f1(x[-1]))  #x k+1
    error=abs(f(x[-1]))
    number_iteration=number_iteration+1
print('牛顿法迭代次数为%f次'%(number_iteration)) #格式化输出
print('方程的根x*为%f'%(x[-1]))
print('f(x*)的值为%f'%(f(x[-1])))