*-- Author :    Stephen Burke   05/05/95
      SUBROUTINE FPPBIN
C----------------------------------------------------                                                         
C                                                                                                             
C---  Initialize an array parray with probabilities                                                           
C---  for chisq versus number of degrees of freedom                                                           
C---  for ndf .lt. 20                                                                                         
C---  first bin n    1-1000 is chisq in n/16ths                                                               
C---  first bin n 1001-2000 is chisq in n/8ths                                                                
C---  for ndf .ge. 20 .and. ndf .le.50                                                                        
C---  first bin n    1-1000 is chisq in n/8ths                                                                
C---  first bin n 1001-2000 is chisq in n/4ths                                                                
C---  second bin is NDF 6-55                                                                                  
C---  if greater than 55 the prob function used directly                                                      
C                                                                                                             
C----------------------------------------------------                                                         
      COMMON/FPBCOM/PARRAY(1800,50)                                     
                                                                        
      DO 300 J=2,51                                                     
         IF (J.LT.20) THEN                                              
            NFACT1 = 16                                                 
            NFACT2 = 8                                                  
         ELSE                                                           
            NFACT1 = 8                                                  
            NFACT2 = 4                                                  
         ENDIF                                                          
         DO 100 I=1,1000                                                
            PARRAY(I,J-1) = PROB(FLOAT(I)/NFACT1,J)                     
 100     CONTINUE                                                       
         DO 200 I=501,1300                                              
            PARRAY(I+500,J-1) = PROB(FLOAT(I)/NFACT2,J)                 
 200     CONTINUE                                                       
 300  CONTINUE                                                          
                                                                        
      RETURN                                                            
      END                                                               
*