学习笔记：定积分的求解(矩形法)

定积分求解

• 矩形法将定积分化为多个函数连续的和。 基本思想是将区间[a，b] 进行n等分 当n越大的时候结果越准确。 图形化成一小块一小块的矩形。底边长都为(b - a) / n. 高为每个等分点的函数值。 然后将每个矩形的面积相加即为所求。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

float integral(float(*func)(float), float a, float b, int n)
{
	int i;
	float x, h, s = 0;
	h = (b - a) / n;//每一个小矩形的底边长
	x = a;
	for ( i = 0; i <n ; i++)
	{
		x = x + h;
		s = s + (*func)(x) * h;//(*func)(x)为高，将每个矩形的面积相加
		//不用函数指针，法2
		//float temp = sin(x);
		//s = s + temp * h;
	}
	return s;
}
float fsin(float x)
{
	//printf("fsinx=%5.2f\n", x);
	return sin(x);
}
float fcos(float x)
{
	//printf("fcosx=%5.2f\n", x);
	return cos(x);
}
float fexp(float x)
{
	//printf("fexp=%5.2f\n", x);
	return exp(x);
}
int main()
{
	float a1, b1, a2, b2, a3, b3;
	float (*pfunc)(float);//声明一个指针p,指向一个形参为float的函数.
	float sum;
	int n = 20;
	//测试1
	//printf("Please enter a1,b1:\n");
	//scanf("%f%f", &a1,&b1);
	//pfunc = fsin;
	//sum = integral(pfunc, a1, b1, n);
	//printf("sinx的定积分为：%f\n", sum);
	//测试2
	//printf("Please enter a2,b2:\n");
	//scanf("%f%f", &a2, &b2);
	//pfunc = fcos;
	//sum = integral(pfunc, a2, b2, n);
	//printf("cosx的定积分为:%f\n", sum);
	//测试3
	printf("Please enter a3,b3:\n");
	scanf("%f%f", &a3, &b3);
	pfunc = fexp;
	sum = integral(pfunc, a3, b3, n);
	printf("cosx的定积分为:%f\n", sum);



	system("pause");
	return 0;
}