prsq(x);

function val = prsq(x);
x = x(:);
val = [ ones(size(x)) zeros(size(x)) ones(size(x)) zeros(size(x))];