clc;
clear;
close all;

%% vstupní parametry

rozsahX = [-5, 5;
        -5,15];
n = 2;
N = 10000;
varX = [1,1];
strX = [0,0];
rho = 0;

%%
covMat = [varX(1)^2, rho*varX(1)*varX(2);rho*varX(1)*varX(2),varX(2)^2];
matY=mvnrnd(strX,covMat,N)';

figure(1);
stem(matY(1,:),matY(2,:),'LineStyle','none');

%%
vektX = nan(n,N);
vektfX = nan(n,N);

covMat = [varX(1)^2, rho*varX(1)*varX(2);rho*varX(1)*varX(2),varX(2)^2];