#include<stdio.h>
int main(){
double a3,a2,a1,a0,a,b,x,fx,fa,fb, bottle;
scanf("%lf %lf %lf %lf\n%lf %lf",&a3,&a2,&a1,&a0,&a,&b);
fa = a3 * a*a*a + a2 * a*a+ a1 *a + a0;
fb = a3 * b*b*b + a2 * b*b+ a1 *b + a0;
do{
if (b - a <= 0.01){
x = (a + b)/ 2;
break;
}else if(fa * fb < 0){
x = (a + b)/ 2;
fx = a3 * x*x*x + a2 * x*x+ a1 *x + a0;
if(fx == 0){
;
}else if(fx * fa > 0){
a = x;
x = (a + b)/ 2;
fa = a3 * a*a*a + a2 * a*a+ a1 *a + a0;
}else{
b = x;
x = (a + b)/ 2;
fb = a3 * b*b*b + a2 * b*b+ a1 *b + a0;
}
}else if(fa == 0){
x = a;
fx = a3 * x*x*x + a2 * x*x+ a1 *x + a0;
}else{
x = b;
fx = a3 * x*x*x + a2 * x*x+ a1 *x + a0;
}
}while(fx != 0);
printf("%.2f",x);
return 0;
}
二分法求多项式单根(未构造函数)
最新推荐文章于 2022-03-07 21:31:21 发布