Esope source files

   1 : C EQPL3P    SOURCE    CHAT      05/01/12    23:44:05     5004           
   2 :       FUNCTION EQPL3P(XX,X1,X2,X3)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C************************************************************************
   6 : C
   7 : C     Calcule la valeur de l'equation d'un plan passant par 3pts A,B,C
   8 : C     au point (x,y,z)
   9 : C************************************************************************
  10 :       DIMENSION XX(3),X1(3),X2(3),X3(3)
  11 : 
  12 :       XA=X1(1)
  13 :       YA=X1(2)
  14 :       ZA=X1(3)
  15 : 
  16 :       XB=X2(1)
  17 :       YB=X2(2)
  18 :       ZB=X2(3)
  19 : 
  20 :       XC=X3(1)
  21 :       YC=X3(2)
  22 :       ZC=X3(3)
  23 : 
  24 :       X=XX(1)
  25 :       Y=XX(2)
  26 :       Z=XX(3)
  27 : 
  28 :       FX=(YB-YA)*(ZC-ZA)-(ZB-ZA)*(YC-YA)
  29 :       FY=(XC-XA)*(ZB-ZA)-(XB-XA)*(ZC-ZA)
  30 :       FZ=(XB-XA)*(YC-YA)-(YB-YA)*(XC-XA)
  31 : 
  32 :       EQPL3P=FX*(X-XA)+FY*(Y-YA)+FZ*(Z-ZA)
  33 :       RETURN
  34 :       END
  35 :  
  36 :  
  37 :