Sommaire


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

c00f.c
/* --------------------------------- */
/* save as c00f.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "ff.h"
/* --------------------------------- */
int main(void)
{
double  M =  simpson_dzdydx(f,
                            uz, vz, LOOP,
                            sy, ty, LOOP,
                            ax, bx, LOOP);
 
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);

 printf(" v : x,y   -> %s  \n", veq); 
 printf(" u : x,y   -> %s\n\n", ueq);

 printf(" t : x     -> %s  \n", teq); 
 printf(" s : x     -> %s\n\n", seq);

 printf(" With the simpson's rule.\n\n\n");
 printf("        (%+.1f  (%s   (%s    \n", bx, teq, veq);
 printf(" A = int(   int(    int(  %s  dz dy dx = %.6f\n",feq, M);
 printf("        (%+.1f  (%s      (%s\n\n\n", ax, seq, ueq);

 printf(" Code Mathematica : \n\n"
        " integral x y   dz dy dx from (x*y) to "
        "(x**2*y**3) from 0 to 1 from x to (x**2)\n\n\n");
 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 -> x*y

 v : x,y   -> x**2*y**3  
 u : x,y   -> x*y

 t : x     -> x**2  
 s : x     -> x

 With the simpson's rule.


        (+1.0  (x**2     (x**2*y**3    
 A = int(   int(    int(      x*y  dz dy dx = 0.010582
        (+0.0  (x     (x*y


 Press return to continue.