  function output=zop(Qy,dx,n)
% function output=zop(Qy,dx,n)
%   Qy=diag(y)-Q  (sparse)

dX=usMat(dx,n);
output=Qy*dX;
output=output(:);
