clc;
format compact; 
wig20_
n0=1800; 
x=wig20((n0+1):n0+200,4); 
y=wig20((n0+1):n0+200,5);
Kowariancja= cov(x);
Korelacja= corrcoef(x);
C=cov([x,y]);
b1=C(1,2)/C(1,1);
b0=mean(y)-b1*mean(x);
x_i=1:200;
y_1=1:200;
plot(x_i,y_i)
hold on
y_regr=b*x_i+b0;
plot(x_i,y_regr)
tytul=['Liczebnosc probki n= ', num2str(n), ' korelacja wynosi: ',num2str(korelacja)];
title(tytul)