牛顿迭代法求方程的根

#include <stdio.h>

#include <math.h>

int

main()

{

double

a,b,c,d,n,m,x,y;

while

(

scanf

(

"%lf,%lf,%lf,%lf"

,&a,&b,&c,&d)!=EOF)

{

x=0;

m=a*x*x*x+b*x*x+c*x+d;

n=3*a*x*x+2*b*x+c;

y=x-m/n;

while

(

fabs

(x-y)>=1e-5)

{

x=y;

m=a*x*x*x+b*x*x+c*x+d;

n=3*a*x*x+2*b*x+c;

y=x-m/n;   

}

printf

(

"%.6lf\n"

,y);

}

return

0;

}

x(n+1)=x(n)－f(x(n))/f'(x(n))