essai03.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 i,nini=100,pas=100,nfin,N,nb,n;
  double *ph,*pError,Exact=M_PI,val;
  printf("nfin=");scanf("%d",&nfin);
  nb=(nfin-nini)/pas+1;
  ph=(double *)Malloc(nb*sizeof(double));
  pError=(double *)Malloc(nb*sizeof(double));
  for (i=0;i<nb;i++){
    N=nini+i*pas;
    val=quadPMAlloc(f,0.,1,N,&ph[i]);
    pError[i]=fabs(val-Exact)/Exact;
    printf("N=%8d, val=%.16lf, error=%.16e\n",N,val,pError[i]);
  }
  saveTwoDoubleDynArray(ph,pError,nb,"essai03.out","%.16e %.16e");
  Free(ph);
  Free(pError);
  return EXIT_SUCCESS; // return 0
}