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);