y' = (x^2 + xy + y^2)/x