added params arg, replaced slow diag mults by bsxfun
only returns full windows; uses bsxfun not repmat
accelerated, including bsxfun
paperfig1 plots a series of pics to illustrate process