Mathc initiation/c36cf
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.