非线性方程的数值解法：正割法 python

厄巧克力

于 2022-08-11 14:46:08 发布

阅读量547

点赞数

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

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

版权

计算方法专栏收录该内容

11 篇文章 2 订阅

订阅专栏

import sympy
from sympy import exp    #用来输入指数函数，详见输入案例
x=sympy.symbols("x")
yuan=input("请输入原函数：")
yuan=sympy.simplify(yuan)
x0=float(input('请输入x0的值='))
x1=float(input('请输入x1的值='))
yibusen=float(input("请输入根的误差限为："))
k=1
xx=[x0,x1]
while True:
    xx.append(0)
    f=yuan.subs(x,xx[k-1])
    f1=yuan.subs(x,xx[k])
    xx[k+1]=xx[k]-(f1/(f1-f))*(xx[k]-xx[k-1])
    if abs(xx[-1]-xx[-2])<yibusen:
        print(xx[-1],'(k=',k+2,')')
        break
    k+=1