function rand_init(seed)

global seed1 seed2 seed3 seed4 seed5

seed1=mod(171*seed,30269);
seed2=mod(172*seed1,30307);
seed3=mod(170*seed2,30323);

k=fix(seed3/52774);
seed5 = 40692*(seed3-k*52774)-k*3791;
if(seed5<0)
    seed5=seed5+2147473399;
end
k=fix(seed5/53668);
seed4 = 40014*(seed5-k*53668)-k*12211;
if(seed4<0)
	seed4=seed4+2147483563;
end