非线性方程的数值解法：Aitken加速 python

最新推荐文章于 2024-08-13 18:29:56 发布

厄巧克力

最新推荐文章于 2024-08-13 18:29:56 发布

阅读量528

点赞数

分类专栏：计算方法文章标签： python 开发语言算法

本文链接：https://blog.csdn.net/weixin_49530111/article/details/126285094

版权

计算方法专栏收录该内容

11 篇文章 2 订阅

订阅专栏

import sympy
from sympy import exp    #用来输入指数函数，详见输入案例
x=sympy.symbols("x")
yuan=input("请输入原函数：")
yuan=sympy.simplify(yuan)
x0=float(input("请输入x0的值："))
N0=int(input('请输入N0(最大迭代次数)：'))
yibusen=float(input("请输入精度要求："))
k=-1  
while True:
    x1=yuan.subs(x,x0)
    x2=yuan.subs(x,x1)
    X=x0-(x1-x0)**2/(x2-2*x1+x0)
    if abs(x1-x0)<yibusen:
        print('x=',X,'(k=',k,')')
        break
    x0=X
    if k>N0:
        print("method failed after N0 iterations,N0=",N0)
        break
    k+=1