C***********************************************************************
C   SUBROUTINE PLOTF OVERWRITES 'PLOT' WITH A PRINTER PLOT OF VARIOUS
C   THEORETICAL DISTRIBUTION FUNCTIONS.  THIS IS FOR COMPARISON TO AN
C   EMPIRICAL CDF
C***********************************************************************
      SUBROUTINE PLOTF (INDEX)
C***********************************************************************
C     DECLARIATIONS
C***********************************************************************
      INTEGER NORM01(100), EXP01(100), HN(100)
      CHARACTER*1  U, N, E, H
      CHARACTER*1 PLOT(51,68)
C***********************************************************************
C     COMMON STATEMENTS
C***********************************************************************
      COMMON /SPLOT/ PLOT
C***********************************************************************
C     DATA STATEMENTS
C***********************************************************************
      DATA U, N, E, H /'u','n','e','h'/
      DATA NORM01 /50,50,49,49,49,49,49,49,48,48,48,48,48,47,47,47,46,46
     1 ,46,45,45,44,44,44,43,43,42,42,41,40,40,39,39,38,37,37,36,35,35,3
     2 4,33,32,32,31,30,29,28,28,27,26,25,24,24,23,22,21,20,20,19,18,17,
     3 17,16,15,15,14,13,13,12,12,11,10,10,9,9,8,8,8,7,7,6,6,6,5,5,5,4,4
     4 ,4,4,4,3,3,3,3,3,3,2,2,2/
      DATA EXP01 /49,47,46,44,43,42,41,39,38,37,36,35,34,33,32,31,30,29,
     1 28,27,27,26,25,24,24,23,22,22,21,20,20,19,19,18,17,17,16,16,16,15
     2 ,15,14,14,13,13,13,12,12,11,11,11,11,10,10,10,9,9,9,9,8,8,8,8,7,7
     3 ,7,7,7,6,6,6,6,6,5,5,5,5,5,5,5,4,4,4,4,4,4,4,4,3,3,3,3,3,3,3,3,3,
     4 3,3,2/
      DATA HN /50,49,49,48,47,46,45,45,44,43,42,42,41,40,39,38,38,37,36,
     1 35,35,34,33,33,32,31,30,30,29,28,28,27,26,26,25,25,24,23,23,22,22
     2 ,21,20,20,19,19,18,18,17,17,16,16,15,15,15,14,14,13,13,13,12,12,1
     3 1,11,11,10,10,10,9,9,9,8,8,8,8,7,7,7,7,6,6,6,6,6,5,5,5,5,5,5,4,4,
     4 4,4,4,4,4,3,3,3/
C
      GO TO (10,30,50,70),INDEX
C   THIS SECTION OVERWRITES 'PLOT' WITH THE PLOT OF THE FITTED NORMAL
C   CDF.  IT IS INTENDED THAT SUBROUTINE 'PCDF' BE CALLED FIRST TO PLOT
C   SOME SAMPLE CDF.  ON THIS CALL TO 'PCDF' THE SCALING LIMITS MUST BE
C   PLTL SMEAN-2.*SD AND PLTU SMEAN+2.*SD WHERE SMEAN IS THE SAMPLE
C   MEAN AND SD IS THE SAMPLE STANDARD DEVIATION.
C   ARRAY NORM01 CONTAINS A GRAPH OF THE NORMAL(0,1) CDF IN THE SENSE
C   THAT PLOT(CDF(J),J) IS THE CORRECT POSITION FOR PLOTING F(J) FOR THE
C   SCALING OF THE X-AXIS OF (-2.,+2.)
C
 10   DO 20 J=1,100
      I=NORM01(J)
      PLOT(I,J*65/100)=N
 20   CONTINUE
      RETURN
C
C   THIS SECTION TAKES THE ARRAY PLOT AND OVERWRITES A PLOT OF THE FITTED
C   EXPONENTIAL CDF.IT IS INTENDED THAT THE SUBROUTINE 'PCDF' BE CALLED
C   FIRST. ON THIS CALL TO 'PCDF' THE SCALING LIMITS MUST BE PLTL O. AND
C   PLTU 3*SMEAN, WHERE SMEAN IS THE SAMPLE MEAN.
C
 30   DO 40 J=1,100
      I=EXP01(J)
      PLOT(I,J*65/100)=E
 40   CONTINUE
      RETURN
C***********************************************************************
C     THIS SECTION OVERWRITES THE PLOT WITH THE UNIFORM DISTRIBUTION
C     CDF.   THE VALUES OF PLTL AND PLTU ARE THE SMALLEST AND LARGEST
C     MEASURMENTS OR THE THEORETICAL LOWER AND UPPER BOUNDS FOR THE
C     DISTRIBUTION.
C***********************************************************************
 50   DO 60 J=2,98,2
      I=51-J/2
 60   PLOT(I,J*65/100)=U
      RETURN
C***********************************************************************
C     THIS SECTION OVERWRITES THE PLOT WITH THE HALF-NORMAL DISTRIBUTION
C     CDF.  THE VALUES OF PLTL AND PLTU SHOULD BE 0. AND 2.*SD WHERE SD
C     IS THE SAMPLE STANDARD DEVIATION.
C***********************************************************************
 70   DO 80 J=1,100
      I=HN(J)
 80   PLOT(I,J*65/100)=H
      RETURN
      END