计算方法——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;
}

计算结果：