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