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()
In [ ]: