用二分法求函数f(x)=x^3-x-1在区间[1,1.5]内的一个零点(精确到0.01)
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/08 17:09:39
用二分法求函数f(x)=x^3-x-1在区间[1,1.5]内的一个零点(精确到0.01)
用二分法求函数f(x)=x^3-x-1在区间[1,1.5]内的一个零点(精确到0.01)
用二分法求函数f(x)=x^3-x-1在区间[1,1.5]内的一个零点(精确到0.01)
f(1)=-10
f(1.25)=-0.2968750
所以又在(1.25,1.375)之间
f(1.3125)=-0.05150
所以在(1.3125,1.34375)之间
f(1.328125)=0.014575>0
所以在(1.3125,1.328125)之间
f(1.3203125)=-0.018701
f(1)=-1<0
f(1.5)=0.875>0
f(1.25)=-0.296875<0
所以在(1.25,1.5)之间
f(1.375)=0.2246>0
所以又在(1.25,1.375)之间
f(1.3125)=-0.0515<0
所以在(1.3125,1.375)之间
f(1.34375)=0.0826>0
...
全部展开
f(1)=-1<0
f(1.5)=0.875>0
f(1.25)=-0.296875<0
所以在(1.25,1.5)之间
f(1.375)=0.2246>0
所以又在(1.25,1.375)之间
f(1.3125)=-0.0515<0
所以在(1.3125,1.375)之间
f(1.34375)=0.0826>0
所以在(1.3125,1.34375)之间
f(1.328125)=0.014575>0
所以在(1.3125,1.328125)之间
f(1.3203125)=-0.018701<0
所以在(1.3203125,1.328125)之间
此时区间间距已小于0.01,所以就是这个区间
收起
#include
#include
double fx(double a,double b,double c,double d,double x1,double x2,double f1,double f2)
{
double x0,f;
while(fabs(f)>1e-6)
{
x0=(x1+x...
全部展开
#include
#include
double fx(double a,double b,double c,double d,double x1,double x2,double f1,double f2)
{
double x0,f;
while(fabs(f)>1e-6)
{
x0=(x1+x2)/2;
f=a*pow(x0,3)+b*x0*x0+c*x0+d;
if(f1*f>0)
{
x1=x0;
f1=f;
}
else
{
x2=x0;
f2=f;
}
}
return x0;
}
void main()
{
double a,b,c,d,x,x1,x2,f1,f2;
printf("请输入方程系数 a,b,c,d:\n");
scanf("%lf,%lf,%lf,%lf",&a,&b,&c,&d);
printf("请输入区间[x1,x2]\n");
scanf("%lf,%lf",&x1,&x2);
f1=a*pow(x1,3)+b*x1*x1+c*x1+d;
f2=a*pow(x2,3)+b*x2*x2+c*x2+d;
if(f1*f2>0)
{
printf("该区间内方程无解\n");
}
else
{
x=fx(a,b,c,d,x1,x2,f1,f2);
printf("方程解为x=%lf\n",x);
}
}
收起