最简单的二分法

区间二分法在非线性方程的数值解法中算是最简单的方法了;我这里给出的是最原始的二分法,没有经过优化的代码。

#include <iostream>

#include <cstdio>
#include <cmath>
#include <algorithm>


using namespace std;


double f(double x){
    return 1-x-sin(x);
}//这是返回要处理的函数的值, 大家可以根据自己的需求改变。
int main()
{
    double a,b,s;
    printf("请输入a,b的值,及精度\n");
    scanf("%lf%lf%lf",&a,&b,&s);
    int k = 0;
    if((f(a)*f(b))>0)
    {
        printf("请输入a,b的值,及精度\n");
        scanf("%lf%lf%lf",&a,&b,&s);
    }
    while(1){
    if((f(a)*f(b))==0)
    {
        if(f(a)==0)
        {
            printf("%lf %d\n",a,k);
            break;
        }
        else
        {
            printf("%lf %d\n",b,k);
            break;
        }
    }
    else
    {
        double m = (a+b)/2;
        double n = a-b;
        if(fabs(n)<s)
        {
            printf("结果:%lf ,迭代次数%d\n",m,k);
            break;
        }
        else
        {
            if((f(a)*f(b))>0)
            {
                a = m;
            }
            else b = m;
        }
    }
    k++;
    }//while循环是整个代码的核心。
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值