C-----------------------------------------------------------------------
C IMSL Name:  PRIME (Single precision version)
C
C Purpose:    Decompose an integer into its prime factors.
C
C Usage:      CALL PRIME (N, NPF, IPF, IEXP, IPW)
C
C Example 1:
      INTEGER    N
      PARAMETER  (N=144)
C
      INTEGER    IEXP(13), IPF(13), IPW(13), NOUT, NPF
      EXTERNAL   PRIME, UMACH
C                                  Get prime factors of 144
      CALL PRIME (N, NPF, IPF, IEXP, IPW)
C                                  Get output unit number
      CALL UMACH (2, NOUT)
C                                  Print results
      WRITE (NOUT,99999) N, IPF(1), IPF(2), IEXP(1), IEXP(2), IPW(1),
     &                  IPW(2), NPF
C
99999 FORMAT ('  The prime factors for', I5, ' are: ', /, 10X, 2I6, //
     &       , '  IEXP =', 2I6, /, '  IPW  =', 2I6, /, '  NPF  =', I6,
     &       /)
      END
C   The prime factors for  144 are:
C                2     3
C
C   IEXP =     4     2
C   IPW  =    16     9
C   NPF  =     2
C