如何用C语言编写一矩形法求定积分的程序？

double intergral(double (* fx)(double x), double a, double b) {
const double dx = 0.000001;
double toReturn = .0;

while(fabs(a - b) >= dx) {
toReturn += dx * fx(a);
a += dx;
}

return toReturn;
}

/* MinGW下通过 */

跟矩形法求定积分的数学思路差不多吧