写出用二分法求方程x^3-x-1=0在区间[1,1.5]上的一个解的算法误差不超过0.001并画出程序框图及程序越...
发布网友
发布时间:2022-05-27 15:40
我来回答
共2个回答
热心网友
时间:2023-10-31 12:59
#include<stdio.h>
#include<math.h>
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);
}
}
程序适合a*x^3+b*x^2+c*x+d=0,精度也很高,修改下就行
热心网友
时间:2023-10-31 12:59
见图。