3.5.3.21 Example 21
yy1+121+(y)2=x

Let y=p and rearranging gives

py=x(1+121+p2)y=x(1p+12p1+p2)=x(22p+12p1+p2)=x(2+1+p22p)(1)=xf+g

Hence

f=2+1+p22pg=0

Since f(p)p then this is d’Almbert ode. Taking derivative of (1) w.r.t. x gives

p=ddx(xf(p)+g(p))=f(p)+xf(p)dpdx+g(p)dpdx

But f(p)=2+1+p22p,f(p)=1p2,g=0,g=0 and the above becomes

p=2+1+p22p+x(121+p221+p22p2)dpdx(2)p+2+1+p22p=x(121+p221+p22p2)dpdx

The singular solution is found by setting dpdx=0 which results in p+2+1+p22p=0. Hence p=±i or y=±i   or y=±ix. But these do not satisfy the ode, hence no singular solutions exist.

The general solution is when dpdx0 in (2). This gives the ode

dpdx=1x(p+2+1+p22p)(121+p221+p22p2)=1x(p3+p)

But this is non-linear in p. Hence inversion is needed. This becomes

dxdp=x(121+p221+p22p2)(p+2+1+p22p)dxdp=xp3+pdxdp1p+p3x=0

Which is now linear in x(p).  The solution is

(3)x=p1+p2c1

We now need to eliminate p. We have two equations to do that, (1) and (3). Here they are side by side

(1)y=x(2+1+p22p)(3)x=p1+p2c1

We can either solve for p from (1) and plugin in the value found into (3). Or we can solve for p from (3) and plugin the value found in (1). Using CAS we can just use the solve command. For an example, using Maple it gives

eq1:=y=-x*(  (2+sqrt(1+p^2))/(2*p)); 
eq2:=x=p/sqrt(1+p^2)*_C1 
sol:=solve([eq1,eq2],[p,y],'allsolutions'); 
[[p = x*RootOf((c__1^2 - x^2)*_Z^2 - 1), y = -(RootOf((c__1^2 - x^2)*_Z^2 - 1)*c__1 + 2)/(2*RootOf((c__1^2 - x^2)*_Z^2 - 1))]]
 

Now we can use allvalues

map(X->allvalues(X),sol) 
[[p = x*sqrt(1/(c__1^2 - x^2)), y = -(sqrt(1/(c__1^2 - x^2))*c__1 + 2)/(2*sqrt(1/(c__1^2 - x^2)))], 
[p = -x*sqrt(1/(c__1^2 - x^2)), y = (-sqrt(1/(c__1^2 - x^2))*c__1 + 2)/(2*sqrt(1/(c__1^2 - x^2)))]]
 

Hence the solutions are

y1=1c12x2c1+221c12x2y2=1c12x2c1+221c12x2

These are verified valid solutions to the ode (had to use assuming positive)