4.45.24 $y^{⁗}(x)-2y^{″}(x)+y(x)=e^x+4$

ODE
$$y^{⁗}(x)-2y^{″}(x)+y(x)=e^x+4$$ ODE Classification

[[_high_order, _with_linear_symmetries]]

Book solution method
TO DO

Mathematica ✓
cpu = 0.0528111 (sec), leaf count = 55

$$\left\{\{y(x)\to \frac{1}{16}{e}^{-x}(e^{2x}(4(4c_4-1)x+16c_3+2x^2+3)+16(c_{2}x+c_1)+64e^x)\}\right\}$$

Maple ✓
cpu = 0.024 (sec), leaf count = 41

$$\{y\left(x\right)=\frac{(16\mathit{\_}\mathit{C}\mathit{4}x+16\mathit{\_}\mathit{C}\mathit{2}){\mathrm{e}}^{-x}}{16}+4+\frac{(2x^2+(16\mathit{\_}\mathit{C}\mathit{3}-4)x+16\mathit{\_}\mathit{C}\mathit{1}+3){\mathrm{e}}^x}{16}\}$$ Mathematica raw input

DSolve[y[x] - 2*y''[x] + y''''[x] == 4 + E^x,y[x],x]

Mathematica raw output

{{y[x] -> (64*E^x + 16*(C[1] + x*C[2]) + E^(2*x)*(3 + 2*x^2 + 16*C[3] + 4*x*(-1 
+ 4*C[4])))/(16*E^x)}}

Maple raw input

dsolve(diff(diff(diff(diff(y(x),x),x),x),x)-2*diff(diff(y(x),x),x)+y(x) = 4+exp(x), y(x),'implicit')

Maple raw output

y(x) = 1/16*(16*_C4*x+16*_C2)*exp(-x)+4+1/16*(2*x^2+(16*_C3-4)*x+16*_C1+3)*exp(x
)