essai02.c

Go to the documentation of this file.

#include <math.h>

#include "mesh1D.h"
#include "quad.h"
#include "utils.h"

double f(double x){
   return 8.*sqrt(x*(1-x)); 
}

int main(){
  int n;
  double h,I;
  printf("n=");scanf("%d",&n);
  I=quadPMAlloc(f,0.,1,n,&h);
  printf("Int=%lf\n",I);
  printf("Error=%.16e\n",fabs(I-M_PI)/M_PI);
  return EXIT_SUCCESS; /* return 0 */
}