牛顿迭代法求方程根

/*
题目：牛顿迭代法求方程根
内容：编写用牛顿迭代法求方程根的函数。方程为ax3+bx2+cx+d = 0,系数a,b,c,d由主函
      数输入。求x在1附近的一个实根。求出根后，由主函数输出。
	  牛顿迭代法的公式是：x = x0 - f(x0)/fd(x0),设迭代到|x-x0|<=1e-5时结束。
作答者：小白杨
收获： 
*/ 

#include <stdio.h>

float fib(float,float,float,float);

int main()
{
	float a,b,c,d,x;
	printf("请输入方程的系数：\n");
	scanf("%f%f%f%f",&a,&b,&c,&d);
	x = fib(a,b,c,d);
	printf("该方程在1附近的根为：%f\n",x);
	
	return 0;
}

float fib(float a,float b,float c,float d)
{
	float x = 1,x0,h,f,fd;
	do{
		x0 = x;
		f = a*x*x*x+b*x*x+c*x+d;
		fd = 3*a*x*x+2*b*x+c;
		h = f/fd;
		x = x-h;
	}while(fabs(x-x0)>=1e-5);
	return x;
}