欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

(C语言)用二分法求方程 2x^3 - 4x^2 + 3x - 6 = 0在(-10, 10)之间的根

程序员文章站 2024-02-02 15:20:52
...

用二分法求方程 2x^3 - 4x^2 + 3x - 6 = 0在(-10, 10)之间的根

二分法,又称分半法,是一种方程式根的近似值求法。对于区间 [a,b] 上连续不断且 f(a) × f(b) < 0 的函数y=f(x),通过不断地把函数f(x)的零点所在的区间一分为二,使区间的两个端点逐步逼近零点,进而得到零点近似值

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define N 50 //二分次数
double fun(double x1, double x2) {
	int count = 0;
	double f1, f2, fm, mid;
	while (1) {
		f1 = 2 * pow(x1, 3) - 4 * pow(x1, 2) + 3 * x1 - 6;
		f2 = 2 * pow(x2, 3) - 4 * pow(x2, 2) + 3 * x2 - 6;
		if (f1 * f2 > 0) {
			return 0;
		}
		++count;
		mid = (x1 + x2) / 2;
		fm = 2 * pow(mid, 3) - 4 * pow(mid, 2) + 3 * mid - 6;
		if (fm * f1 == 0 || fm * f2 == 0 || count == N) {//循环终止条件,(找出根的准确值或者已达二分次数)
			break;
		}
		if (fm * f1 < 0) {
			x2 = mid;
		}
		else if (fm * f2 < 0) {
			x1 = mid;
		}
	}
	return mid;
}
void main() {
	int x1 = -10;
	int x2 = 10;
	if (fun(x1, x2) != 0) {
		printf("方程2x^3 - 4x^2 + 3x - 6 = 0在区间(%d,%d)之间的根为:%f\n", x1, x2, fun(x1, x2));
	}
	else {
		printf("方程2x^3 - 4x^2 + 3x - 6 = 0在区间(%d,%d)之间无根\n",x1, x2);
	}
	system("pause");
}

(C语言)用二分法求方程 2x^3 - 4x^2 + 3x - 6 = 0在(-10, 10)之间的根