DSolve[a*y[x] + (b + 2*a*x)*Derivative[1][y][x] + 3*(-1 + 2*x)*Derivative[2][y][x] + 2*(-1 + x)*x*Derivative[3][y][x] == 0,y[x],x]