function z = proj(x,y)
    z = x'*y / (y' * y) * y;
end