Mathc initiation/c36ca


Sommaire


Installer et compiler ces fichiers dans votre répertoire de travail.

Crystal Clear mimetype source c.png c17a.c
'
/* ---------------------------------- */
/* save as c17a.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "fa.h"
/* --------------------------------- */
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.