计算方法——C语言实现——二分法求解非线性方程

版权声明:本文为博主原创文章,遵循 CC 4.0 by-sa 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/you_big_father/article/details/83015165

最近在上计算方法这门课,要求是用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;
}


计算结果:二分法

展开阅读全文

没有更多推荐了,返回首页