ODE
\[ y''''(x)+2 y''(x)+y(x)=\cos (x) \] ODE Classification
[[_high_order, _linear, _nonhomogeneous]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.18343 (sec), leaf count = 43
\[\left \{\left \{y(x)\to \left (-\frac {x^2}{8}+c_2 x+\frac {5}{16}+c_1\right ) \cos (x)+\frac {1}{4} (x+4 c_4 x+4 c_3) \sin (x)\right \}\right \}\]
Maple ✓
cpu = 1.143 (sec), leaf count = 38
\[\left [y \left (x \right ) = \left (-\frac {x^{2}}{8}+\frac {1}{4}\right ) \cos \left (x \right )+\frac {3 x \sin \left (x \right )}{8}+\textit {\_C1} \cos \left (x \right )+\textit {\_C2} \sin \left (x \right )+\textit {\_C3} x \cos \left (x \right )+\textit {\_C4} x \sin \left (x \right )\right ]\] Mathematica raw input
DSolve[y[x] + 2*y''[x] + y''''[x] == Cos[x],y[x],x]
Mathematica raw output
{{y[x] -> (5/16 - x^2/8 + C[1] + x*C[2])*Cos[x] + ((x + 4*C[3] + 4*x*C[4])*Sin[x
])/4}}
Maple raw input
dsolve(diff(diff(diff(diff(y(x),x),x),x),x)+2*diff(diff(y(x),x),x)+y(x) = cos(x), y(x))
Maple raw output
[y(x) = (-1/8*x^2+1/4)*cos(x)+3/8*x*sin(x)+_C1*cos(x)+_C2*sin(x)+_C3*x*cos(x)+_C
4*x*sin(x)]