Esope source files

   1 : C FZ2       SOURCE    CHAT      05/01/13    00:14:15     5004           
   2 :       FUNCTION FZ2(D,A,U,SQA,SQU,RSINT2)                                        
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 :       SQAD=SQA**2-D**2                                                          
   6 :       SQUD=SQU**2-D**2                                                          
   7 :       IF(SQAD.EQ.0.D0) GOTO 2                                                   
   8 : 1     IF(SQUD.EQ.0.D0) GOTO 3                                                   
   9 :       GOTO 4                                                                    
  10 : 2     FZ2=-SIGN(A,D)                                                            
  11 :       RETURN                                                                    
  12 : 3     FZ2=+SIGN(U,D)                                                            
  13 :       RETURN                                                                    
  14 : 4     FZ2=-D*RSINT2*(A/(SQAD*SQA)-U/(SQUD*SQU))                                 
  15 :       RETURN                                                                    
  16 :       END                                                                       
  17 :