Installer et compiler ces fichiers dans votre répertoire de travail. .
 |
c18e.c' |
|---|
/* ---------------------------------- */
/* save as c18e.c */
/* ---------------------------------- */
#include "x_hfile.h"
#include "fe.h"
/* ---------------------------------- */
int main(void)
{
double ax = 0;
double bx = 2.;
int nx = 2*50;
int nz = 2*50;
double m = 0;
clrscrn();
m = surface_dzdx( g,h,
u,v,
nz,
ax,bx,nx,
H);
printf(" Evaluate this surface integrals equation.\n\n");
printf(" // \n");
printf(" || \n");
printf(" || g(x,h(x,z),z)* sqrt[ h_x(x,z)^2 + h_y(x,z)^2 + 1 ] dA = %.3f\n",m);
printf(" || \n");
printf(" // \n");
printf(" Rzx \n\n\n");
printf(" With.\n\n\n");
printf(" g : (x,y,z)-> %s\n\n", geq);
printf(" h : (x,z)-> %s\n\n", feq);
printf(" u : (x)-> %s \n", ueq);
printf(" v : (x)-> %s\n\n", veq);
printf(" a = %+.1f b = %+.1f \n",ax,bx);
stop();
return 0;
}
/* ---------------------------------- */
/* ---------------------------------- */
L'algorithme consiste à adapter la fonction qui calcule les intégrales doubles au calcul des surfaces en 3d.
Exemple de sortie écran :
Evaluate this surface integrals equation.
//
||
|| g(x,h(x,z),z)* sqrt[ h_x(x,z)^2 + h_y(x,z)^2 + 1 ] dA = -181.417
||
//
Rzx
With.
g : (x,y,z)-> x**2-2*y+z
h : (x,z)-> 8-4*x
u : (x)-> 0
v : (x)-> 6
a = +0.0 b = +2.0
Press return to continue.
Sortie écran d'un CAS :
> restart:
> g := (x,y,z)-> x**2-2*y+z:
> h := (x,z)-> 8-4*x:
> Int(
> Int(
> g(x,h(x,z),z)*
> sqrt( diff(h(x,z),x)^2 + diff(h(x,z),z)^2 + 1),
>
> z = 0 .. 6),
> x = 0 .. 2);
2 6
/ /
| | 2
| | (x - 16 + 8 x + z) sqrt(17) dz dx
| |
/ /
0 0
> evalf(value(%));
>
-181.4166475
>
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.