Esope source files

   1 : C SURF3P    SOURCE    CHAT      05/01/13    03:27:49     5004           
   2 :       SUBROUTINE SURF3P(P1,P2,P3,S)
   3 : C=============================================================
   4 : C                                                            =
   5 : C     INPUT:                                                 =
   6 : C     P1,P2 ET P3 SONT TROIS POINTS DANS UN PLAN             =
   7 : C                                                            =
   8 : C     OUTPUT:                                                =
   9 : C      S : SURFACE DU TRIANGLE                               =
  10 : C                                                            =
  11 : C     G. L. JIANG                                            =
  12 : C                                                            =
  13 : C=============================================================
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)
  16 :       DIMENSION P1(*),P2(*),P3(*)
  17 :       P=P2(1)-P1(1)
  18 :       Q=P2(2)-P1(2)
  19 :       R=P2(3)-P1(3)
  20 :       A=P3(1)-P1(1)
  21 :       B=P3(2)-P1(2)
  22 :       C=P3(3)-P1(3)
  23 :       C1=C*Q-B*R
  24 :       C2=A*R-C*P
  25 :       C3=B*P-A*Q
  26 :       S=0.5D0*SQRT(C1*C1+C2*C2+C3*C3)
  27 :       END
  28 :  
  29 :  
  30 :