定积分计算（谭浩强c语言第5版p272）

    本算法用C语言实现定积分的计算：利用“指向函数的指针”作为integral函数的形参，在调用integral函数时，将f函数送给指针形参，通过指针调用f函数并输出结果，计算过程通过f函数实现，计算结果通过"return()"返回到integral函数。

#include <stdio.h>
#include <math.h>
int main()
{
	int integral(double, double, double (*p)(double,double));
	double f1(double, double);
	double f2(double, double);
	double f3(double, double);
	double f4(double, double);
	double f5(double, double);
	int i;
	double a,b;
	printf("你想计算以下哪个函数的定积分？\n1.1+x  2.2x+3  3.e^x+1  4.(1+x)^2  5.x^3\n");
	scanf_s("%d", &i);
	printf("请输入积分下限a和积分上限b:\n");
	scanf_s("%lf,%lf", &a, &b);
	switch (i)
	{
	case 1:
		integral(a, b, f1);
		break;
	case 2:
		integral(a, b, f2);
		break;
	case 3:
		integral(a, b, f3);
		break;
	case 4:
		integral(a, b, f4);
		break;
	case 5:
		integral(a, b, f5);
		break;
	}
	return 0;
}
int integral(double a, double b, double (*p)(double, double))
{
	double result;
	result = (*p)(a, b);
	printf("%f", result);
	return 0;
}
double f1(double a, double b)
{
	double c=(b*b/2.0+b)-(a*a/2.0+a);
	printf("函数“1+x”在区间(%f,%f)的定积分为：",a,b);
	return(c);
}
double f2(double a, double b)
{
	double c = (b*b + 3.0*b) - (a*a + 3.0*a);
	printf("函数“2x+3”在区间(%f,%f)的定积分为：", a, b);
	return(c);
}
double f3(double a, double b)
{
	double c = (exp(b)+b) - (exp(a)+a);
	printf("函数“e^x+1”在区间(%f,%f)的定积分为：", a, b);
	return(c);
}
double f4(double a, double b)
{
	double c = pow(b+1.0,3)/3.0 - pow(a + 1.0, 3) / 3.0;
	printf("函数“(1+x)^2”在区间(%f,%f)的定积分为：", a, b);
	return(c);
}
double f5(double a, double b)
{
	double c = pow(b,4)/4.0 - pow(a, 4) / 4.0;
	printf("函数“x^3”在区间(%f,%f)的定积分为：", a, b);
	return(c);
}

使用头文件<math.h>以调用函数：

1."exp(x)"表示"e^x";

2."pow(x,y)"表示"x^y"。

运行结果：