I want to use the mersenne twister to generate 'N' random numbers between 10 to 50. I want to be able to generate the same sequence over and over again.
I wrote the following code: (seed = 50, a = 10, b = 50, N = number of required random numbers)
s = rng(seed, 'twister');
r = a + (b-a)*rand(N,1);
rng(s);
r1 = a + (b-a)*rand(N,1);
Now even I print
r1 - r
I don't get zero. I expect to get zero as I have reset the random number generator to it's initial state in the third line of my code.
My question is where am I going wrong?
Aucun commentaire:
Enregistrer un commentaire