非线性方程的数值解法：二分法 python

最新推荐文章于 2023-10-15 08:59:31 发布

厄巧克力

最新推荐文章于 2023-10-15 08:59:31 发布

阅读量603

点赞数 1

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

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

版权

计算方法专栏收录该内容

11 篇文章 2 订阅

订阅专栏

import sympy
import math
x=sympy.symbols("x")
yuan=input("请输入原函数：")
duishu=sympy.log(x)
a=float(input('请输入a的值='))
b=float(input('请输入b的值='))
k=1
yaoqiu=int(input('请输入采用第几种精度要求：'))
yibusen=float(input('精度为='))
N=(sympy.log(b-a)-sympy.log(yibusen))/sympy.log(2)
N=math.ceil(N)

while True:
    X=(a+b)/2
    h=(b-a)/2
    fx=sympy.sympify(yuan).subs(x,X)
    if yaoqiu == 1 :
        if abs(fx)<yibusen:
            print("二分次数%s\n零点x=%s\n函数值f(x)=%s\n"%(k,X,fx))
            break
    else:
        if h<yibusen:
            print("二分次数%s\n零点x=%s\n函数值f(x)=%s\n"%(k,X,fx))
            break
    fa=sympy.sympify(yuan).subs(x,a)
    if fa*fx<0:
         b=X	
    else:
         a=X
    if k>=N: 
        print("分半N0次还没有达到精度要求信息")
        break       
    k+=1