Integration.c

#include "mesh1D.h"

/* Sur une discretisation reguliere */
double QuadPMreg(double (*f)(double), double *X, int N, double h){
  int i;
  double S=0.;
  for (i=0;i<N-1;i++)
    S+=f((X[i]+X[i+1])/2.);
  return S*h;
}