function uzero=zerolev(u,nlevel)
%function uzero=zerolev(u,nlevel)
%   subtracts u(:,nlevel) from u
[m,n]=size(u);
uzero=u-u(:,nlevel)*ones(1,n);
