Esope source files

   1 : C FONCF     SOURCE    CHAT      05/01/13    00:06:02     5004           
   2 : C
   3 : C*********************************************************************
   4 : C
   5 :       FUNCTION FONCF(Q,X,Y,Z,X0)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 :       EPSIL2=1.D-8
   9 :       Q2=SQRT(Q*Q+1.)
  10 :       D=Z+Q*X0
  11 :       IF(D.EQ.0.) THEN
  12 :          Z=Z+EPSIL2
  13 :       ENDIF
  14 :       DELTA=SQRT((X-X0)**2+Y**2+(Q*X+Z)**2)
  15 : C****
  16 :       DELTA1=MAX(X*Q2+(Q*Z-X0)/Q2+DELTA,EPSIL2)
  17 :       DELTA2=MAX(Y+DELTA,EPSIL2)
  18 :       FONCF=Y*LOG(DELTA1)
  19 :      &     +(X*Q2+(Q*Z-X0)/Q2)*LOG(DELTA2)
  20 :      &     -(Z+Q*X0)/Q2*ATAN(Y*(X*(Q*Q+1.)+Q*Z-X0)/(Z+Q*X0)/DELTA)
  21 : C****
  22 :       RETURN
  23 :       END
  24 :  
  25 :