新手求教一个求SIN函数定积分的问题 谢谢各位了

请问下  我只是想求 0到π这个范围下sinx的定积分  这是我的算法
 double integral(double a,double b,double (*fun)(double))
{
const double n=1000.0;
double l,sum=0.0;
l=(b-a)/n;
for(int i=1;i<=n;i++)
sum+=l*fun(a+i*l);
return sum;
}   请问错在哪   谢谢了   a是函数的下限 b是函数的上限 l是分成1000等分 每一等分的长度  我是把图形分成1000个矩形 然后算矩形的面积和的方式求解的