Mathc initiation/c36ba
Installer ce fichier dans votre répertoire de travail.
/* save as fa.h */
double f(
double x,
double y,
double z)
 return( (y*z) );
char feq[] = "y*z";
double ux(
double y,
double z)
 return(0);
char ueq[] = "0";
double vx(
double y,
double z)
 return((1-y-z));
char veq[] = "1-y-z";
double sy(
double z)
 return((0));
char seq[] = "0";
double ty(
double z)
 return((1-z));
char teq[] = "1-z";
Dans ce fichier il y a les fonctions pour le fichier c17a.c

Mathc initiation/c36bb
Installer ce fichier dans votre répertoire de travail.
/* save as fb.h */
double f(
double x,
double y,
double z)
 return( (1) );
char feq[] = "1";
double ux(
double y,
double z)
 return(0);
char ueq[] = "0";
double vx(
double y,
double z)
 return((2-y-z));
char veq[] = "2-y-z";
double sz(
double y)
 return((-sqrt(1-y*y)));
char seq[] = "-sqrt(1-y**2)";
double tz(
double y)
 return((sqrt(1-y*y)));
char teq[] = "sqrt(1-y**2)";
Dans ce fichier il y a les fonctions pour le fichier c17b.c

Mathc initiation/c36bc
Installer ce fichier dans votre répertoire de travail.
/* save as fc.h */
double f(
double x,
double y,
double z)
 return( (1) );
char feq[] = "1";
double uy(
double x,
double z)
 return(0);
char ueq[] = "0";
double vy(
double x,
double z)
 return( (z*z+x*x+1) );
/* --------------------------------- */
char veq[] = "z**2+x**2+1";
double sx(
double z)
 return((0));
char seq[] = "0";
double tx(
double z)
 return( (2-2*z) );
/* --------------------------------- */
char teq[] = "2-2*z";

Mathc initiation/c36bd
Installer ce fichier dans votre répertoire de travail.
/* save as fd.h */
double f(
double x,
double y,
double z)
 return( (1) );
char feq[] = "1";
double uy(
double x,
double z)
 return(0);
char ueq[] = "0";
double vy(
double x,
double z)
 return((6-z));
char veq[] = "6-z";
double sz(
double x)
 return((3*x*x));
char seq[] = "3*x**2";
double tz(
double x)
 return((4-x*x));
char teq[] = "4-x**2";
Dans ce fichier il y a les fonctions pour le fichier c17d.c

Mathc initiation/c36be
Installer ce fichier dans votre répertoire de travail.
/* save as fe.h */
double f(
double x,
double y,
double z)
 return( (1) );
char feq[] = "1";
double uz(
double x,
double y)
 return(0);
char ueq[] = "0";
double vz(
double x,
double y)
 return((4-y));
char veq[] = "4-y";
double sx(
double y)
 return((-sqrt(y)));
char seq[] = "-sqrt(y)";
double tx(
double y)
 return((sqrt(y)));
char teq[] = "sqrt(y)";
Dans ce fichier il y a les fonctions pour le fichier c17e.c

Mathc initiation/c36bf
Installer ce fichier dans votre répertoire de travail.
/* save as ff.h */
double f(
double x,
double y,
double z)
 return(x*y);
char feq[] = "x*y";
/* --------------------------------- */
double uz(
double x,
double y)
 return(x*y);
char ueq[] = "x*y";
double vz(
double x,
double y)
 return((x*x*y*y*y));
char veq[] = "x**2*y**3";
/* --------------------------------- */
double sy(
double x)
 return((x));
char seq[] = "x";
double ty(
double x)
 return((x*x));
char teq[] = "x**2";
/* --------------------------------- */
Dans ce fichier il y a les fonctions pour le fichier c17f.c

Mathc initiation/c36ca
Installer et compiler ces fichiers dans votre répertoire de travail.
/* save as c17a.c */
int main(void)
double az,bz,nz, 
 nx, ny, M;
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);
 printf(" v : y,z -> %s \n", veq);
 printf(" u : y,z -> %s\n\n", ueq);
 printf(" t : z -> %s \n", teq);
 printf(" s : z -> %s\n\n", seq);
 az = 0;
 bz = 1;
 nz = 2*50; nx = 2*50; ny = 2*50;
 M = sympson_dxdydz( f,
 az, bz, nz,
 sy, ty, ny,
 ux, vx, nx);
 printf(" With the sympson's rule.\n\n\n");
 printf(" (%+.1f (%s (%s \n", bz, teq, veq);
 printf(" A = int( int( int( %s dx dy dz = %.6f\n",feq, M);
 printf(" (%+.1f (%s (%s\n\n\n", az, seq, ueq);
 stop();
 return 0;
/* --------------------------------- */
 Même résultat que la version avec l'algorithme des trapèzes à quelques décimales près.
