function z = rndLCG(T)
  m    = 2^32;
  a    = 1664525;
  c    = 1013904223; 
  k    = 100;
  x(1) = 327680; 
  x(1) = time*10^5;

  for i = 2:T+k;
    x(i) = mod ( (a*x(i-1) + c), m );
  end
  z = x(k+1:T+k) / m;
  x = [];
end