c语言非线性方程根迭代法,计算方法——C语言实现——牛顿迭代法求解非线性方程...

最近在上计算方法这门课，要求是用MATLAB做练习题，但是我觉得C语言也很棒棒啊~

问题：利用牛顿迭代法计算X3+2X2+10X-20=0的根，精度为0.000001

使用VS2017，代码如下：

/*牛顿迭代法求解非线性方程的根*/

#include "stdafx.h"

#include "math.h"

#define tol 0.000001 //精度

/*由方程导出的函数在此，返回计算值*/

double fun(double fun_a)

{

fun_a = fun_a * fun_a*fun_a + 2*fun_a* fun_a + 10* fun_a - 20;

return fun_a;

}

/*由方程导出的函数的导数在此，返回计算值*/

double fun1(double fun_a)

{

fun_a = 3 * fun_a* fun_a + 4 * fun_a + 10;

return fun_a;

}

int main()

{

unsigned int c_times = 0;

double c_result,p;

p = 1.5;//给一个摸约是根的值，值约接近，迭代越快

printf("以下利用牛顿迭代法计算X^3+2X^2+10X-20=0的根\t精度为%g\n",tol);

printf("取初值p为：%g\n",p);

for (c_times = 1; c_times < 1000; c_times++) // 迭代最多1000次，不行就退出

{

c_result = p - fun(p)/fun1(p);//迭代方程

printf("当前计算次数为：%d\t计算结果为：%g\n", c_times, c_result);

if (fabs(p - c_result) < tol) break;

p = c_result;

if (c_times == 999)

{

printf("计算方法失效！");

goto _end;

}

}

printf("最终计算得结果为：%g\t所需计算次数为：%d\n",c_result,c_times);

_end:

printf("计算完成，按回车退出程序。");

getchar();

return 0;

}

运行结果如下：