二分法求解方程的解(CUMTB)

意大利数学家斐波那契(1175~1250)研究过一个一元三次方程 x3 + 2x2 + 10x - 20 = 0的求解问题,并给出了一个精度非常高的近似解x=1.368808108,但是没有人知道他是怎么计算出来的。 
一般地,已知一元三次方程 x3 + 2x2 + 10x - c = 0 当15<=c<=35时方程在区间(1, 2)之间肯定有一个解,请你写一个程序用二分法求出此方程的解,并输出。 
要求c从键盘输入,输出此方程的解 
要求近似计算的精度达到|f(x)|<1e-9

#include<iostream>
#include<cstdio>
#include<cmath>
const double Max=1e-9;
using namespace std;
double sb(double a,double b,double c)
{
    double mid=(a+b)/2;
    double fa=a*a*a+2*a*a+10*a-c;
    double f=mid*mid*mid+2*mid*mid+10*mid-c;
    if(fabs(f)<Max) return mid;
    if(fa*f<0)
    {
       return sb(a,mid,c);
    }
    else
    {
       return sb(mid,b,c);
    }

}
int main()
{
    double c;
    cin>>c;
    if(c<15||c>35)
    {
        cout<<"input error!"<<endl;
        return 0;
    }
    printf("%.9f\n",sb(1,2,c));

}


版权声明:本文为qq_33588413原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明。