x=normrnd(0,1,300,1);
y=normrnd(0,1,300,1);
subplot(2,2,1);plot(y,x,'.'); 
set(gca,'xtick',[],'ytick',[],'xlim',[-3 3],'ylim',[-3 3]);xlabel('x');ylabel('y');
Korelacja=corrcoef(x,y); %macierz korelacji
text(-2,3.1,sprintf('\\rho_{x, y} ~ 0  (%.2f)', Korelacja(1,2)));box off;
y=.7*x+normrnd(0,1,300,1);
subplot(2,2,2);plot(y,x,'.');
set(gca,'xtick',[],'ytick',[],'xlim',[-3 3],'ylim',[-3 3]);xlabel('x');ylabel('y');
Korelacja=corrcoef(x,y); %macierz korelacji
text(-2,3.1,sprintf('\\rho_{x, y} > 0  (%.2f)', Korelacja(1,2)));box off
y=x;
subplot(2,2,3);plot(y,x,'.');
set(gca,'xtick',[],'ytick',[],'xlim',[-3 3],'ylim',[-3 3]);xlabel('x');ylabel('y');
Korelacja=corrcoef(x,y); %macierz korelacji
text(-2,3.1,sprintf('\\rho_{x, y} = 1  (%.2f)', Korelacja(1,2)));box off;
y=-x+0.3*normrnd(0,1,300,1);
subplot(2,2,4);plot(y,x,'.');
set(gca,'xtick',[],'ytick',[],'xlim',[-3 3],'ylim',[-3 3]);xlabel('x');ylabel('y');
Korelacja=corrcoef(x,y); %macierz korelacji
text(-2,3.1,sprintf('\\rho_{x, y} << 0  (%.2f)', Korelacja(1,2)));box off;
