Installer et compiler ces fichiers dans votre répertoire de travail.
![]() |
c18c.c ' |
---|
/* ---------------------------------- */
/* save as c18c.c */
/* ---------------------------------- */
#include "x_hfile.h"
#include "fc.h"
/* ---------------------------------- */
int main(void)
{
double az = -1;
double bz = 1.;
int nz = 2*100;
int ny = 2*100;
double m= 0;
/* ---------------------------------- */
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 g(x,y,z) = x - (%s)\n\n",keq);
printf(" n = grad(g(x,y,z)) / ||grad(g(x,y,z))||\n\n\n");
m = flux_dydz( M,N,P,
k, u,v,ny,
az,bz,nz,
H);
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(" 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(" u : (z)-> %s \n", ueq);
printf(" v : (z)-> %s \n\n", veq);
printf(" a = %+.1f b = %+.1f \n",az,bz);
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(" u : (z)-> %s \n", ueq);
printf(" v : (z)-> %s \n\n", veq);
printf(" a = %+.1f b = %+.1f \n",az,bz);
stop();
return 0;
}
/* ---------------------------------- */
/* ---------------------------------- */
L'algorithme consiste à adapter la fonction qui calcule les intégrales doubles au calcul des flux.
Remarque :
Dans cette version nous utilisons cet algorithme simplifié pour l'intégrale.
/ b / v(z)
| |
| | F.(+i-k_yj-k_zk) dy dz =
| |
/ a / u(z)
Dans la précédente version nous avons utilisé cet algorithme.
/ b / v(z)
| |
| | F.(+i-k_yj-k_zk) [1+k_y^2+k_z^2]^1/2 dy dz =
| | -----------
| | [1+k_y^2+k_z^2]^1/2
| |
/ a / u(z)
Exemple de sortie écran :
Let S be the part of the graph of x = 1.
If F(x,y,z) = x+yi + zj + xzk, find the flux of F through S
Consider g(x,y,z) = x - (1)
n = grad(g(x,y,z)) / ||grad(g(x,y,z))||
The flux of F through S is
//
||
|| F.n dS = 4.000
||
//
S
Press return to continue.
Exemple de sortie écran :
/ 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.
Cet article est issu de Wikibooks. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.