INTEGER FUNCTION LCM(J,K)
C
C	This integer function returns the least common multiple (LCM)
C	of two positive integers (J & K).  Both J & K must not have
C	any prime factors greater than 241.
C
	DIMENSION PRIME(53)
	DATA PRIME/
	1  2., 3., 5., 7., 11., 13., 17.,
	2  19., 23., 29., 31., 37., 41., 43.,
	3  47., 53., 59., 61., 67., 71., 73.,
	4  79., 83., 89., 97., 101., 103., 107.,
	5  109., 113., 127., 131., 137., 139., 149.,
	6  151., 157., 163., 167., 173., 179., 181.,
	7  191., 193., 197., 199., 211., 223., 227.,
	8  229., 233., 239., 241./
C
	LCM = 1
	DO 1000 I=1,53
	DIVSOR = PRIME(I)
C
100	JREM = J - IFIX( DIVSOR*IFIX(J/DIVSOR) )
	KREM = K - IFIX( DIVSOR*IFIX(K/DIVSOR) )
	IF( JREM*KREM .NE. 0 ) GO TO 1000
	LCM = LCM*PRIME(I)
	DIVSOR = DIVSOR*PRIME(I)
	GO TO 100
C
1000	CONTINUE
	RETURN
	END