Mathc initiation/c35cd

Sommaire

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

	c00d.c

/* ---------------------------------- */
/* save as c00d.c                    */
/* --------------------------------- */
#include  "x_hfile.h"
#include       "fd.h"
/* --------------------------------- */
int main(void)
{
double  ax = -1,
        bx =  1,
        nx = 2*50,   ny = 2*50,  nz = 2*50;

double  M =  trapezoid_dydzdx(f,
                              uy,vy,ny,
                              sz,tz,nz,
                              ax,bx,nx);
 
 clrscrn();
 printf(" f : x,y,z -> %s\n\n", feq);
 
 printf(" v : x,z   -> %s  \n", veq);
 printf(" u : x,z   -> %s\n\n", ueq);

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

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

 printf(" Code Mathematica : \n\n"
        " integral 1 dy dz dx from (0) to"
        " (6-z) from -1 to 1 from (3*x**2) to (4-x**2)\n\n\n");
 stop();
 
 return 0;
}
/* --------------------------------- */

Exemple de sortie écran :

 f : x,y,z -> 1

 v : x,z   -> 6-z  
 u : x,z   -> 0

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

 With the trapezoidal's rule.


        (+1.0  (4-x**2     (6-z    
 A = int(   int(    int(      1  dy dz dx = 20.265067
        (-1.0  (3*x**2     (0


 Press return to continue.