function seed_rand()
% Initialize random number generator to a fixed seed
% This ensures that results are reproducible

% Try to make different versions of matlab happy
try
  RandStream.setDefaultStream(RandStream('mt19937ar','seed',3));
catch
  rand('twister',3);
end