4.39.35 \((a+y(x)) y''(x)+y'(x)^2=b\)

ODE
\[ (a+y(x)) y''(x)+y'(x)^2=b \] ODE Classification

[[_2nd_order, _missing_x], [_2nd_order, _exact, _nonlinear], [_2nd_order, _reducible, _mu_x_y1], [_2nd_order, _reducible, _mu_xy]]

Book solution method
TO DO

Mathematica
cpu = 0.484945 (sec), leaf count = 75

\[\left \{\left \{y(x)\to -\frac {a b+\sqrt {b^3 (x+c_2){}^2-b e^{2 c_1}}}{b}\right \},\left \{y(x)\to \frac {-a b+\sqrt {b^3 (x+c_2){}^2-b e^{2 c_1}}}{b}\right \}\right \}\]

Maple
cpu = 0.413 (sec), leaf count = 53

\[\left [y \left (x \right ) = -a -\sqrt {b \,x^{2}-2 \textit {\_C1} x +a^{2}+2 \textit {\_C2}}, y \left (x \right ) = -a +\sqrt {b \,x^{2}-2 \textit {\_C1} x +a^{2}+2 \textit {\_C2}}\right ]\] Mathematica raw input

DSolve[y'[x]^2 + (a + y[x])*y''[x] == b,y[x],x]

Mathematica raw output

{{y[x] -> -((a*b + Sqrt[-(b*E^(2*C[1])) + b^3*(x + C[2])^2])/b)}, {y[x] -> (-(a*
b) + Sqrt[-(b*E^(2*C[1])) + b^3*(x + C[2])^2])/b}}

Maple raw input

dsolve((a+y(x))*diff(diff(y(x),x),x)+diff(y(x),x)^2 = b, y(x))

Maple raw output

[y(x) = -a-(b*x^2-2*_C1*x+a^2+2*_C2)^(1/2), y(x) = -a+(b*x^2-2*_C1*x+a^2+2*_C2)^
(1/2)]