第八章 例8.6

/* 用弦截法求方程f(x)=x^3-5x^2+16x-80=0的根 */

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

float f(float x)					/* 定义f函数，实现f(x)=x^3-5x^2+16x-80 */
{
	float y;
	y = ((x - 5.0) * x + 16.0) * x -80;

	return y;
}

float xpoint(float x1, float x2)	/* 定义xpoint函数，求弦与x轴的交点 */
{
	float y;
	y = (x1 * f(x2) - x2 * f(x1)) / (f(x2) - f(x1));

	return y;
}

float root(float x1, float x2)		/* 定义root函数，求近似根 */
{
	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.000001);

	return x;
}

int main()
{
	float x1, x2, f1, f2, x;
	do
	{
		printf("Please input x1, x2:\n");
		scanf_s("%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\n\n", x);
	
	return 0;
}