Mathc initiation/c36a4

Installer ce fichier dans votre répertoire de travail.

	x_xzy.h

/* ---------------------------------- */
/* save as x_xzy.h                   */
/* --------------------------------- */
double intx_dxdzdy(

double (*P_f)(double x, double y, double z),

double (*Px0)(double y, double z),
double (*Px1)(double y, double z),
   int nx,
   
double z,

double y
)
{
   int i = 0;
double m = 0.;
double M = 0.;

 for(i = 0; i <= nx; i++)
 {
       if(i ==0 || i== nx){m = 1.;}
  else if(fmod(i,2) ==  0){m = 2.;}
  else                    {m = 4.;}

  M += m * (*P_f)( ((*Px0)(y,z)) + i*(((*Px1)(y,z))-((*Px0)(y,z)))/nx, y, z);
 }

  return( ((((*Px1)(y,z)) -((*Px0)(y,z)))*M) / (3*nx) );
}
/* --------------------------------- */
double intz_dxdzdy(
double (*P_f)(double x, double y, double z),
   
double (*Px0)(double y, double z),
double (*Px1)(double y, double z),
   int    nx,

double (*Pz0)(double x),
double (*Pz1)(double x),
   int    nz,
   
double y
)
{
   int i = 0;
double m = 0.;
double M = 0.;

 for(i = 0; i <= nz; i++)
 {
       if(i ==0 || i== nz){m = 1.;}
  else if(fmod(i,2) ==  0){m = 2.;}
  else                    {m = 4.;}

  M += m *  intx_dxdzdy(
       (*P_f),    
       (*Px0),(*Px1),nx,          
     (((*Pz0)(y))+i*(((*Pz1)(y))-((*Pz0)(y)))/nz),          
          y);
 }

 return( ((((*Pz1)(y)) -((*Pz0)(y)))*M) / (3*nz) );
}
/* --------------------------------- */
double simpson_dxdzdy(
double (*P_f)(double x, double y, double z), 
   
double (*Px0)(double y, double z),
double (*Px1)(double y, double z),
int    nx,

double (*Pz0)(double x),
double (*Pz1)(double x),
   int    nz,


double y0,
double y1,
   int ny
)
{
   int i = 0;
double m = 0.;
double M = 0.;

 for(i = 0; i <= ny; i++)
 {
       if(i ==0 || i== ny){m = 1.;}
  else if(fmod(i,2) ==  0){m = 2.;}
  else                    {m = 4.;}

  M += m * intz_dxdzdy(
       (*P_f),         
       (*Px0),(*Px1),nx,       
       (*Pz0),(*Pz1),nz, 
       (y0+i*(y1-y0)/ny));
 }

  return( ((y1 -y0)*M) / (3*ny) );
}
/* --------------------------------- */
/* --------------------------------- */

Les fonctions pour calculer les intégrales de la forme dxdzdy