下面是一个书上的例子 全是照抄的 但运行居然有错（有时候会是个死循环）

题目：用弦截法求方程f(x)=x*x*x-5*x*x+16*x-80=0的解．

书上给的答案是：（看看那里有问题）

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

float f(float x)
{
float y;
y=((x-5.0)*x+16.0)*x-80.0;
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;
}

int main(void)
{
float x1,x2,f1,f2,x;
do
{
printf("input x1.x2:");
scanf("%f,%f",&x1,&x2);
f1=f(x1);
f2=f(x2);
}
while(f1*f2>=0);
x=root(x1,x2);

printf("A root of equation is %8.4f.",x);

system("pause");
return 0;
}