8.12 用牛頓迭代法求根。方程為：ax^3+bx^2+cx+d=0 ，系數a，b，c，d由主函數輸入。求X在1附近的一個實根。求出后由主函數輸出.

 

 

//https://baike.baidu.com/item/%E7%89%9B%E9%A1%BF%E8%BF%AD%E4%BB%A3%E6%B3%95/10887580?fr=aladdin#4
//百度牛頓迭代法    
    #include <stdio.h>
    #include <math.h>

    double solut(double a,double b,double c,double d)
    {
        double x1=1,x,f,f1;
        
        //迭代
        do
        {   x=x1;
            f=((a*x+b)*x+c)*x+d;//原函數
            f1=(3*a*x+2*b)*x+c;//求導之后的函數
               x1=x-f/f1;        
        } while(fabs(x1-x)>=0.00000001);
        return x1;
    }

    int main()
    {
        double solut(double ,double ,double ,double );
        double a,b,c,d;
        scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
        printf("%.2f\n",solut(a,b,c,d));
        return 0;
    }

結果：