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 ax = -2.;
double bx = 2.;
int nx = 2*50;
int ny = 2*50;
int nz = 2*50;
double m = 0;
/* --------------------------------- */
clrscrn();
printf(" Use the divergence theorem to find,\n\n");
printf(" the flux of F through S.\n\n");
printf(" // /// \n");
printf(" || ||| \n");
printf(" || F.n dS = ||| div F dV \n");
printf(" || ||| \n");
printf(" // /// \n");
printf(" S Q \n\n\n");
printf(" If F = Mi + Nj + Pk \n\n\n");
printf(" /// /// \n");
printf(" ||| ||| \n");
printf(" ||| div F dV = ||| M_x + N_y + P_z dV \n");
printf(" ||| ||| \n");
printf(" /// /// \n");
printf(" Q Q \n\n\n");
stop();
/* --------------------------------- */
clrscrn();
m = flux_dydzdx(
M,N,P,
ax,bx,nx,
s, t, nz,
u,v,ny,
H);
printf(" / b / t(x) / v(x, z) \n");
printf(" | | | \n");
printf(" | | | M_x N_y P_z dydzdx = %.3f \n",m);
printf(" | | | \n");
printf(" / a / s(x) / u(x, z) \n\n\n");
printf(" With.\n\n\n");
printf(" F : (x,y,z)-> %si %sj %sk \n\n",Meq,Neq,Peq);
printf(" u : (x,z)-> %s \n", ueq);
printf(" v : (x,z)-> %s \n\n", veq);
printf(" s : (x)-> %s \n", seq);
printf(" t : (x)-> %s \n\n", teq);
printf(" ax = %+.1f bx = %+.1f\n\n",ax,bx);
printf(" dV = dy dz dx \n\n");
stop();
return 0;
}
/* --------------------------------- */
/* --------------------------------- */
Ce travail consiste à adapter l'intégrale triple au calcul du flux en 3d par le théorème de la divergence : (M_x + N_y + P_z)
Exemple de sortie écran :
Use the divergence theorem to find,
the flux of F through S.
// ///
|| |||
|| F.n dS = ||| div F dV
|| |||
// ///
S Q
If F = Mi + Nj + Pk
/// ///
||| |||
||| div F dV = ||| M_x + N_y + P_z dV
||| |||
/// ///
Q Q
Press return to continue.
Exemple de sortie écran :
/ b / t(x) / v(x, z)
| | |
| | | M_x N_y P_z dydzdx = 131.657
| | |
/ a / s(x) / u(x, z)
With.
F : (x,y,z)-> + x^3 + sin(z)i + x^2*y + cos(z)j + exp(x^2 + y^2)k
u : (x,z)-> +0
v : (x,z)-> 5-z
s : (x)-> +0
t : (x)-> 4-x^2
ax = -2.0 bx = +2.0
dV = dy dz dx
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.