c语言实现牛顿迭代法

 1 #include<stdio.h>
 2 #include<stdlib.h>  
 3 #include<math.h>  
 4 #include<float.h>
 5 #include<time.h>
 6 
 7 #define PI 3.14159265358979323846  /* pi */
 8 #define ε 1.0e-12
 9 int main()
10 {
11     double x0 = PI;//取的初始值
12     double x1 = 0.0;//有x0算出的x1，初始值先给定0
13     double fx = 0.0;//f（x）
14     double fxp = 0.0;//f（x）的导数
15     double faix = 0.0;//计算结果,牛顿迭代格式 faix =x-(fx/fxp)
16     int i = 0;//迭代计算次数
17     while (fabs((x0 - x1) / x1) > ε)//判断两次迭代变量之间相对误差与精度比较
18     {
19         x1 = x0;//用x1存放上次的x0
20         fx = sin(x0) - x0 / 2;
21         fxp = cos(x0) - 0.5;
22         faix = x0 - fx / fxp;
23         x0 = faix;//将迭代后的结果赋给上次代入值变量，供下一次代入使用
24         i++;//计算次数
25         printf("第%d次迭代，迭代结果为：,%+.12e  \n", i, x1);
26     }
27 }

题目：计算sinx=x/2的根。

分析：newton法在大范围的收敛定理：

函数f（x）在区间[a，b]上存在二阶连续导数，且满足4个条件：

1.　　f（a）*f（b)<0

2.　　当x属于[a,b]时，函数的导数值不等于零。

3.　　当x属于[a,b]时，函数的二阶导数值保号。

4.　　a-f（a）/f'（a）<=b,且b-f（b）/f'（b）<=a

 

 

 

 计算结果：

 

转载于:https://www.cnblogs.com/zhubinglong/p/6044763.html