Fichier Gibiane

OPTI DIME 2 ELEM QUA4 ;
*
* Variables globales
*
* LX0 : Dimension suivant x
* LY0 : Dimension suivant y
* XMIN,YMIN : Coordonnées du coin inférieur
* NX0 : Nombre de points suivant x
* NY0 : Nombre de points suivant y
*
LX0  = 2.D0 ; LY0 = 1.D0 ;
XMIN = 0.D0 ; YMIN = 0.D0 ;
NX0  = 3 ; NY0 = 5 ;
XMAX = XMIN + LX0 ; YMAX = YMIN + LY0 ;
YMID = YMAX + YMIN / 2.D0 ;
*
* Création des points
P1 = XMIN YMIN ;
P2 = XMAX YMIN ;
P3 = XMAX YMID ; 
P4 = XMAX YMAX ;
P5 = XMIN YMAX ;
P6 = XMIN YMID ;
*
* Création des droites
P1P2 = P1 DROI NX0 P2 ;
P2P3 = P2 DROI NY0 P3 ;
P3P4 = P3 DROI NY0 P4 ;
P4P5 = P4 DROI NX0 P5 ;
P5P6 = P5 DROI NY0 P6 ;
P6P1 = P6 DROI NY0 P1 ;
P3P6 = P3 DROI NX0 P6 ;
P6P3 = INVE P3P6 ;
*
* Création des surfaces
SURF1 = DALL P1P2 P2P3 P3P6 P6P1 ;
'OPTION' 'ELEM' TRI3 ;
SURF2 = DALL P6P3 P3P4 P4P5 P5P6 ;
SURF0 = SURF1 ET SURF2 ;
CONT0 = CONT SURF0 ;
*
* Création d'un CHAMPOIN contenant la fonction z = ax + by + c
A1 = 2.D0 ; B1 = 3.D0 ; C0 = 5.D0 ;
X0 Y0 = COOR SURF0 ;
Z0 = A1*X0 + (B1*Y0) + C0 ;
*
* Evolution de z = ax + by + c le long de l'axe médian
ZP6P3 = EVOL CHPO Z0 P6P3 ;
*
* Tracé
*
TRAC SURF0 ;
TRAC Z0 SURF0 CONT0 ;
DESS ZP6P3 ;
*
* ARRET
FIN ;

Figure 1: Maillage

Figure 2: CHAMPOIN bilinéaire f(x,y)=2x+3y+5

Figure 3: Evolution du CHAMPOIN f(x,y)=2x+3y+5 le long de l'axe y=0.5