Интегральный множитель получается равен m= exp ^ [int [[(dP/dy)-(dQ/dx)] / Q]dx] =
= x^[(1-4*y)/(2*y-1)]



А возможно его вычислить другим способом, чтобы он не был таким громоздким?