Sources Esope | Cast3M

foncf
C FONCF     SOURCE    CB215821  17/07/21    21:15:10     9513           
C
C*********************************************************************
C
      FUNCTION FONCF(Q,X,Y,Zin,X0)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      EPSIL2=1.D-8
      Z     = Zin
 
      Q2=SQRT(Q*Q+1.)
      D=Z+Q*X0
 
      IF(D.EQ.0.) THEN
         Z=Z+EPSIL2
      ENDIF
      DELTA=SQRT((X-X0)**2+Y**2+(Q*X+Z)**2)
C
      DELTA1=MAX(X*Q2+(Q*Z-X0)/Q2+DELTA,EPSIL2)
      DELTA2=MAX(Y+DELTA,EPSIL2)
      FONCF=Y*LOG(DELTA1)
     &     +(X*Q2+(Q*Z-X0)/Q2)*LOG(DELTA2)
     &     -(Z+Q*X0)/Q2*ATAN(Y*(X*(Q*Q+1.)+Q*Z-X0)/(Z+Q*X0)/DELTA)
C
      RETURN
      END