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