function ma=moving(y,wlen)
% ma=MOVING(y,wlen)
%
% Moving average routine of signal 'y' with window length 'wlen'
%
% Written by fjsimons-at-alum.mit.edu, 05/18/2010

y = y(:);

[m n] = size(y);
if m < wlen
  disp('Window wider than sample');
  ma=0;
else
  % Don't attempt to fix for NaN's
  ma=cumsum([sum(y(1:wlen)); y(wlen+1:m)-y(1:m-wlen)])./wlen;
end