DSolve[y''[x] - 4 y'[x] + 4 y[x] ==Exp[x] x^(-4), x, y[x]]