Télécharger foncf.eso

Retour à la liste

Numérotation des lignes :

  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.  
  26.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales