function y = sem(x,varargin)
if isempty(varargin);
  dim=1;
elseif size(nargin,2)==1;
  dim=varargin{1};
end
n = size(x);
n = n(dim);
y = nanstd(x,[],dim) / sqrt(n);