Télécharger foncg.eso

Retour à la liste

Numérotation des lignes :

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

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