✓ Mathematica : cpu = 0.0187738 (sec), leaf count = 18
DSolve[(1 + 2*x)*y[x] - x*(1 + 2*x)*Derivative[1][y][x] + x^2*(1 + x)*Derivative[2][y][x] == 0,y[x],x]
✓ Maple : cpu = 0.03 (sec), leaf count = 15
dsolve(x^2*(1+x)*diff(diff(y(x),x),x)-x*(2*x+1)*diff(y(x),x)+(2*x+1)*y(x)=0,y(x))