C语言——矩形法求定积分的原理和实现

一、前言：

即将其元素化，依据积分的几何意义将其变成若干小矩形，

被积函数：f ( x )

积分区间：【a,b】

每个小矩形的面积：
S i = ( b − a ) / N ∗ ∣ f ( a + ( b − a ) / N ∗ i ) ∣ Si = (b - a )/ N * |f(a + (b - a )/ N * i)|Si=(b−a)/N∗∣f(a+(b−a)/N∗i)∣

二、代码实现：

#include<stdio.h>
#include<math.h>
#define N 1000000
double get_integral_sin(double a, double b)
{
	double x = (b - a) / N;
	double ans = 0;
	for (int i = 1; i <= N; ++i)
	{
		ans += x * sin(a + x * i);
	}
	return ans;
}
double get_integral_cos(double a, double b)
{
	
	double x = (b - a) / N;
	double ans = 0;
	for (int i = 1; i <= N; ++i)
	{
		ans += x * cos(a + x * i);
	}
	return ans;
}
double get_integral_e(double a, double b)
{
	
	double x = (b - a) / N;
	double ans = 0;
	for (int i = 1; i <= N; ++i)
	{
		ans += x * exp(a + x * i);
	}
	return ans;
}
int main() {
	printf("%lf\n", get_integral_cos(0, 2));
	printf("%lf\n", get_integral_sin(0, 2));
	printf("%lf\n", get_integral_e(0, 2));

	return 0;
}

三、运行结果：

四、注意事项：

1.不需要使用abs绝对值，因为x轴上方面积为正，下方面积为负，其和才是正确积分结果

2.abs绝对值函数为int型，对于微分矩形，算出来的值都为0