% Q1
V1 = randn(1000,1);
mean(V1)
var(V1)
std(V1)

% Q2
M1 = 3 + sqrt(5)*randn(10,10);
mean(M1(:))  % use M1(:) to extract all the elements of the matrix M1 and display as a vector
var(M1(:))
cov(M1)
corr(M1)

% Q3
M2 = 5*randn(10,8);
corr(M1,M2)

% Q4
plot(V1)
hist(V1,20)
hist(M2(:),20)

