4.46.22 y′′′′′′(x)+ay(x)=0

ODE
y′′′′′′(x)+ay(x)=0 ODE Classification

[[_high_order, _missing_x]]

Book solution method
TO DO

Mathematica
cpu = 0.0083526 (sec), leaf count = 98

{{y(x)c2e(1)5/6a6x+c3e16a6x+c5e(1)5/6a6x+c6e16a6x+c4sin(a6x)+c1cos(a6x)}}

Maple
cpu = 0.013 (sec), leaf count = 102

{y(x)=_C1e(i31)x2a6+_C2e(i3+1)x2a6+_C3e(i31)x2a6+_C4e(i3+1)x2a6+_C5ea6x+_C6ea6x} Mathematica raw input

DSolve[a*y[x] + y''''''[x] == 0,y[x],x]

Mathematica raw output

{{y[x] -> E^((-1)^(5/6)*a^(1/6)*x)*C[2] + C[3]/E^((-1)^(1/6)*a^(1/6)*x) + C[5]/E
^((-1)^(5/6)*a^(1/6)*x) + E^((-1)^(1/6)*a^(1/6)*x)*C[6] + C[1]*Cos[a^(1/6)*x] + 
C[4]*Sin[a^(1/6)*x]}}

Maple raw input

dsolve(diff(diff(diff(diff(diff(diff(y(x),x),x),x),x),x),x)+a*y(x) = 0, y(x),'implicit')

Maple raw output

y(x) = _C1*exp(-1/2*(I*3^(1/2)-1)*(-a)^(1/6)*x)+_C2*exp(-1/2*(I*3^(1/2)+1)*(-a)^
(1/6)*x)+_C3*exp(1/2*(I*3^(1/2)-1)*(-a)^(1/6)*x)+_C4*exp(1/2*(I*3^(1/2)+1)*(-a)^
(1/6)*x)+_C5*exp(-(-a)^(1/6)*x)+_C6*exp((-a)^(1/6)*x)