Mathc initiation/c36ca
Installer et compiler ces fichiers dans votre répertoire de travail.
c00a.c |
---|
/* --------------------------------- */
/* save as c00a.c */
/* --------------------------------- */
#include "x_hfile.h"
#include "fa.h"
/* --------------------------------- */
int main(void)
{
double M = simpson_dxdydz(f,
ux, vx, LOOP,
sy, ty, LOOP,
az, bz, LOOP);
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);
printf(" With the simpson'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);
printf(" Code Mathematica : \n\n"
" integral y z dx dy dz from 0 to "
"(1-y-z) from 0 to 1 from 0 to (1-z)\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 -> y*z
v : y,z -> 1-y-z
u : y,z -> 0
t : z -> 1-z
s : z -> 0
With the simpson'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.