Télécharger pnormj.eso

Retour à la liste

Numérotation des lignes :

  1. C PNORMJ SOURCE CHAT 05/01/13 02:15:32 5004
  2. FUNCTION PNORMJ (N,ALPHA,BETA)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. ONE = 1.D0
  6. TWO = 2.D0
  7. DN = DBLE(N)
  8. CONST = ALPHA+BETA+ONE
  9. IF (N.LE.1) THEN
  10. PROD = GAMMAF(DN+ALPHA)*GAMMAF(DN+BETA)
  11. PROD = PROD/(GAMMAF(DN)*GAMMAF(DN+ALPHA+BETA))
  12. PNORMJ = PROD * TWO**CONST/(TWO*DN+CONST)
  13. RETURN
  14. ENDIF
  15. PROD = GAMMAF(ALPHA+ONE)*GAMMAF(BETA+ONE)
  16. PROD = PROD/(TWO*(ONE+CONST)*GAMMAF(CONST+ONE))
  17. PROD = PROD*(ONE+ALPHA)*(TWO+ALPHA)
  18. PROD = PROD*(ONE+BETA)*(TWO+BETA)
  19. DO 100 I=3,N
  20. DINDX = DBLE(I)
  21. FRAC = (DINDX+ALPHA)*(DINDX+BETA)/(DINDX*(DINDX+ALPHA+BETA))
  22. PROD = PROD*FRAC
  23. 100 CONTINUE
  24. PNORMJ = PROD * TWO**CONST/(TWO*DN+CONST)
  25. RETURN
  26. END
  27.  
  28.  
  29.  

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