MAX
      SUBROUTINE MAX(X,N,IWRITE,XMAX)
  
PURPOSE--This subroutine computes the sample maximum of the data
         in the input vector X.
INPUT  ARGUMENTS--X      = The single precision vector of
                           (unsorted or sorted) observations.
                --N      = The integer number of observations
                           in the vector X.
                --IWRITE = An integer flag code which (if set
                           to 0) will suppress the printing of
                           the sample maximum as it is computed;
                           or (if set to some integer value not
                           equal to 0, like, say, 1) will cause
                           the printing of the sample maximum
                           at the time it is computed.
OUTPUT ARGUMENTS--XMAX   = The single precision value of the
                           computed sample maximum.
OUTPUT--The computed single precision value of the sample
        maximum.
PRINTING--None, unless IWRITE has been set to a non-zero
          integer, or unless an input argument error
          condition exists.
RESTRICTIONS--There is no restriction on the maximum value of N
              for this subroutine.
OTHER DATAPAC   SUBROUTINES NEEDED--None.
FORTRAN LIBRARY SUBROUTINES NEEDED--None.
MODE OF INTERNAL OPERATIONS--Single precision.
LANGUAGE--ANSI FORTRAN.
REFERENCES--David, Order Statistics, 1970, Page 7.
WRITTEN BY--James J. Filliben
            Statistical Engineering Laboratory (205.03)
            National Bureau of Standards
            Gaithersburg, MD  20899
            Phone--  301-921-2315
ORIGINAL VERSION--June      1972.
UPDATED         --September 1975.
UPDATED         --November  1975.