function m = MCanl(data)
%MCanl(data) returns the sample mean and prints a simple
%summary of the data and 95% confidence error.
N=length(data);
m=mean(data);
v=var(data);
fprintf('\nSample variance=%g\n',v)
fprintf('Probable error=%g\n',.6745*sqrt(v/N))
fprintf('95%% interval=[%g, %g]\n',1.96*sqrt(v/N)*[-1,1]+m)
end

