Mathc initiation/003t

Sommaire

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

	c00b2.c

/* --------------------------------- */
/* save as c00b2.c                   */
/* --------------------------------- */
#include  "x_afile.h"
#include      "fb2.h"
/* --------------------------------- */
int main(void)
{
double M =  J_wvu(             f,
                           X,  Y,  Z,
                           r0, r1, LOOP,
                           b0, b1, LOOP,
                           a0, a1, LOOP); 
 clrscrn();
 
 printf("  f :  x,y,z -> %s\n\n", feq);
 
 printf(" X  :  r,k,z -> %s       \tThe change of variables.\n",Xeq);
 printf(" Y  :  r,k,z -> %s                                 \n",Yeq);
 printf(" Z  :  r,k,z -> %s                               \n\n",Zeq);

 printf(" r1 :    r,k -> %s \t\t\t\tThe spherical coordinates.\n",r1eq);
 printf(" r0 :    r,k -> %s                                 \n\n",r0eq);

 printf(" b1 :    r   -> %s  \n", b1eq);
 printf(" b0 :    r   -> %s\n\n", b0eq);

 printf(" a1 :        -> %s  \n", a1eq); 
 printf(" a0 :        -> %s\n\n", a0eq);
                          
 printf("        (%s (%s  (%s      \n",               a1eq,b1eq,r1eq);
 printf(" A = int(  int(   int( f(X,Y,Z) |J| drdbda = %.6f\n", M);
 printf("        (%s    (%s     (%s\n\n",               a0eq,b0eq,r0eq);

 printf(" Code Mathematica : \n"
        " integral (r*cos(b))*(r**2 *sin(b)) dr db da"
        " from 0 to (2*Pi) from 0 to Pi/2 from 0 to (1)\n\n");                           
 stop();

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

a = xy;... ... b = zr;... ... r = r

Exemple de sortie écran :

  f :  x,y,z -> (z)

 X  :  r,k,z -> (r*sin(b)*cos(a))       	The change of variables.
 Y  :  r,k,z -> (r*sin(b)*sin(a))                                 
 Z  :  r,k,z -> (r*cos(b))                               

 r1 :    r,k -> 1 				The spherical coordinates.
 r0 :    r,k -> 0                                 

 b1 :    r   -> PI/2  
 b0 :    r   -> 0

 a1 :        -> 2*Pi  
 a0 :        -> 0

        (2*Pi (PI/2  (1      
 A = int(  int(   int( f(X,Y,Z) |J| drdbda = 0.785398
        (0    (0     (0

 Code Mathematica : 
 integral (r*cos(b))*(r**2 *sin(b)) dr db da from 0 to (2*Pi) from 0 to Pi/2 from 0 to (1)

 Press return to continue.