编写一个求定积分的通用函数用C++,用指向函数的指针作函数参数

float integral(float (*fun)(float),float a,float b);
其中,a、b表示积分区间,fun是函数指针.
float f1(float x)
{
float f;
f = 1 + x*x;
return f;
}
float f2(float x)
{
float f;
f = 1 + x + x*x + x*x*x;
return f;
}
float f3(float x)
{
float f;
f = x / (1 + x*x);
return f;
}
float integral(float (*fun)(float),float a,float b)
{
float s,h,y;
int n,i;
s = ( (*fun)(a) + (*fun)(b) ) /2.0;
n = 100;
h = (b-a)/n;
for(i=1; i