[[_high_order, _quadrature]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.0351655 (sec), leaf count = 32
Maple ✓
cpu = 0.023 (sec), leaf count = 29
DSolve[y''''[x] == x*Cos[x],y[x],x]
Mathematica raw output
{{y[x] -> C[1] + x*C[2] + x^2*C[3] + x^3*C[4] + x*Cos[x] - 4*Sin[x]}}
Maple raw input
dsolve(diff(diff(diff(diff(y(x),x),x),x),x) = x*cos(x), y(x),'implicit')
Maple raw output
y(x) = -4*sin(x)+x*cos(x)+1/6*x^3*_C1+1/2*x^2*_C2+_C3*x+_C4