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