*
* $Id: u100pr.F,v 1.1.1.1 1996/04/01 15:01:46 mclareni Exp $
*
* $Log: u100pr.F,v $
* Revision 1.1.1.1  1996/04/01 15:01:46  mclareni
* Mathlib gen
*
*
#include "gen/pilot.h"
      SUBROUTINE U100PR (I,K,M)
      DIMENSION IPRIME(40),K(40)
      DATA IPRIME/2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,
     + 71,73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,
     + 157,163,167,173/
      II=I
      M=1
      IF(I.LT.1)  GO TO 6
      CALL VZERO (K, 40)
      IF(I.EQ.1)  RETURN
      DO 2 KL=1,40
    1 IF(MOD(I,IPRIME(KL)).NE.0)   GO TO 2
      K(KL)=K(KL)+1
      I=I/IPRIME(KL)
      IF(I.EQ.1)  GO TO 3
      GO TO 1
    2 CONTINUE
      IF(I.NE.1)    WRITE(6,4)    I
    4 FORMAT(' PRIME NO. LARGER THAN THE LIST,THE FACTOR LEFT=',I15)
      M=I
    3 I=II
      RETURN
    6 WRITE(6,7)
    7 FORMAT(' PLEASE OBSERVE THE CONVENTION OF YOUR ASKING NO. BEING',
     F'LARGER THAN +0')
      RETURN
      END