[[_2nd_order, _missing_x], [_2nd_order, _reducible, _mu_x_y1]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.209895 (sec), leaf count = 111
Maple ✓
cpu = 0.389 (sec), leaf count = 54
DSolve[y[x]*y''[x] == a,y[x],x]
Mathematica raw output
{{y[x] -> E^(-(C[1] + 2*a*InverseErf[(-I)*Sqrt[2/Pi]*Sqrt[a*E^(C[1]/a)*(x + C[2]
)^2]]^2)/(2*a))}, {y[x] -> E^(-(C[1] + 2*a*InverseErf[I*Sqrt[2/Pi]*Sqrt[a*E^(C[1
]/a)*(x + C[2])^2]]^2)/(2*a))}}
Maple raw input
dsolve(y(x)*diff(diff(y(x),x),x) = a, y(x),'implicit')
Maple raw output
Intat(1/(2*a*ln(_a)-2*_C1*a)^(1/2),_a = y(x))-x-_C2 = 0, Intat(-1/(-2*a*(_C1-ln(
_a)))^(1/2),_a = y(x))-x-_C2 = 0