意大利数学家斐波那契(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));
}