UNISF
      SUBROUTINE UNISF(P,SF)
  
PURPOSE--This subroutine computes the sparsity function value
         for the uniform (rectangular) distribution on the unit
         interval (0,1).  This distribution has MEAN = 0.5 and
         STANDARD DEVIATION = SQRT(1/12) = 0.28867513.  This
         distribution has the probability density function
         F(X) = 1.  Note that the sparsity function of a
         distribution is the derivative of the percent point
         function, and also is the reciprocal of the
         probability function density (but in units of P rather
         than X).
INPUT  ARGUMENTS--P      = The single precision value (between
                           0.0 and 1.0) at which the sparsity
                           function is to be evaluated.
OUTPUT ARGUMENTS--SF     = The single precision sparsity
                           function value.
OUTPUT--The single precision sparsity function value SF>
PRINTING--None unless an input argument error condition exists.
RESTRICTIONS--P should be between 0.0 and 1.0, inclusively.
OTHER DATAPAC   SUBROUTINES NEEDED--None.
FORTRAN LIBRARY SUBROUTINES NEEDED--None.
MODE OF INTERNAL OPERATIONS--Single precision.
LANGUAGE--ANSI FORTRAN.
REFERENCES--Filliben, Simple and Robust Linear Estimation
            of the Location Parameter of a Symmetric
            Distribution (unpublished PH.D. Dissertation,
            Princeton University), 1969, Pages 21-44, 229-231.
          --Filliben, 'The Percent Point Function',
            (unpublished manuscript), 1970, Pages 28-31.
          --Johnson and Kotz, Continuous Univariate
            Distributions--2, 1970, Pages 57-74.
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.