DSolve[-y[x] + x*Sin[x*y[x]]*y[x]^2 + (-4*x + x^2*Sin[x*y[x]]*y[x])*Derivative[1][y][x] == 0,y[x],x]