(xdy)/(dx)=y+sqrt(x^2-y^2)
⚠ (xdy)/(dx)=y+sqrt(x^2-y^2) is read as (x)*y'=y+sqrt(x^2-y^2): use dif(expression, variable) for a derivative with respect to that variable, use (numerator)/(denominator) for a quotient
1) Let's solve the differential equation x·y' = y + x2 -y2, where y is a function of x; on an interval where x > 0, isolate y': y' = y + x2 -y2x
On this interval, the right side only depends on yx: the equation is homogeneous. With v = yx, it reads y' = v + 1 -v2
Let y = vx, where v is a function of x: then y' = v + x·v', and the equation becomes v + x·v' = v + 1 -v2, that is x·v' = 1 -v2
The right side is defined for v ∈ [-1, 1]; this condition remains necessary for each branch
Before dividing, 1 -v2 = 0 gives the equilibrium solutions v = { -1, 1 }, that is (y = -x ∨ y = x)
On intervals where 1 -v2 ≠ 0, separate the variables: dv1 -v2 = dxx
⇥1.1) Let's integrate: ∫11 -v2dv
⇥ⓘ The derivative of arcsin(v) is 11 -v2, so arcsin(v) is an antiderivative of 11 -v2: ∫11 -v2dv = arcsin(v)
⇥ arcsin(v)
Integrating both sides, with ∫1xdx = ln|x| since x ≠ 0: arcsin(v) = ln|x| + C, where C is an arbitrary real constant, chosen independently on each interval
Isolate v: v = sin(ln|x| + C). Since y = vx, the solution branches are:
These branches hold on intervals where x > 0. At points where 1 -v2 = 0, a join with an equilibrium must preserve the value and derivative and satisfy the original equation
On an interval where x < 0, isolate y': y' = y + x2 -y2x
On this interval, the right side only depends on yx: the equation is homogeneous. With v = yx, it reads y' = v -1 -v2
Let y = vx, where v is a function of x: then y' = v + x·v', and the equation becomes v + x·v' = v -1 -v2, that is x·v' = -1 -v2
The right side is defined for v ∈ [-1, 1]; this condition remains necessary for each branch
Before dividing, -1 -v2 = 0 gives the equilibrium solutions v = { -1, 1 }, that is (y = -x ∨ y = x)
On intervals where -1 -v2 ≠ 0, separate the variables: -dv1 -v2 = dxx
⇥1.2) Let's compute the integral: ∫(-11 -v2)dv
⇥ⓘ The derivative of arccos(v) is -11 -v2, so arccos(v) is an antiderivative of -11 -v2: ∫(-11 -v2)dv = arccos(v)
⇥ arccos(v)
Integrating both sides, with ∫1xdx = ln|x| since x ≠ 0: arccos(v) = ln|x| + C, where C is an arbitrary real constant, chosen independently on each interval
Isolate v: v = cos(ln|x| + C). Since y = vx, the solution branches are:
These branches hold on intervals where x < 0. At points where -1 -v2 = 0, a join with an equilibrium must preserve the value and derivative and satisfy the original equation
▶Solution:   ((x > 0 ∧ (y = -x ∨ y = x ∨ y = x·sin(ln|x| + C))) ∨ (x < 0 ∧ (y = -x ∨ y = x ∨ y = x·cos(ln|x| + C))))