c语言刷题 牛顿迭代法

用牛顿迭代法求根。方程为ax^3+bx^2+cx+d=0，系数a、b、c、d由主函数输入，求x在1附近的一个实根。求出根后，由主函数输出。
相邻两次迭代结果之差的绝对值<1e-5就结束迭代。
程序的运行示例如下：
输入方程的系数a、b、c、d：
1 2 3 4
方程是： 1.00x^3+ 2.00x^2+ 3.00x+ 4.00=0
X=-1.6506292
 ***输入数据提示信息***："\n输入方程的系数a、b、c、d：\n"
 ***输入数据格式***："%f%f%f%f"
 ***输出数据格式***："\n方程是：%5.2fx^3+%5.2fx^2+%5.2fx+%5.2f=0" 
                   "\nX=%10.7f\n"

注意牛顿迭代法和二分法是两种逼近方法！

f(Xn)=a*x^3+b*x^2+c*x+d

f ' (Xn)=3a*x^2+2b*x+c

Xn+1=Xn - f(Xn)/f ' (Xn)

代码：

#include <stdio.h>
#include <math.h>
double f(double x, double a,double b, double c,double d);
double f1(double x, double a,double b, double c);
int main()
{
    double a, b, c, d, x1 ,x=1;
    printf("\n输入方程的系数a、b、c、d：\n");
    scanf("%lf%lf%lf%lf", &a, &b, &c, &d);
    printf("\n方程是：%5.2fx^3+%5.2fx^2+%5.2fx+%5.2f=0" , a, b, c, d);
    do
    {
       x=x-f(x,a,b,c,d)/f1(x,a,b,c);
       x1=x-f(x,a,b,c,d)/f1(x,a,b,c);
    }
    while (fabs(x1-x) >= 1e-5);
    printf("\nX=%10.7f\n", x1);
    return 0;
}
double f(double x, double a,double b, double c,double d)
{
    double n;
    n = a * x * x * x + b * x * x + c * x + d;
    return n;
}
double f1(double x, double a,double b, double c)
{
    double n;
    n=3*a*x*x+2*b*x+c;
    return n;
}