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