function y=omatrix(x,i,j)
%makes a matrix with zeros everywhere except for a one at the i,j entry
y=zeros(size(x));
y(i,j)=1;
