第八周上机实践项目：矩形法求定积分

*All rights reserved.
*文件名称：ipj.cpp
*作    者：赵子琳
*完成日期：2016年4月24日
*版 本 号：v1.0
*
*问题描述：写一个用矩形法求定积分的通用函数
*/

#include<stdio.h>
#include<cmath>

float integral(float (*p)(float),float a,float b,int n)
{
float x,h,s;
h=(b-a)/n;
x=a;
s=0;
int i;
for(i=1;i<=n;i++)
{
x=x+h;
s=s+(*p)(x)*h;
}
return s;
}

float fsin(float x)
{
return sin(x);
}
float fcos(float x)
{
return cos(x);
}
float fexp(float x)
{
return exp(x);
}

int main()
{
float integral(float (*p)(float),float a,float b,int n);
float a1,b1,a2,b2,a3,b3,c,(*p)(float);
float fsin(float);
float fcos(float);
float fexp(float);
int n=20;
scanf("%f%f",&a1,&b1);
scanf("%f%f",&a2,&b2);
scanf("%f%f",&a3,&b3);
p=fsin;
c=integral(p,a1,b1,n);
printf("The integral of sin(x) is :%.2f\n",c);
p=fcos;
c=integral(p,a2,b2,n);
printf("The integral of cos(x) is :%.2f\n",c);
p=fexp;
c=integral(p,a3,b3,n);
printf("The integral of exp(x) is :%.2f\n",c);
return 0;
}