% Random number generator from Clocksin and Mellish (p. 149)

seed(13).

random(R, N) :-
	retract(seed(S)),
	N is (S mod R) + 1,
	NewSeed is (125 * S + 1) mod 4096,
	asserta(seed(NewSeed)), !.
