用弦解法求方程y=a*x*x*x+b*x*x+c*x+d根

#include
#include
int a,b,c,d;
float f(float x)
{float y;
y=((a*x+b)*x+c)*x+d;
return(y);
}
float xpoint(float x1,float x2)
{float y;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
return(y);
}
float root(float x1,float x2)
{float x,y,y1;
y1=f(x1);
do{x=xpoint(x1,x2);
y=f(x);
if(y*y1>=0)
{y1=y;
x1=x;
}
else
x2=x;
}while(fabs(y)>=0.0001);
return(x);
}
void main()
{float x1,x2,f1,f2,x;
printf("input a,b,c,d");
scanf("%d%d%d%d",&a,&b,&c,&d);
do{ptintf("input x1,x2");
scanf("%f%f",&x1,&x2);
f1=f(x1);
f2=f(x2);
} while(f1*f2>=0);
x=root(x1,x2);
printf("the num is %f",x);
}
不知道有没有错误地方
有谁看出来请多多指教啊