4.37.36 \(y''(x)=a y(x) \left (\left (b-y'(x)\right )^2+1\right )^{3/2}\)

ODE
\[ y''(x)=a y(x) \left (\left (b-y'(x)\right )^2+1\right )^{3/2} \] ODE Classification

[[_2nd_order, _missing_x], [_2nd_order, _reducible, _mu_x_y1]]

Book solution method
TO DO

Mathematica
cpu = 13.3632 (sec), leaf count = 1391

\[\left \{\left \{y(x)\to \text {InverseFunction}\left [-\frac {i \left (\sqrt {a^2 \left (b^2+1\right )}-a c_1\right ) \sqrt {\frac {-2 a^2 \text {$\#$1}^2-4 a c_1+4 \sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}} \sqrt {\frac {a^2 \text {$\#$1}^2+2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \sqrt {-\left (a \text {$\#$1}^2+2 c_1\right ){}^2 \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )} E\left (i \sinh ^{-1}\left (\sqrt {\frac {a^2}{2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}} \text {$\#$1}\right )|-\frac {a c_1+\sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}\right )-a b \sqrt {\frac {a^2}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \text {$\#$1} \left (a \text {$\#$1}^2+2 c_1\right ) \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )-i \sqrt {2} \sqrt {a^2 \left (b^2+1\right )} F\left (i \sinh ^{-1}\left (\sqrt {\frac {a^2}{2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}} \text {$\#$1}\right )|-\frac {a c_1+\sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}\right ) \sqrt {\frac {-a^2 \text {$\#$1}^2-2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}} \sqrt {\frac {a^2 \text {$\#$1}^2+2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \sqrt {-\left (a \text {$\#$1}^2+2 c_1\right ){}^2 \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )}}{a \left (b^2+1\right ) \sqrt {\frac {a^2}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \left (a \text {$\#$1}^2+2 c_1\right ) \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )}\& \right ][x+c_2]\right \},\left \{y(x)\to \text {InverseFunction}\left [\frac {i \left (\sqrt {a^2 \left (b^2+1\right )}-a c_1\right ) \sqrt {\frac {-2 a^2 \text {$\#$1}^2-4 a c_1+4 \sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}} \sqrt {\frac {a^2 \text {$\#$1}^2+2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \sqrt {-\left (a \text {$\#$1}^2+2 c_1\right ){}^2 \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )} E\left (i \sinh ^{-1}\left (\sqrt {\frac {a^2}{2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}} \text {$\#$1}\right )|-\frac {a c_1+\sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}\right )+a b \sqrt {\frac {a^2}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \text {$\#$1} \left (a \text {$\#$1}^2+2 c_1\right ) \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )-i \sqrt {2} \sqrt {a^2 \left (b^2+1\right )} F\left (i \sinh ^{-1}\left (\sqrt {\frac {a^2}{2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}} \text {$\#$1}\right )|-\frac {a c_1+\sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}\right ) \sqrt {\frac {-a^2 \text {$\#$1}^2-2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}{\sqrt {a^2 \left (b^2+1\right )}-a c_1}} \sqrt {\frac {a^2 \text {$\#$1}^2+2 a c_1+2 \sqrt {a^2 \left (b^2+1\right )}}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \sqrt {-\left (a \text {$\#$1}^2+2 c_1\right ){}^2 \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )}}{a \left (b^2+1\right ) \sqrt {\frac {a^2}{a c_1+\sqrt {a^2 \left (b^2+1\right )}}} \left (a \text {$\#$1}^2+2 c_1\right ) \left (a^2 \text {$\#$1}^4+4 a c_1 \text {$\#$1}^2-4 b^2+4 c_1{}^2-4\right )}\& \right ][x+c_2]\right \}\right \}\]

Maple
cpu = 1.067 (sec), leaf count = 431

\[\left [\int _{}^{y \left (x \right )}\frac {2 b}{\frac {2 \textit {\_a}^{2} \left (-a \,\textit {\_a}^{2}-2 \textit {\_C1} a +\sqrt {-\textit {\_a}^{4} a^{2} b^{2}-4 \textit {\_C1} \,\textit {\_a}^{2} a^{2} b^{2}-4 \textit {\_C1}^{2} a^{2} b^{2}+4 b^{4}+4 b^{2}}\right ) a}{\textit {\_a}^{4} a^{2}+4 \textit {\_C1} \,\textit {\_a}^{2} a^{2}+4 a^{2} \textit {\_C1}^{2}-4 b^{2}}+\frac {4 \textit {\_C1} \left (-a \,\textit {\_a}^{2}-2 \textit {\_C1} a +\sqrt {-\textit {\_a}^{4} a^{2} b^{2}-4 \textit {\_C1} \,\textit {\_a}^{2} a^{2} b^{2}-4 \textit {\_C1}^{2} a^{2} b^{2}+4 b^{4}+4 b^{2}}\right ) a}{\textit {\_a}^{4} a^{2}+4 \textit {\_C1} \,\textit {\_a}^{2} a^{2}+4 a^{2} \textit {\_C1}^{2}-4 b^{2}}+2 b^{2}+2}d \textit {\_a} -x -\textit {\_C2} = 0, \int _{}^{y \left (x \right )}\frac {2 b}{-\frac {2 \textit {\_a}^{2} \left (a \,\textit {\_a}^{2}+2 \textit {\_C1} a +\sqrt {-\textit {\_a}^{4} a^{2} b^{2}-4 \textit {\_C1} \,\textit {\_a}^{2} a^{2} b^{2}-4 \textit {\_C1}^{2} a^{2} b^{2}+4 b^{4}+4 b^{2}}\right ) a}{\textit {\_a}^{4} a^{2}+4 \textit {\_C1} \,\textit {\_a}^{2} a^{2}+4 a^{2} \textit {\_C1}^{2}-4 b^{2}}-\frac {4 \textit {\_C1} \left (a \,\textit {\_a}^{2}+2 \textit {\_C1} a +\sqrt {-\textit {\_a}^{4} a^{2} b^{2}-4 \textit {\_C1} \,\textit {\_a}^{2} a^{2} b^{2}-4 \textit {\_C1}^{2} a^{2} b^{2}+4 b^{4}+4 b^{2}}\right ) a}{\textit {\_a}^{4} a^{2}+4 \textit {\_C1} \,\textit {\_a}^{2} a^{2}+4 a^{2} \textit {\_C1}^{2}-4 b^{2}}+2 b^{2}+2}d \textit {\_a} -x -\textit {\_C2} = 0\right ]\] Mathematica raw input

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

Mathematica raw output

{{y[x] -> InverseFunction[-((-(a*b*Sqrt[a^2/(Sqrt[a^2*(1 + b^2)] + a*C[1])]*#1*(
2*C[1] + a*#1^2)*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]*#1^2 + a^2*#1^4)) + I*(Sqrt[a
^2*(1 + b^2)] - a*C[1])*EllipticE[I*ArcSinh[Sqrt[a^2/(2*Sqrt[a^2*(1 + b^2)] + 2*
a*C[1])]*#1], -((Sqrt[a^2*(1 + b^2)] + a*C[1])/(Sqrt[a^2*(1 + b^2)] - a*C[1]))]*
Sqrt[(4*Sqrt[a^2*(1 + b^2)] - 4*a*C[1] - 2*a^2*#1^2)/(Sqrt[a^2*(1 + b^2)] - a*C[
1])]*Sqrt[(2*Sqrt[a^2*(1 + b^2)] + 2*a*C[1] + a^2*#1^2)/(Sqrt[a^2*(1 + b^2)] + a
*C[1])]*Sqrt[-((2*C[1] + a*#1^2)^2*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]*#1^2 + a^2*
#1^4))] - I*Sqrt[2]*Sqrt[a^2*(1 + b^2)]*EllipticF[I*ArcSinh[Sqrt[a^2/(2*Sqrt[a^2
*(1 + b^2)] + 2*a*C[1])]*#1], -((Sqrt[a^2*(1 + b^2)] + a*C[1])/(Sqrt[a^2*(1 + b^
2)] - a*C[1]))]*Sqrt[(2*Sqrt[a^2*(1 + b^2)] - 2*a*C[1] - a^2*#1^2)/(Sqrt[a^2*(1 
+ b^2)] - a*C[1])]*Sqrt[(2*Sqrt[a^2*(1 + b^2)] + 2*a*C[1] + a^2*#1^2)/(Sqrt[a^2*
(1 + b^2)] + a*C[1])]*Sqrt[-((2*C[1] + a*#1^2)^2*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[
1]*#1^2 + a^2*#1^4))])/(a*(1 + b^2)*Sqrt[a^2/(Sqrt[a^2*(1 + b^2)] + a*C[1])]*(2*
C[1] + a*#1^2)*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]*#1^2 + a^2*#1^4))) & ][x + C[2]
]}, {y[x] -> InverseFunction[(a*b*Sqrt[a^2/(Sqrt[a^2*(1 + b^2)] + a*C[1])]*#1*(2
*C[1] + a*#1^2)*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]*#1^2 + a^2*#1^4) + I*(Sqrt[a^2
*(1 + b^2)] - a*C[1])*EllipticE[I*ArcSinh[Sqrt[a^2/(2*Sqrt[a^2*(1 + b^2)] + 2*a*
C[1])]*#1], -((Sqrt[a^2*(1 + b^2)] + a*C[1])/(Sqrt[a^2*(1 + b^2)] - a*C[1]))]*Sq
rt[(4*Sqrt[a^2*(1 + b^2)] - 4*a*C[1] - 2*a^2*#1^2)/(Sqrt[a^2*(1 + b^2)] - a*C[1]
)]*Sqrt[(2*Sqrt[a^2*(1 + b^2)] + 2*a*C[1] + a^2*#1^2)/(Sqrt[a^2*(1 + b^2)] + a*C
[1])]*Sqrt[-((2*C[1] + a*#1^2)^2*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]*#1^2 + a^2*#1
^4))] - I*Sqrt[2]*Sqrt[a^2*(1 + b^2)]*EllipticF[I*ArcSinh[Sqrt[a^2/(2*Sqrt[a^2*(
1 + b^2)] + 2*a*C[1])]*#1], -((Sqrt[a^2*(1 + b^2)] + a*C[1])/(Sqrt[a^2*(1 + b^2)
] - a*C[1]))]*Sqrt[(2*Sqrt[a^2*(1 + b^2)] - 2*a*C[1] - a^2*#1^2)/(Sqrt[a^2*(1 + 
b^2)] - a*C[1])]*Sqrt[(2*Sqrt[a^2*(1 + b^2)] + 2*a*C[1] + a^2*#1^2)/(Sqrt[a^2*(1
 + b^2)] + a*C[1])]*Sqrt[-((2*C[1] + a*#1^2)^2*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]
*#1^2 + a^2*#1^4))])/(a*(1 + b^2)*Sqrt[a^2/(Sqrt[a^2*(1 + b^2)] + a*C[1])]*(2*C[
1] + a*#1^2)*(-4 - 4*b^2 + 4*C[1]^2 + 4*a*C[1]*#1^2 + a^2*#1^4)) & ][x + C[2]]}}

Maple raw input

dsolve(diff(diff(y(x),x),x) = a*y(x)*(1+(b-diff(y(x),x))^2)^(3/2), y(x))

Maple raw output

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