C语言算法---牛顿迭代法求三次方程的根

一、牛顿迭代法
在方程代表的函数上随便取一个点，过该点做函数的切线，切线交x轴于Xn，过Xn做x轴的垂线交f（x）于另一点，再过该点做f（x）的切线，并重复以上过程，即有公式Xn+1=Xn-f（xn）/f '(xn);当Xn+1于Xn差距在0.001范围时，认为其已足够接近函数的根，取Xn作为根的近似。

二、代码
//a,b,c,d分别为三次方程从高阶到低阶的系数

float fun(int a,int b,int c,int d){
do{
float x0=1,f,f1,x; //取初始随机点为1，f为该点的函数值即f（x），f1为函数在该点的导数,x为下一近似点
f=a*x0*x0*x0+b*x0*x0+c*x0+d;  //求函数值
f1=3*a*x0*X0+2*b*x0+c;  //求导数
x=x0-f/f1;
}while(fabs(x-x0)>=1e-3)//在指定精度内求解
return x;
}