Mathc initiation/a467
Installer ce fichier dans votre répertoire de travail.
x_jwuv.h |
---|
/* --------------------------------- */
/* save as x_jwuv.h */
/* --------------------------------- */
double J_z_wuv(
double (*P_f)(double x, double y, double z),
double (*P_X)(double x, double y, double z),
double (*P_Y)(double x, double y, double z),
double (*P_Z)(double x, double y, double z),
double (*P_az)(double x, double y),
double (*P_bz)(double x, double y),
int nz,
double x,
double y
)
{
pt3d p;
double z = 0;
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.;}
z = ((*P_az)(x,y)) + i*(((*P_bz)(x,y))-((*P_az)(x,y)))/nz;
p.x = x;
p.y = y;
p.z = z;
M += m * (*P_f)( ((*P_X)(x,y,z)),
((*P_Y)(x,y,z)),
((*P_Z)(x,y,z)))
*
/* Jacobian */
fabs(
fxyz_x((*P_X),H,p)
*
( fxyz_y((*P_Y),H,p)*fxyz_z((*P_Z),H,p)
-
fxyz_z((*P_Y),H,p)*fxyz_y((*P_Z),H,p) )
-
fxyz_y((*P_X),H,p)
*
( fxyz_x((*P_Y),H,p)*fxyz_z((*P_Z),H,p)
-
fxyz_z((*P_Y),H,p)*fxyz_x((*P_Z),H,p) )
+
fxyz_z((*P_X),H,p)
*
( fxyz_x((*P_Y),H,p)*fxyz_y((*P_Z),H,p)
-
fxyz_y((*P_Y),H,p)*fxyz_x((*P_Z),H,p) )
);
}
return( ((((*P_bz)(x,y)) -((*P_az)(x,y)))*M) / (3*nz) );
}
/* --------------------------------- */
double J_x_wuv(
double (*P_f)(double x, double y, double z),
double (*P_X)(double x, double y, double z),
double (*P_Y)(double x, double y, double z),
double (*P_Z)(double x, double y, double z),
double (*P_az)(double x, double y),
double (*P_bz)(double x, double y),
int nz,
double (*P_ax)(double y),
double (*P_bx)(double y),
int nx,
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 * J_z_wuv( (*P_f),
(*P_X), (*P_Y),(*P_Z),
(*P_az),(*P_bz), nz,
(((*P_ax)(y))+i*(((*P_bx)(y))-((*P_ax)(y)))/nx),
y);
}
return( ((((*P_bx)(y)) -((*P_ax)(y)))*M) / (3*nx) );
}
/* --------------------------------- */
double J_wuv(
double (*P_f)(double x, double y, double z),
double (*P_X)(double x, double y, double z),
double (*P_Y)(double x, double y, double z),
double (*P_Z)(double x, double y, double z),
double (*P_az)(double x, double y),
double (*P_bz)(double x, double y),
int nz,
double (*P_ax)(double y),
double (*P_bx)(double y),
int nx,
double ay,
double by,
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 * J_x_wuv( (*P_f),
(*P_X), (*P_Y),(*P_Z),
(*P_az),(*P_bz),nz,
(*P_ax),(*P_bx),nx,
(ay+i*(by-ay)/ny));
}
return( ((by -ay)*M) / (3*ny) );
}
/* --------------------------------- */
/* --------------------------------- */
----