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 }