MINICERN/mathlib/gen/u/u100pr.F.ori

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