牛顿迭代法

算法描述如下:

给定初始值x0，ε为根的容许误差，η为 | f(x) | 的容许误差，N为迭代次数的容许值。

1.如果 f'(x0)=0 或迭代次数大于N，则算法失败，结束；否则执行2。

2.计算x1=x0-f(x0)/f'(x0)。

3.若| x1-x0 |<ε 或| f(x1) |<η ，则输出x1，程序结束；否则执行4。

4.令x0=x1，转向1

 

程序实例(c语言实现)

求方程xe(x) - 1 = 0 的根

1. /*牛顿迭代法*/
2. #include <stdio.h>
3. #include <math.h>
4. #define N 100
5. #define eps 1e-6
6. #define eta 1e-8
8. float Newton(float(*f)(float),float(*fl)(float),float x0)
9. {
10.     float x1,d;
11.     int k=0;
12.     do
13.     {
14.         x1=x0-(*f)(x0)/(*fl)(x0);
15.         if(k++>N || fabs((*fl)(x1))<eps)
16.         {
17.             printf("/n Newton 迭代发散");
18.             break;
19.         }
20.         d=fabs(x1)<1?(x1-x0):(x1-x0)/x1;  //计算误差值
21.         x0=x1;
22.         printf("x(%d)=%f/t",k,x0);
23.     }
24.     while(fabs(d)>eps && fabs((*f)(x1))>eta);
25.     return x1;
26. }
28. float f(float x)
29. {
30.     return x*exp(x)-1;
31. }
33. float fl(float x)
34. {
35.     return exp(x)+x*exp(x);
36. }
38. void main()
39. {
40.     float x0,y0;
41.     printf("请输入迭代初值 x0/n");
42.     scanf("%f",&x0);
43.     printf("x0=%f/n",x0);
44.     y0=Newton(f,fl,x0);
45.     printf("方程的根为: %f/n",y0);
46. }

输出:

请输入迭代初值 x0

1.0

x0=1.000000

x(1)=0.683940   x(2)=0.577454   x(3)=0.567230   x(4)=0.567143   x(5)=0.567143

方程的根为: 0.567143