Since the equation is not exact, we need to find a function µ(x, y), such that
µ(3x2y + 2xy + y3) + µ(x2 + y2)y’ = 0 is exact.
M = µ(3x2y + 2xy + y3) N = µ(x2 + y2)
Then My = µy’(3x2y + 2xy + y3) + µ(3x2 + 2x + 3y2)
Nx = µx’(x2 + y2) + µ(2x)
Let My = Nx , we get µy’(3x2y + 2xy + y3) + µ(3x2 + 3y2) = µx’(x2 + y2)
First suppose µ is a function of x only.
Then µy’ = 0, we get
µ(3x2 + 3y2) = µx’(x2 + y2)
3µ = ∂µ/∂x
By separable equation,
∫(1/µ) ∂µ = ∫3 ∂x
ln(µ) = 3x
µ = e3x (the integration factor)
By theorem, there exist ϕ(x,y) , such that , ϕx = M, ϕy = N
ϕ = ∫e3x(3x2y + 2xy + y3) = e3xx2y + (e3xy3)/3 + h(y)
ϕy = e3xx2 + h'(y) = e3xx2 + e3xy2
Thus, h'(y) = e3xy2, h(y) = (e3xy3)/3 + C
Therefore, the general solution is ϕ = e3xx2y + (e3xy3)/3 = C