%% Task 3.0
clear all
load deciles

%% Task 3.1
help std
doc std

%% Task 3.2
sort(dec10)

%% Task 3.2.1
[dec10_sorted,dec10_indices]=sort(dec10)

%% Task 3.2.2
[dec01_sorted,dec01_indices]=sort(dec01)

%% Task 3.3
linsapce(0,1,11)

%% Task  3.4
deciles=[dec01 dec02 dec03 dec04 dec05 dec06 dec07 dec08 dec09 dec10];
std(deciles)
std(deciles,[],1)

%% Task  3.4.1
std(deciles,[],2)

%% Task 3.5
% Contents of file normal_likelihood
%
%  function f=normal_likelihood(mu,sigma2,x)
%
%  f=1/(sqrt(2*pi*sigma2))*exp(-(x-mu)^2/(2*sigma2));

normal_likelihood(0,1,0)