最近在上计算方法这门课,要求是用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;
}
运行结果如下: