请设计二分法算法，求方程f(x)=x^3-x-1=0在区间【1,1.5】内的解（精确到0.01）

二分法的计算过程：
1.x=1 时，f(1)=-1<0,f(1.5)=0.875>0
2.[1,1.5]的中点为1.25，f(1.25)=-0.297<0
3.那么区间变为[1.25,1.5],|1.5-1.25|=0.25>0.01,继续
4.[1.25,1.5]的中点为1.375，f(1.375)=0.224>0
5.那么区间变为[1.25,1.375],|1.375-1.25|=0.13>0.01，继续
6.[1.25,1.375]的中点为1.313，f(1.313)=-0.050<0
7.那么区间变为[1.313,1.375],|1.375-1.313|=0.06>0.01，继续
8.[1.313,1.375]的中点为1.344,f(1.344)=0.084>0
9.那么区间变为[1.313,1.344],|1.344-1.313|=0.03>0.01，继续
10.[1.313,1.344]的中点为1.329,f(1.329)=0.018>0
11.那么区间变为[1.313,1.329],|1.329-1.313|=0.02>0.01，继续
12.[1.313,1.329]的中点为1.321,f(1.321)=-0.016<0
13.那么区间变为[1.321,1.329],|1.321-1.329|=0.01
14.得出方程的解为（1.321+1.329）/2=1.325

#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);
}
}

1.324718