✓ Mathematica : cpu = 0.0237825 (sec), leaf count = 36
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]
✓ Maple : cpu = 0.04 (sec), leaf count = 21
dsolve(x*diff(diff(y(x),x),x)+(2*a*x*ln(x)+1)*diff(y(x),x)+(a^2*x*ln(x)^2+a*ln(x)+a)*y(x)=0,y(x))