>>21
計算精度を考慮に入れなければ…

double S = 0.0;
double x;

x1 x2 は領域の端点
f は積分する関数

for (x = x1; x < x2; x += (x2-x1)/100)
　S += f(x) * ((x2-x1)/100);