Xhungab

|+ Calcul de int(f(x)) en fonction des bornes a et b avec n itérations . C'est la deuxième partie du calcul.

double trapezoidsimpson(

if(i ==0 || i== n){m = 1.;}

else if(fmod(i,2) == 0){m = 2.;}

return( ((b -a)*M) / (23*n) );

double trapezoid_dydxsimpson_dydx(

int ny)

if(i ==0 || i== n){m = 1.;}

else if(fmod(i,2) == 0){m = 2.;}

return( ((b-a)*M) / (23*n) );

Dans les deux premières colonnes, il y a la fonction de référence pour calculer une intégrale simple par la méthode desde trapèzesSympson. Dans les deuxièmes colonnes il y a deux fonctions pour calculer l'intégrale double. L'intégrale double est calculé par la fonction trapezoid_dydxsympson_dydx(); qui appelle la fonction int_dy();

'''int_dy();''' Cette fonction applique la méthode desde trapèzesSympson pour la variable y.

'''trapezoid_dydx();''' Cette fonction applique la méthode desde trapèzesSympson pour la variable x.

|+ Calcul de int(f(y)) en fonction des bornes ay et by avec ny itérations

double trapezoidsimpson(

if(i ==0 || i== n){m = 1.;}

else if(fmod(i,2) == 0){m = 2.;}

return( ((b -a)*M) / (23*n) );

if(i ==0 || i== ny){m = 1.;}

else if(fmod(i,2) == 0){m = 2.;}

M += m * (*P_f)( x, ay+i*(by-ay)/ny);

return( ((by-ay)*M) / (23*ny) );