FUNCTION TEST0(X,Y)
C  CETTE FONCTION DECRIT UN CUBE  DE COTE 10
 X1=0
 X2=10
 Y1=0
 Y2=10
 H=10
 PX=(X-X1)*(X-X2)
 PY=(Y-Y1)*(Y-Y2)
 IF(PX.GT.0.) TX=0
 IF(PX.LT.0.) TX=1
 IF(PY.GT.0.) TY=0
 IF(PY.LT.0.) TY=1
 A=TX*TY
 IF(A.EQ.0) T=0
 IF(A.EQ.1) T=H
 TEST0=T
 END
CCCCCCCCCCCCCCCCCCCCCCCC
 FUNCTION HA(X,Y,X1,X2,Y1,Y2,H)
C LA FCT VAUT H DS LE RECTANGLE X1,X2,Y1,Y2 
C ET 0 AILLEURS
C POUR AVOIR H IL FAUT DONC QUE:
C X1<X<X2 ET Y1<Y<Y2
 PX=(X-X1)*(X-X2)
 PY=(Y-Y1)*(Y-Y2)
 IF(PX.GT.0.) TX=0
 IF(PX.LT.0.) TX=1
 IF(PY.GT.0.) TY=0
 IF(PY.LT.0.) TY=1
 A=TX*TY
 IF(A.EQ.0) HA=0
 IF(A.EQ.1) HA=H
C WRITE(6,100) X,Y,X1,X2,Y1,Y2,HA
100  FORMAT(' X,Y=',2F4.1,' X12,Y12=',4F4.1,' FA=',F3.0)
 END