function [mean,stdev] = stat(x)
% STAT Interesting statistics. This documentaion is displayed by the command
% help stat
  n = length(x);
  mean = sum(x)/n;
  stdev = sqrt(sum((x-mean).^2)/(n-1));
