Mathc initiation/Fichiers c : c62cc

Sommaire

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

	c00c.c

/* ---------------------------------- */
/* save as c00c.c                     */
/* ---------------------------------- */
#include "x_hfile.h"
#include      "fc.h"
/* ---------------------------------- */
int main(void)
{
double  m = IntFlux_simpson_dydz( M,N,P,
                                  k,
                                  u, v,LOOP,
                                 az,bz,LOOP);
 clrscrn();

 printf(" Let S be the part of the graph of x = %s.  \n\n", keq);
 printf(" If F(x,y,z) = %si %sj %sk, find the flux of F through S\n\n\n",
          Meq,Neq,Peq);
 printf(" Consider k(y,z) = (%s)\n\n",keq);
 printf("          n      = grad(k(y,z)) / ||grad(k(y,z))||\n\n\n");

 
                           
 printf(" The flux of F through S is \n\n");
 printf("     //                  \n");
 printf("    ||                   \n");
 printf("    || F.n  dS = %.3f\n",m);
 printf("    ||                   \n");
 printf("   //                    \n");
 printf("   S                 \n\n\n");
 stop();

 clrscrn();
 printf("     / b   / v(z)\n");
 printf("    |     |      \n");
 printf("    |     |     F.(+i-k_yj-k_zk)    [1+k_y**2+k_z**2]**1/2 dy dz = %.3f\n",m);
 printf("    |     |        -----------     \n");
 printf("    |     |      [1+k_y**2+k_z**2]**1/2\n");
 printf("    |     |      \n");
 printf("   /  a  /   u(z)\n\n\n");

 printf("     / b   / v(z)\n");
 printf("    |     |      \n");
 printf("    |     |     F.(+i-k_yj-k_zk)   dy dz = %.3f\n",m);
 printf("    |     |      \n");
 printf("   /  a  /   u(z)\n\n\n");
 stop();
 
 clrscrn();
 printf("     / b   / v(z)\n");
 printf("    |     |      \n");
 printf("    |     |     F.(+i-k_yj-k_zk)  dy dz = %.3f\n",m);
 printf("    |     |      \n");
 printf("   /  a  /   u(z)\n\n\n");

 printf(" With.\n\n\n");
 printf(" F : (x,y,z)-> %si %sj %sk \n\n", Meq,Neq,Peq);
 printf(" k :   (y,z)-> %s          \n\n", keq);

 printf(" v :     (z)-> %s          \n",  veq);
 printf(" u :     (z)-> %s        \n\n",  ueq); 
 printf(" b = %+.1f\n a = %+.1f \n\n\n", bz,az);
 stop();

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

Dans cette version nous utilisons l'algorithme simplifié qui calcule les intégrales de flux de surface.

     / b   / v(z)
    |     |      
    |     |     F.(+i-k_yj-k_zk)  dy dz = 4.000
    |     |      
   /  a  /   u(z)


 With.


 F : (x,y,z)-> x+yi + zj  + xzk 

 k :   (y,z)-> 1    

 u :     (z)-> -1      
 v :     (z)-> +1    

 a = -1.0   b = +1.0 

 Press return to continue.