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