Answer on Question #69645 – Math – Differential Equations
Question
Solve the Riccati's equation
d y d x = − 1 − x 2 + y 2 . \frac{dy}{dx} = -1 - x^2 + y^2. d x d y = − 1 − x 2 + y 2 . Solution
One searches for the general solution in the following form
y ( x ) = z ( x ) + y p ( x ) , y(x) = z(x) + y_p(x), y ( x ) = z ( x ) + y p ( x ) ,
where y p ( x ) y_p(x) y p ( x ) is a particular solution of the given equation.
This substitution transforms the equation into Bernoulli's one of the unknown function z ( x ) z(x) z ( x ) .
Let the particular solution be
y p ( x ) = a x . y_p(x) = ax. y p ( x ) = a x .
Substituting it into the differential equation we get
d d x ( a x ) = − 1 − x 2 + ( a x ) 2 \frac{d}{dx}(ax) = -1 - x^2 + (ax)^2 d x d ( a x ) = − 1 − x 2 + ( a x ) 2
or
a = − 1 + ( a 2 − 1 ) x 2 . a = -1 + (a^2 - 1)x^2. a = − 1 + ( a 2 − 1 ) x 2 .
Value
a = − 1 a = -1 a = − 1
satisfies the previous equation.
The particular solution is
y p ( x ) = − x , y_p(x) = -x, y p ( x ) = − x ,
and the general solution may be found as
y ( x ) = z ( x ) − x . y(x) = z(x) - x. y ( x ) = z ( x ) − x .
Substituting it into the Riccati's equation we get
d d x ( z ( x ) − x ) = − 1 − x 2 + ( z ( x ) − x ) 2 , \frac{d}{dx}(z(x) - x) = -1 - x^2 + (z(x) - x)^2, d x d ( z ( x ) − x ) = − 1 − x 2 + ( z ( x ) − x ) 2 , d z d x − 1 = − 1 − x 2 + z 2 − 2 x z + x 2 , \frac {d z}{d x} - 1 = - 1 - x ^ {2} + z ^ {2} - 2 x z + x ^ {2}, d x d z − 1 = − 1 − x 2 + z 2 − 2 x z + x 2 , d z d x = z 2 − 2 x z . \frac {d z}{d x} = z ^ {2} - 2 x z. d x d z = z 2 − 2 x z .
This is Bernoulli's equation. It can be solved using the substitution z ( x ) = u ( x ) v ( x ) z(x) = u(x)v(x) z ( x ) = u ( x ) v ( x ) :
d ( u v ) d x = ( u v ) 2 − 2 x u v , \frac {d (u v)}{d x} = (u v) ^ {2} - 2 x u v, d x d ( uv ) = ( uv ) 2 − 2 xuv , v d u d x + u d v d x + 2 x u v − u 2 v 2 = 0 , v \frac {d u}{d x} + u \frac {d v}{d x} + 2 x u v - u ^ {2} v ^ {2} = 0, v d x d u + u d x d v + 2 xuv − u 2 v 2 = 0 , ( d u d x + 2 x u ) v + ( u d v d x − u 2 v 2 ) = 0 , \left(\frac {d u}{d x} + 2 x u\right) v + \left(u \frac {d v}{d x} - u ^ {2} v ^ {2}\right) = 0, ( d x d u + 2 xu ) v + ( u d x d v − u 2 v 2 ) = 0 ,
We get the system
{ d u d x + 2 x u = 0 u d v d x − u 2 v 2 = 0 . \left\{ \begin{array}{l} \frac {d u}{d x} + 2 x u = 0 \\ u \frac {d v}{d x} - u ^ {2} v ^ {2} = 0 \end{array} \right.. { d x d u + 2 xu = 0 u d x d v − u 2 v 2 = 0 .
We need a particular solution of the first equation of the system. Check that it is u ( x ) = b e − x 2 u(x) = be^{-x^2} u ( x ) = b e − x 2 :
d ( b e − x 2 ) d x + 2 b x e − x 2 = 0 , \frac {d \left(b e ^ {- x ^ {2}}\right)}{d x} + 2 b x e ^ {- x ^ {2}} = 0, d x d ( b e − x 2 ) + 2 b x e − x 2 = 0 , − 2 x b e − x 2 + 2 b x e − x 2 = 0 , - 2 x b e ^ {- x ^ {2}} + 2 b x e ^ {- x ^ {2}} = 0, − 2 x b e − x 2 + 2 b x e − x 2 = 0 , 0 = 0. 0 = 0. 0 = 0.
Let b = 1 b = 1 b = 1 , so a particular solution is u ( x ) = e − x 2 u(x) = e^{-x^2} u ( x ) = e − x 2 .
As u ( x ) = e − x 2 ≠ 0 u(x) = e^{-x^2} \neq 0 u ( x ) = e − x 2 = 0 , the second equation may be reduced to
d v d x − u v 2 = 0 \frac {d v}{d x} - u v ^ {2} = 0 d x d v − u v 2 = 0
or
d v d x − e − x 2 v 2 = 0. \frac {d v}{d x} - e ^ {- x ^ {2}} v ^ {2} = 0. d x d v − e − x 2 v 2 = 0.
This is the separable DE.
∫ d v v 2 = ∫ e − x 2 d x , \int \frac {d v}{v ^ {2}} = \int e ^ {- x ^ {2}} d x, ∫ v 2 d v = ∫ e − x 2 d x , − 1 v ( x ) = π 2 erf x − C 2 , - \frac {1}{v (x)} = \frac {\sqrt {\pi}}{2} \operatorname {erf} x - \frac {C}{2}, − v ( x ) 1 = 2 π erf x − 2 C ,
where C C C is an integration constant.
Then
v ( x ) = − 1 π 2 erf x − C 2 = 2 C − π erf x v (x) = - \frac {1}{\frac {\sqrt {\pi}}{2} \operatorname {erf} x - \frac {C}{2}} = \frac {2}{C - \sqrt {\pi} \operatorname {erf} x} v ( x ) = − 2 π erf x − 2 C 1 = C − π erf x 2
and
z ( x ) = u ( x ) v ( x ) = 2 e − x 2 C − π erf x . z (x) = u (x) v (x) = \frac {2 e ^ {- x ^ {2}}}{C - \sqrt {\pi} \operatorname {erf} x}. z ( x ) = u ( x ) v ( x ) = C − π erf x 2 e − x 2 .
Finally, the general solution of the Riccati's equation is
y ( x ) = z ( x ) − x = 2 e − x 2 C − π erf x − x . y (x) = z (x) - x = \frac {2 e ^ {- x ^ {2}}}{C - \sqrt {\pi} \operatorname {erf} x} - x. y ( x ) = z ( x ) − x = C − π erf x 2 e − x 2 − x .
Answer: y ( x ) = 2 e − x 2 C − π erf x − x . y(x) = \frac{2e^{-x^2}}{C - \sqrt{\pi}\operatorname{erf}x} - x. y ( x ) = C − π erf x 2 e − x 2 − x .
Answer provided by https://www.AssignmentExpert.com