samplevar.m

function s2 = samplevar(x)
%Input
% x data set (vector)
%
%Output
% s2 variance (scaler)
% s2 = 1/(n-1) * sum((x_i - x_bar)^2)

n = length(x);

% Calculates mean
x_bar = (1/n)*sum(x);

% Creates an empty vector to later be summed together
s2_vec = zeros(1,n);

for i=1:n
    
    %Calculates each component that needs to be summed
    s2_vec(i) = (x(i)-x_bar)^2;
end

%Calculates variance
s2 = (1/(n-1))*sum(s2_vec);
end