Mathc initiation/c36ca
Installer et compiler ces fichiers dans votre répertoire de travail.
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.