eq1: (4 * x^(2)) + (4 * (x * y)) + y^(2) - (24 * x) + (38 * y) = 139
f = FirstAxis(eq1)
H = Intersect(eq1, f)
O = Intersect(xAxis, yAxis)
u = Vector(H, O)
eq1' = Translate(eq1, u)
α: atand(Slope[f])
β: 180° - α
eq1'_{1} = Rotate(eq1', β)