%std dev
DATA2 = importdata('1_2Reference_DOT_DECIMAL.txt');
i = DATA2(:,1);
T_ref = DATA2(:,2);

Tmean = mean(T_ref)

tsum = 0;
for elm = i' + 1
    tsum = tsum + ((T_ref(elm, 1) - Tmean).^2);
end
tsum = sum(tsum)
N = size(T_ref);
N = N(:,1);
stdavv = tsum / (N - 1)
stdavv = sqrt(stdavv)



std(T_ref)

%pd = fitdist(T_ref, 'Normal')
%paramci(pd)