DSolve[x y'[x] == y[x] + 2 x^3*Sin[y[x]/x]^2, y[x], x]