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;
}

轉載https://blog.csdn.net/you_big_father/article/details/83015165