4.23.36 ay(x)y(x)ax+y(x)y(x)2+1=0

ODE
ay(x)y(x)ax+y(x)y(x)2+1=0 ODE Classification

[[_homogeneous, `class A`], _rational, _dAlembert]

Book solution method
No Missing Variables ODE, Solve for y

Mathematica
cpu = 0.306254 (sec), leaf count = 148

{{y(x)e2c1((a21)3(e2c1)x2+2(a21)xe(a2+1)c1+e2a2c1)(a21)3},{y(x)e2c1((a21)3(e2c1)x2+2(a21)xe(a2+1)c1+e2a2c1)(a21)3}}

Maple
cpu = 0.262 (sec), leaf count = 129

{[x(_T)=ea(_T1_T2+1+a)(a_T_T2+1)1(a_T2_T_T2+1+a)1d_T_C1,y(_T)=_C1aea(_T1_T2+1+a)(a_T_T2+1)1(a_T2_T_T2+1+a)1d_T(a_T+_T2+1)1]} Mathematica raw input

DSolve[-(a*x) - a*y[x]*y'[x] + y[x]*Sqrt[1 + y'[x]^2] == 0,y[x],x]

Mathematica raw output

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

Maple raw input

dsolve(y(x)*(1+diff(y(x),x)^2)^(1/2)-a*y(x)*diff(y(x),x)-a*x = 0, y(x),'implicit')

Maple raw output

[x(_T) = exp(Int(a*(-1/(_T^2+1)^(1/2)*_T+a)/(a*_T-(_T^2+1)^(1/2))/(a*_T^2-_T*(_T
^2+1)^(1/2)+a),_T))*_C1, y(_T) = a*exp(Int(a*(-1/(_T^2+1)^(1/2)*_T+a)/(a*_T-(_T^
2+1)^(1/2))/(a*_T^2-_T*(_T^2+1)^(1/2)+a),_T))*_C1/(-a*_T+(_T^2+1)^(1/2))]