Here,
dy/dx + xy = x3y3
Therefore
1/y3(dy/dx) + x.1/y2 = x3 (1)
Put 1/y2 = z. Then
-2/y3(dy/dx) = dz/dx
Therefore, Eq. (1) becomes
- 1/2(dz/dx) + xz = x3
⇒ dz/dx - 2xz = -2x3
It is in the linear form. So, the integrating factor e∫-2xdx = e-x^2.
Multiplying by it,

