Esope source files

   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 :