Mathc initiation/003g

Installer ce fichier dans votre répertoire de travail.

	x_sphbar.h

/* --------------------------------- */
/* save as x_sphbar.h                */
/* --------------------------------- */
double spheric_x_dbdadr(

double (*P_f)(double a, double b, double r),

double (*Pb0)(double a, double r),
double (*Pb1)(double a, double r),
   int   nb,
   
double    a,
   
double    r
)
{
double b = 0;
   int i = 0;
double m = 0.;
double M = 0.;

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

  b = ((*Pb0)(a,r)) + i*(((*Pb1)(a,r))-((*Pb0)(a,r)))/nb;
  
  M += m * (*P_f)(a,b,r) * r*r * sin(b);
 }

  return( ((((*Pb1)(a,r)) -((*Pb0)(a,r)))*M) / (3*nb) );
}
/* --------------------------------- */
double spheric_b_dbdadr(
double (*P_f)(double a, double b, double r),
   
double (*Pb0)(double a, double r),
double (*Pb1)(double a, double r),
   int    nb,

double (*Pa0)(double r),
double (*Pa1)(double r),
   int    na,
   
double     r
)
{
   int i = 0;
double m = 0.;
double M = 0.;

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

  M += m * spheric_x_dbdadr(  (*P_f),
                           
                             (*Pb0),
                             (*Pb1),
                               nb,
                               
                           (((*Pa0)(r))+i*(((*Pa1)(r))-((*Pa0)(r)))/na),
                               
                              r);
 }

 return( ((((*Pa1)(r)) -((*Pa0)(r)))*M) / (3*na) );
}
/* --------------------------------- */
double spheric_dbdadr(
double (*P_f)(double a, double b, double r),
   
double (*Pb0)(double a, double r),
double (*Pb1)(double a, double r),
   int    nb,
   
double (*Pa0)(double r),
double (*Pa1)(double r),
   int    na,

double r0,
double r1,
   int nr
)
{
   int i = 0;
double m = 0.;
double M = 0.;

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

  M += m *  spheric_b_dbdadr((*P_f),
                               
                             (*Pb0),
                             (*Pb1),
                                nb,
                             
                             (*Pa0),
                             (*Pa1),
                                na,
                                
                             (r0+i*(r1-r0)/nr));
 }

  return( ((r1 -r0)*M) / (3*nr) );
}
/* --------------------------------- */
/* --------------------------------- */