Mathc initiation/001Q

Sommaire

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 = 5 * simpson_dzdydx(f,
                               uz, vz, LOOP,
                               sy, ty, LOOP,
                               ax, bx, LOOP);
                               
double  N =     simpson_dzdydx(g,
                               uz, vz, LOOP,
                               sy, ty, LOOP,
                               ax, bx, LOOP);
                               
 clrscrn();
 printf("First property of linearity                         \n\n" 
 
        "        (b (d (f                (b (d (f              \n"
        " (A) int(  f(x,y,z) dzdydx = int( (A) f(x,y,z) dzdydx \n"
        "        (a (c (e                (a (c (e          \n\n\n");       
 stop();
   
 clrscrn();
 printf(" I multiply int(f(x,y,z) dz dy dx by 5 \n\n");
 
 printf(" f : x,y,z ->   %s\n\n", feq);

 printf("             (%+.1f  (%s (%s    \n", bx, teq, veq);
 printf(" M = 5 * [int(   int(  int(    %s  dz dy dx] = %.6f\n",feq, M);
 printf("             (%+.1f  (%s    (%s\n\n\n", ax, seq, ueq);


 printf(" I multiply f(x,y,z) by 5 \n\n");
 
 printf(" g : x,y,z -> %s\n\n", geq);
 
 printf("             (%+.1f  (%s (%s    \n", bx, teq, veq);
 printf(" N =     [int(   int(  int( %s   dz dy dx  = %.6f\n",geq, N);
 printf("             (%+.1f  (%s    (%s\n\n\n", ax, seq, ueq);
 
  stop();

 return 0;
}
/* --------------------------------- */
/* --------------------------------- */

Verifier avec mathematica : Free: wolframalpha.com

integral 5 * (x**2 + y + z**2  +1)  dz dy dx from 0 to 1 from x to (x**2) from (x*y) to (x**2*y**3)

Exemple de sortie écran :

First property of linearity                         

        (b (d (f                (b (d (f              
 (A) int(  f(x,y,z) dzdydx = int( (A) f(x,y,z) dzdydx 
        (a (c (e                (a (c (e          


 Press return to continue. 


 I multiply int(f(x,y,z) dz dy dx by 5 

 f : x,y,z ->   (x**2+y+z**2+1)

             (+1.0  (x**2 (x**2*y**3    
 M = 5 * [int(   int(  int(    (x**2+y+z**2+1)  dz dy dx] = 0.295483
             (+0.0  (x    (x*y


 I multiply f(x,y,z) by 5 

 g : x,y,z -> 5*(x**2+y+z**2+1)

             (+1.0  (x**2 (x**2*y**3    
 N =     [int(   int(  int( 5*(x**2+y+z**2+1)   dz dy dx  = 0.295483
             (+0.0  (x    (x*y


 Press return to continue.