function [Px P]=ZeroOut(x,index)

y=zeros(size(x));
y(index)=norm(x);

if nargout == 1   % Don't waste time computing P uless the user asks for it

  Px=reflect(x,(x-y));
    
else 
    
  [Px P]=reflect(x,(x-y));

end