In [20]:
import numpy as np
import matplotlib.pyplot as plt

konstanty pro správný čistě multiplikativní generátor

In [22]:
a=16807
m=2147483647

semínko

In [24]:
i_seed=12345
In [25]:
ndata=100000
In [26]:
x = np.zeros(ndata, dtype=float) #deklarace pole x-souradnic
y = np.zeros(ndata, dtype=float) #deklarace pole y-souradnic
color=np.zeros([ndata,3],dtype=float) #deklarace pole barva RGB
In [27]:
i_old=i_seed
for i in range(0,ndata):
   i_next=(a*i_old) % m
   i_old=i_next
   x[i]=i_next/m
   i_next=(a*i_old) % m
   i_old=i_next
   y[i]=i_next/m
   i_next=(a*i_old) % m
   i_old=i_next
   color[i,0]=i_next/m
   i_next=(a*i_old) % m
   i_old=i_next
   color[i,1]=i_next/m
   i_next=(a*i_old) % m
   i_old=i_next
   color[i,2]=i_next/m
In [28]:
plt.scatter(x,y,c=color,s=1)
plt.xlim(0,1)
plt.ylim(0,1)
plt.show()
No description has been provided for this image
In [ ]: