Esope source files

   1 : C FPLQ4     SOURCE    CHAT      05/01/13    00:07:37     5004           
   2 :       FUNCTION FPLQ4(A,Y,Z)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 :       A2=A**2
   6 :       Y2=Y**2
   7 :       Z2=Z**2
   8 :       C1=SQRT(A2+Y2+Z2)
   9 :       APC1=A+C1
  10 :       C2=A2+Z2
  11 :       C3=SQRT(C2)
  12 :       IF(Z.NE.0.) THEN
  13 :         SINUS=(A*C1+C2)/(C3*APC1) 
  14 :       ELSE 
  15 :         SINUS=SIGN(1.D0,A)
  16 :       ENDIF
  17 :       FPLQ4=Y*LOG(APC1)-Y+A*LOG(Y+C1)
  18 :      &     -Z*ASIN(SINUS)
  19 :       RETURN
  20 :       END
  21 :  
  22 :