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

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

[[_1st_order, _with_linear_symmetries], _dAlembert]

Book solution method
No Missing Variables ODE, Solve for x

Mathematica
cpu = 96.7453 (sec), leaf count = 407

{Solve[a2(x2+y(x)2)1+14ilog((a2x21)y(x)22xy(x)a2x2+a2y(x)21+a2y(x)4+x2)+14ilog((a2x21)y(x)2+2xy(x)a2x2+a2y(x)21+a2y(x)4+x2)+tan1(1a2(x2+y(x)2)1)12ilog(1a2y(x)2)12ilog(x2+y(x)2)+tan1(xy(x))=c1,y(x)],Solve[4a2x2+a2y(x)21+ilog(a2x2y(x)22xy(x)a2(x2+y(x)2)1+a2y(x)4+x2y(x)2)+ilog(a2x2y(x)2+2xy(x)a2(x2+y(x)2)1+a2y(x)4+x2y(x)2)+4tan1(1a2(x2+y(x)2)1)2ilog(a2y(x)21)+4c12ilog(x2+y(x)2)=4tan1(xy(x)),y(x)]}

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

{[x(_T)=_T(arctan(_T)a+1_Ta+_C1)1_T2+1,y(_T)=_C1a_T+arctan(_T)a1_T2+1]} Mathematica raw input

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

Mathematica raw output

{Solve[ArcTan[x/y[x]] + ArcTan[1/Sqrt[-1 + a^2*(x^2 + y[x]^2)]] - (I/2)*Log[x^2 
+ y[x]^2] - (I/2)*Log[1 - a^2*y[x]^2] + (I/4)*Log[x^2 + (-1 + a^2*x^2)*y[x]^2 + 
a^2*y[x]^4 - 2*x*y[x]*Sqrt[-1 + a^2*x^2 + a^2*y[x]^2]] + (I/4)*Log[x^2 + (-1 + a
^2*x^2)*y[x]^2 + a^2*y[x]^4 + 2*x*y[x]*Sqrt[-1 + a^2*x^2 + a^2*y[x]^2]] + Sqrt[-
1 + a^2*(x^2 + y[x]^2)] == C[1], y[x]], Solve[4*ArcTan[1/Sqrt[-1 + a^2*(x^2 + y[
x]^2)]] + 4*C[1] - (2*I)*Log[x^2 + y[x]^2] - (2*I)*Log[-1 + a^2*y[x]^2] + I*Log[
x^2 - y[x]^2 + a^2*x^2*y[x]^2 + a^2*y[x]^4 - 2*x*y[x]*Sqrt[-1 + a^2*(x^2 + y[x]^
2)]] + I*Log[x^2 - y[x]^2 + a^2*x^2*y[x]^2 + a^2*y[x]^4 + 2*x*y[x]*Sqrt[-1 + a^2
*(x^2 + y[x]^2)]] + 4*Sqrt[-1 + a^2*x^2 + a^2*y[x]^2] == 4*ArcTan[x/y[x]], y[x]]
}

Maple raw input

dsolve((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) = 1/(_T^2+1)^(1/2)*_T*(1/a*arctan(_T)+1/_T/a+_C1), y(_T) = -(_C1*a-_T+arc
tan(_T))/(_T^2+1)^(1/2)/a]