1:编写程序,用二分法求方程2x^3-4x^2+3x-6=0在[-10,10]之间的根,^符号表示乘幂,如x^3表示x的立方,使

问题描述:

1:编写程序,用二分法求方程2x^3-4x^2+3x-6=0在[-10,10]之间的根,^符号表示乘幂,如x^3表示x的立方,使

#include
#include
void main()
{double f(double x);/*对关于x的函数进行声明*/
void erfen(double a,double b,double h,double ep);/*对二分法的子程序进行声明*/
erfen(-10,10.0,0.100,0.00001); /*对二分法的子程序进行调用*/
}
double f(double x)/* 定义一个关于x的函数*/
{double y;
y=2*x*x*x-4*x*x+3*x-6;/*将函数值赋给y*/
return y;} /*返回函数值*/
void erfen(double a,double b,double h,double ep)/*定义二分法的子程序*/
{double x[10],a0,b0,c0,a1,b1;/*将得到的根存放在x[]中,定义区间的端点为浮点型变量*/
int k,i,j;
k=0;/*k为循环的控制变量*/
a0=a; /*将a赋给第一个区间的起点*/
while(a0{for(j=0;;j++) /*循环每执行一次,则向后移动一个区间,区间的长度为步长*/
{b0=a0+h; /*b0为第一个区间的末端*/
if(fabs(f(a0)){x[k]=a0; /*将a0的值放入数组x中*/
a0=b0; /*指向下一个区间*/
k++; /*x的地址指向下一个*/
break;} /* 跳出循环*/
if(fabs(f(b0)){a0=b0;
break;} /*如果b0的函数值小于给定精度,则跳出循环,指向下一个区间*/
if(f(a0)*f(b0)>0)
{a0=b0;
break;} /*如果a0与b0的函数值乘积大于0,则跳出循环,指向下一个区间*/
if(f(a0)*f(b0){
a1=a0;b1=b0;/*使a1始终b1为区间的端点*/
for(i=0;;i++)
{ c0=(a1+b1)/2;/*使区间缩为一半*/
if(fabs(f(c0)){x[k]=c0;
a0=b0;
k++;
break;}/* 跳出循环,指向下一个区间*/
if(f(c0)*f(a1)>0)
a1=c0;/*c0不是根,使c0的值赋给左端点*/
else
b1=c0;}/*c0不是根,是c0的值赋给右端点*/
break;}/* 跳出循环,指向下一个区间*/
}}
if(k!=0) /* k不等于0时,输出实根的值*/
for(i=0;iprintf("%f,%f\n",x[i],f(x[i]));
else
printf("no root"); /*k等于0时,输出无实根*/
}