[[_2nd_order, _missing_x]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.155287 (sec), leaf count = 101
Maple ✓
cpu = 0.293 (sec), leaf count = 92
DSolve[y''[x]^2 == a + b*y'[x]^2,y[x],x]
Mathematica raw output
{{y[x] -> -(a*b*E^(2*Sqrt[b]*x) + E^(2*Sqrt[b]*C[1]))/(2*b^(3/2)*E^(Sqrt[b]*(x +
C[1]))) + C[2]}, {y[x] -> (a*b + E^(2*Sqrt[b]*(x + C[1])))/(2*b^(3/2)*E^(Sqrt[b
]*(x + C[1]))) + C[2]}}
Maple raw input
dsolve(diff(diff(y(x),x),x)^2 = a+b*diff(y(x),x)^2, y(x),'implicit')
Maple raw output
y(x) = -1/b*(-a*b)^(1/2)*x+_C1, y(x) = 1/b*(-a*b)^(1/2)*x+_C1, y(x) = _C1+_C2*ex
p(b^(1/2)*x)+1/4*a/b^2/_C2*exp(-b^(1/2)*x), y(x) = _C1+1/4*a/b^2/_C2*exp(b^(1/2)
*x)+_C2*exp(-b^(1/2)*x)