landmark_gzz 2022-04-25 11:25 采纳率: 100%
浏览 244
已结题

Python关于不动点迭代法相关问题

这是原题
用不动点迭代法求f(x)=e^x+x-2=0在[0,1]内的近似根,要求误差小于10-6.
这是代码部分
 from math import exp
   xl = 0  #区间下限
   xu = 1  #区间上限
   x = (xl+xu)/2  #迭代初始值
   x_list = [x]
   i = 0

 while True:
x = exp(x) + x - 2 + x
x_list.append(x)
if len(x_list) > 1:
    i += 1
    error = abs((x_list[-1] - x_list[-2]) / x_list[-1])
    if error < 10**(-6):
        print(f'迭代第{i}次后,误差小于10^-6')
        break
else:
    pass
这是他的报错
 OverflowError: math range error

也就是溢出,数学范围错误

代码部分我是参考的Fo*(Bi)的
 https://blog.csdn.net/weixin_48615832/article/details/109470488
不知该如何解决
  • 写回答

2条回答 默认 最新

  • landmark_gzz 2022-04-25 13:12
    关注

    以下为我修改后成功运行的代码

    from math import exp,log
    xl = 0  #区间下限
    xu = 1  #区间上限
    x = (xl+xu)/2  #迭代初始值
    x_list = [x]
    i = 0
    
    while True:
        x = log(2-x)
        x_list.append(x)
        if len(x_list) > 1:
            i += 1
            error = abs((x_list[-1] - x_list[-2]) / x_list[-1])
            if error < 10**(-6):
                print(f'迭代第{i}次后,误差小于10^-6')
                break
        else:
            pass
    print(f'所求方程式的根为{x_list[-1]}') 
    
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论 编辑记录
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 5月10日
  • 已采纳回答 5月2日
  • 创建了问题 4月25日

悬赏问题

  • ¥15 echarts动画效果失效的问题。官网下载的例子。
  • ¥60 许可证msc licensing软件报错显示已有相同版本软件,但是下一步显示无法读取日志目录。
  • ¥15 Attention is all you need 的代码运行
  • ¥15 一个服务器已经有一个系统了如果用usb再装一个系统,原来的系统会被覆盖掉吗
  • ¥15 使用esm_msa1_t12_100M_UR50S蛋白质语言模型进行零样本预测时,终端显示出了sequence handled的进度条,但是并不出结果就自动终止回到命令提示行了是怎么回事:
  • ¥15 前置放大电路与功率放大电路相连放大倍数出现问题
  • ¥30 关于<main>标签页面跳转的问题
  • ¥80 部署运行web自动化项目
  • ¥15 腾讯云如何建立同一个项目中物模型之间的联系
  • ¥30 VMware 云桌面水印如何添加