c语言写方程次数,计算方法——C语言实现——二分法求解非线性方程

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

问题：利用二分法计算X^3-X-1=0的根，精度为0.001

使用VS2017，代码如下：

/*二分法求解非线性方程的根*/

#include "stdafx.h"

#include "math.h"

#define tol 0.001 //精度

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

double fun(double fun_a)

{

fun_a = fun_a * fun_a*fun_a - fun_a - 1;//X^3-X-1

return fun_a;

}

int main()

{

unsigned int c_times = 0;

double c_result = 0,c_result1,a,b;

a = 1;

b = 2;//大致含*根区间在1到2

printf("以下利用二分法计算X^3-X-1=0的根\t精度为%g\n\n",tol);

for (c_times = 1; c_times < 1000; c_times++)//二分1000次，得不到结果就退出

{

c_result1 = c_result;

c_result = (a + b) / 2;//二分

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

if (0 == fun(c_result) || fabs(c_result-c_result1) < tol) break;//满足精度，退出

if (fun(a)*fun(c_result) < 0) b = c_result;//重新划定边界

else a = c_result;

if (c_times == 999)

{

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

goto _end;

}

}

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

_end:

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

getchar();

return 0;

}

计算结果：