double integ(double a , double b , int n)//積分区間１、積分区間２、分割数
{
int i;
double s,xl,xr,dx;
s=0,0; //積分区間の面積の総和
dx=(b-a)/(double)n; //分割区間の長さ
for(i=0;i<n;++i){
xl = a + dx * (double)i; //y[n]のx座標
xr = xl + dx; 　//y[n+1]のx座標
s = s + (fn(xl)+fn(xr)) * dx / 2.0; // (上底+下底)*高さ/2
}
return(s);
}
計算部分のみ載せました。変なとこあったらご指導お願いします