用弦截法求解方程的根

/* 用弦截法求解方程的根 算法思想 三元方程的函数增减单调特性 定义X1,X2两个坐标
使得F(X1) 与F(X2) 的值相反 这样X1,X2之间必有一跟 由下面的弦截法公式求焦点坐标
当F(X)与F(X1)同符号时，证明F(X)>0则新区间为 X1--X 如果 F(X)与F(X2)同符号，则新
区间为X--X2 就这样无限趋近下去 最后到一个数量级很小时，那个X就近似为一元三次
方程的根 */

/*这个程序是很早很早前写的，只是觉得这个方法很不错，呵呵*/
#include"stdio.h"
#include"math.h"
float f(float x)
{
	return (x*(x-8)+12)*x-30;  //求X点出的函数值
}
float xpoint(float x1,float x2)
{
	return (x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));  //求F(X1),F(X2)连线与X轴的交点
}
float root(float x1,float x2)
{
	float x,y,y1;
	y1=f(x1);
	while(fabs(y)>=0.0001)       //当Y的绝对值足够小，即X足够接近目标根
	{
	x=xpoint(x1,x2);
	y=f(x);
	if(y1*y>0)                  //当F(X1)*F(X)>0 则新区间为X--X2，否则新区间X1--X
	{
		y1=y;
		x1=x;
	}
	else 
	    x2=x;
	}
	return x;
}
void main()
{
	float x1,x2,y;
	printf("please input x1,x2!\n");
input:	sca