import numpy as np #knihovna numpy
import matplotlib.pyplot as plt #knihovna matplot.lib.pyplot
n=50000
seed=48541
a=16807
m=2147483647
#a=65539
#m=2147483648

x=np.empty(n)           #pole nahodnych hodnot x
y=np.empty(n)           #pole nahodnych hodnot y
color=np.empty([n,3])   #pole nahodnych barev RGB

#generator vsech nahodnych promennych dohromady
old=seed
for i in range(0,n):
    #generator nahodnych promennych x
    new=a*old % m
    x[i]=new/m
    old=new
    #generator nahodnych promennych x
    new=a*old % m
    y[i]=new/m
    old=new
    #generator nahodnych barev
    new=a*old % m
    color[i,0]=new/m
    old=new
    new=a*old % m
    color[i,1]=new/m
    old=new
    new=a*old % m
    color[i,2]=new/m
    old=new

plt.scatter(x,y,s=0.1,c='black')   #nahodna cisla cernobile
plt.savefig("generator3-BW.png",dpi=300)
plt.show()

plt.scatter(x,y,s=1,c=color)    #nahodna cisla + nahodne barvy
plt.savefig("generator3-RGB.png",dpi=300)
plt.show()