第33题
有如下数学公式:
已知梯形法求积分公式为:
(其中n为积分区间的等分数)函数trap是一个利用梯形法求定积分的通用求积分函数。double pexp( )是计算公式一的函数,double poly( )是计算公式二的函数,请根据以下调用语句,完成trap函数中的填空。
调用语句:
y1=trap(pexp,0.0,1.0)/sqrt(2.0*3.1416);
y2=trap(sin,0.0,3.1416/2.0)/2.0;
y3=trap(polyt,0.0,4.0);
double trap(__1___,double a,double b)
{
double t,h;
int i,n=1000;
t=((*fun)(a)+(*fun)(b))/2.0;
h=fabs(a-b)/(double)(n);
for(i=1;i<=n-1;i++)
t+=__2__;
t*=h;
return t;
}