Z Skrypty dla studentów Ekonofizyki UPGOW

Generacja liczb losowych

\(\int_0^1 \displaystyle\frac{1}{(1+x^a)dx}\)

function y=myran(x); 
  a=1664525;
  b=1013904223;
  m=2^32;
  y=mod(a*x+b,m);
  return; 
end

x(1)=123;
for i=2:10; 
	x(i)=myran(x(i-1));
	disp(x(i)/2^32);
end

Liczby o zadanym rozkładzie

Rozkład Gaussa

Algorytm Boxa-Mullera