function [y]=congruential1(x)

b=0;
m=8192;
x0=1;