Télécharger intfr2.eso

Retour à la liste

Numérotation des lignes :

  1. C INTFR2 SOURCE CB215821 17/11/30 21:16:29 9639
  2. SUBROUTINE INTFR2(CC,SS,XF,XD)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. C =====================================================================
  7. C SOUS-PROGRAMME FORTRAN APPELE PAR ELFGR1 POUR CALCULER LES INTEGRALES
  8. C
  9. C / XF
  10. C CC=1/SQRT(2PI)/ COS(T)/SQRT(T) DT
  11. C / XD
  12. C
  13. C / XF
  14. C SS=1/SQRT(2PI)/ SIN(T)/SQRT(T) DT
  15. C / XD
  16. C
  17. C AVEC 1 < XF ET XD < 10
  18. C
  19. C PROGRAMMEUR : GUILBAUD
  20. C CREATION : 13/01/87
  21. C =====================================================================
  22. C
  23. -INC CCREEL
  24. DEUPI=2.D0*XPI
  25. DX2=XPI/2000.D0
  26. C
  27. C INTEGRATION PAR LA METHODE DES TRAPEZES
  28. C
  29. KMAX=INT((XF-XD)/DX2)
  30. IF(KMAX.EQ.0) KMAX=1
  31. DX=(XF-XD)/(KMAX+1)
  32. SSX=SIN(XD)*0.5D0/SQRT(XD)
  33. SCX=COS(XD)*0.5D0/SQRT(XD)
  34. XX=XD
  35. DO 20 K=1,KMAX
  36. XX=XX+DX
  37. SSX=SSX+SIN(XX)/SQRT(XX)
  38. SCX=SCX+COS(XX)/SQRT(XX)
  39. 20 CONTINUE
  40. SS=(SSX+0.5D0*SIN(XF)/SQRT(XF))*DX/SQRT(DEUPI)
  41. CC=(SCX+0.5D0*COS(XF)/SQRT(XF))*DX/SQRT(DEUPI)
  42. RETURN
  43. END
  44.  
  45.  
  46.  
  47.  
  48.  

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