Exemple de sortie écran :
 f : x,y,z -> y*z
 v : y,z -> 1-y-z 
 u : y,z -> 0
 t : z -> 1-z 
 s : z -> 0
 With the sympson's rule.
 (+1.0 (1-z (1-y-z 
 A = int( int( int( y*z dx dy dz = 0.008333
 (+0.0 (0 (0
 Press return to continue.

Mathc initiation/c36cb
Installer et compiler ces fichiers dans votre répertoire de travail.
/* save as c17b.c */
int main(void)
double ay,by,ny,
 nx,nz,M;
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);
 printf(" v : y,z -> %s \n", veq);
 printf(" u : y,z -> %s\n\n", ueq);
 printf(" t : y -> %s \n", teq);
 printf(" s : y -> %s\n\n", seq);
 ay = -1;
 by = 1;
 ny = 2*150; nx = 2*150; nz = 2*150;
 M = sympson_dxdzdy( f,
 ay,by,ny,
 sz,tz,nz,
 ux,vx,nx);
 printf(" With the sympson's rule.\n\n\n");
 printf(" (%+.1f (%s (%s \n", by, teq, veq);
 printf(" A = int( int( int( %s dx dz dy = %.6f\n",feq, M);
 printf(" (%+.1f (%s (%s\n\n\n", ay, seq, ueq);
 stop();
 return 0;
/* --------------------------------- */
 Même résultat que la version avec l'algorithme des trapèzes à quelques décimales près.
Exemple de sortie écran :
 f : x,y,z -> 1
 v : y,z -> 2-y-z 
 u : y,z -> 0
 t : y -> sqrt(1-y**2) 
 s : y -> -sqrt(1-y**2)
 With the sympson's rule.
 (+1.0 (sqrt(1-y**2) (2-y-z 
 A = int( int( int( 1 dx dz dy = 6.282685
 (-1.0 (-sqrt(1-y**2) (0
 Press return to continue.

Mathc initiation/c36cc
Installer et compiler ces fichiers dans votre répertoire de travail.
/* save as c17c.c */
int main(void)
double az,bz,nz,
 nx,ny,M;
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);
 printf(" v : x,z -> %s \n", veq);
 printf(" u : x,z -> %s\n\n", ueq);
 printf(" t : z -> %s \n", teq);
 printf(" s : z -> %s\n\n", seq);
 az = 0;
 bz = 1;
 nz = 2*50;nx = 2*50;ny = 2*50;
 M = sympson_dydxdz( f,
 az, bz, nz,
 sx, tx, nx,
 uy, vy, ny);
 printf(" With the sympson's rule.\n\n\n");
 printf(" (%+.1f (%s (%s \n", bz, teq, veq);
 printf(" A = int( int( int( %s dy dx dz = %.6f\n",feq, M);
 printf(" (%+.1f (%s (%s\n\n\n", az, seq, ueq);
 stop();
 return 0;
/* --------------------------------- */
 Même résultat que la version avec l'algorithme des trapèzes à quelques décimales près.
Exemple de sortie écran :
 f : x,y,z -> 1
 v : x,z -> z**2+x**2+1 
 u : x,z -> 0
 t : z -> 2-2*z 
 s : z -> 0
 With the sympson's rule.
 (+1.0 (2-2*z (z**2+x**2+1 
 A = int( int( int( 1 dy dx dz = 1.833333
 (+0.0 (0 (0
 Press return to continue.

Mathc initiation/c36cd
Installer et compiler ces fichiers dans votre répertoire de travail.
/* save as c17d.c */
int main(void)
double ax,bx,nx, 
 ny,nz,M;
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);
 printf(" v : x,z -> %s \n", veq);
 printf(" u : x,z -> %s\n\n", ueq);
 printf(" t : x -> %s \n", teq);
 printf(" s : x -> %s\n\n", seq);
 ax = -1;
 bx = 1;
 nx = 2*50;ny = 2*50;nz = 2*50;
 M = sympson_dydzdx( f,
 ax,bx,nx,
 sz,tz,nz,
 uy,vy,ny);
 printf(" With the sympson's rule.\n\n\n");
 printf(" (%+.1f (%s (%s \n", bx, teq, veq);
 printf(" A = int( int( int( %s dy dz dx = %.6f\n",feq, M);
 printf(" (%+.1f (%s (%s\n\n\n", ax, seq, ueq);
 stop();
 return 0;
/* --------------------------------- */
 Même résultat que la version avec l'algorithme des trapèzes à quelques décimales près.
Exemple de sortie écran :
 f : x,y,z -> 1
 v : x,z -> 6-z 
 u : x,z -> 0
 t : x -> 4-x**2 
 s : x -> 3*x**2
 With the sympson's rule.
 (+1.0 (4-x**2 (6-z 
 A = int( int( int( 1 dy dz dx = 20.266667
 (-1.0 (3*x**2 (0
 Press return to continue.