二分法求根

基于二分法求方程的根,主要思路是通过给定两个初始值a和b,使f(a)f(b)<0。如果f(x)是连续的,则必然在(a,b)内存在一个值c,使f(c)=0。

       首先计算c=(a+b)/2;

       如果f(c)=0,则直接返回c即可,否则,进行下面的判断:

       如果f(a)f(c)>0,则令a=c,否则令b=c。

参考代码如下:

#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <cmath>
#include <math.h>
#include <cstring>
#include <string>
#include <queue>
#include <deque>
#include <stack>
#include <stdlib.h>
#include <list>
#include <map>
#include <utility>
#include <set>
#include <bitset>
#include <vector>
#define pi acos(-1.0)
#define inf 0x3f3f3f3f
#define linf 0x3f3f3f3f3f3f3f3fLL
#define ms(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long ll;
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}

int main()
{
//2x3-4x2+3x-6
    float a,b,mid;
    float temp1,temp2,temp3;
    a=-10.;
    b=10.;
    mid=(a+b)/2;

    temp3=2*mid*mid*mid-4*mid*mid+3*mid-6;
    while(fabs(temp3)>1e-6)
    {
        temp1=2*a*a*a-4*a*a+3*a-6;
        temp2=2*b*b*b-4*b*b+3*b-6;
        if(temp1*temp3>0)
            a=mid;
        else
            b=mid;
        mid=(a+b)/2;
        temp3=2*mid*mid*mid-4*mid*mid+3*mid-6;
    }
    printf("the solution of the equation is %.4f.\n",mid);

    return 0;
}



  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值