二分法解一元三次方程的C++程序自动出结果改进自我检测功能.docVIP

#includestdio.h #includemath.h void main() { double x0,x1,xm,f0,f1,fm,x2,x3,g=10;//x2,x3是驻点，x0,x1,xm,f0,x1是二分法求根的工具。g代替人手工输入。 double a[3],r[3]; int i,j=0; printf(input 3 coefficients:\n); for(i=0;i3;i++) { printf(a[%d]=,i); scanf(%lf,a+i); } printf(the function is:\n); printf(x*x*x+%5.2f*x*x+%5.2f*x+%5.2f=0\n,a[0],a[1],a[2]); //y=3*x*x+2*a[0]*x+a[1] if(4*a[0]*a[0]-12*a[1]0)//方程单调递增，与横轴只有一个交点。 { x0=-g; x1=g; f0=x0*x0*x0+a[0]*x0*x0+a[1]*x0+a[2]; f1=x1*x1*x1+a[0]*x1*x1+a[1]*x1+a[2]; if(f0*f1==0)//正好-g和g有一个是根 { if(f0==0)//x0是根，x1不是。 { xm=x0; printf(the only root of the function is %.6f\n,xm); } else//x1是根，x0不是。 { xm=x1; printf(the only root of the function is %.6f\n,xm); } if(fabs((xm*xm*xm+a[0]*xm*xm+a[1]*xm+a[2])/(3*xm*xm+2*a[0]*xm+a[1]))1e-6) printf(comgratulations! the answer is right!\n); else printf(sorry, you should try again.\n); } else { while(f0*f1=0) { g*=10; x0=-g; x1=g; f0=x0*x0*x0+a[0]*x0*x0+a[1]*x0+a[2]; f1=x1*x1*x1+a[0]*x1*x1+a[1]*x1+a[2]; } do { xm=(x0+x1)/2; fm=xm*xm*xm+a[0]*xm*xm+a[1]*xm+a[2]; if(f0*fm0) x0=xm; else x1=xm; }while(fabs(x0-x1)1e-6); xm=(x0+x1)/2; printf(the only root of the function is %.6f\n,xm); if(fabs((xm*xm*xm+a[0]*xm*xm+a[1]*xm+a[2])/(3*xm*xm+2*a[0]*xm+a[1]))1e-6) printf(congratulations! the answer is right!\n); else printf(sorry, you should try again.\n); } } else//方程有增有减，但与横轴的交点不确定。 { x2=(-2*a[0]-sqrt(4*a[0]*a[0]-12*a[1]))/6; x3=(-2*a[0]+sqrt(4*a[0]*a[0]-12*a[1]))/6; printf(the stagnation of the function are:\n); printf(x2=%.6f\nx3=%.6f\n,x2,x3); if((x2*x2*x2+a[0]*x2*x2+a[1]*x2+a[2