The question is to solve $x\frac{\mathrm{d} y}{\mathrm{d} x}+3y=x^3y^2$.
I have found the homogenous solution $y_h = c_1 x^{-3}$
I am stuck at finding the particular solution. I am familiar with variation of parameters(which involve Wronskian and just $r(x)$ in RHS ) and solution by undetermined coefficients
Please do as well suggest any necessary reading required for the same
Soham