一元高次方程求解,二分法

没啥好说的,就是从估计区间-100到100之间,以1为单位,一个一个从左到右拓展,遇到正整数解就直接输出,遇到在1个单位内存在实根(就  f(l) * f(r)<0 发生了,在l和r之间函数线穿越X轴 )

那么就准备用二分法求近似值,近似值准确性由while语句控制,二分法没啥好说的,根在mid左边,那么mid就是新的右边界,反过来,根在mid的右边,那么mid就是新的左边界。

新的左右边界继续二分,直到左右边界的差值小于精度就停止循环,此时认为l或r无限接近,l或r就是方程的一个近似根

for循环继续搜索1个单位,一直搜到100为止。

如果你大概预估你的根区域不止-100,100区间,就改一下for语句

如果你要调整根的精度,就改while语句

#include <bits/stdc++.h>
using namespace std;
double f(double x)
{
    return 3*pow(x,4)-5*pow(x,3)+2*x*x+7*x-120;     //根据要求改函数即可。 

void solve()
{
    for(int i=-100;i<=100;i++)  //盲猜根的区间在-100到100之间,查询每一个格子是否有实根 
    {
        double l=i,r=i+1;
        if(f(l)==0) printf("%.3lf\n",l); //发现实根就是整数
        if(f(l)*f(r)<0){        //发现在lr之间存在实根,就用二分法查到近似值 
            while(r-l>=0.0001)  //精确到第几位,小数点后就几个0 
            {
                double mid=(l+r)/2;
                if(f(l)*f(mid)<0){
                    r=mid;
                }else{
                    l=mid;
                }
            }
            printf("%.3lf\n",l);
        }
    }
    
}
int main()
{
    solve();
    return 0;
}

  • 10
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值