essai01.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 *X,h;
  printf("n=");scanf("%d",&n);
  DisRegAlloc(0,1,n,&X,&h);
  printf("h = %lf\n",h);
  printDoubleDynArray(X,n,"X","%8.5lf",5);
  printf("Int=%lf\n",quadPM(f,X,n,h));
  Free(X);
  return EXIT_SUCCESS; // return 0
}