C语言编程-8_8 公式求解 牛顿迭代法求解方程

编写函数，用牛顿迭代法求方程f(x)=2x3-4x2+3x-6=0在1.5附近的根。

牛顿迭代公式为：xn+1=xn-f(xn)/f'(xn)

其中，f'(xn) 是f在xn处的导数。

结束条件：|f(xn+1)|< eps与|xn+1-xn|< eps同时成立（eps是一个很小的正数，从键盘输入）

同时编写主函数，在主函数中调用并输出函数值。

函数原型如下：

double root(double (*f)(double),double (*f1)(double),double x,double eps);

其中，形参f接收原函数f(x)，形参f1接收导函数f'(x)

如输入：

1e-8

如输出结果为：

2.000000

注意 浮点数的绝对值函数是 fabs,整数的是abs，不能混用。

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
// f（x） =  2x^3 -4x^2 +3x - 6 = 0 在1.5处的根
// f'(x) =  6x^2 - 8x + 3

double fun(double x)
{
	x = 2 * x*x*x - 4 * x*x + 3 * x - 6;
	return x;
}

double fun1(double x)
{
	x = 6 * x*x - 8 * x + 3;
	return x;
}

double root(double(*f)(double), double(*f1)(double),double x, double eps)
{
	double x0;
	
	while (1)
	{
		x0 = x;
		x = x0 - (*f)(x) / (*f1)(x);
		if (fabs((*f)(x)) < eps && fabs(x - x0) < eps)
		{
			break;
		}
	}
	return x;
}

int main()
{
	double x, eps;
	x = 1.5;
	scanf_s("%lf", &eps);
	x = root(fun,fun1,x, eps);
	printf("%lf", x);
	return 0;
}