Question

I am trying to solve the system of differential equations dx/dt = (y-x)/(x+y), dy/dt = -y/(x+y), where x and y are functions of t, and x(0)=0, y(0)=c (positive constant). I would like to find x and y explicitly in terms of t.

This problem is related to my previous question "Probability problem with solution involving e."

Any help would be appreciated. Thanks.

Answer 1

Continuing Michael's answer: Solve $dy/dx = y/(x-y), y(0)=c$ to get (in terms of the Lambert W function) $y(x) = \operatorname{e} ^{W(-x/c) + \mathrm{log} (c)}$. Substitute this into the original equations to get a single equation: $$ \frac{d x (t)}{d t} = \frac{\operatorname{e} ^{\Bigl(W \Bigl(-\frac{x (t)}{c}\Bigr) + \operatorname{log} (c)\Bigr)} - x (t)}{\operatorname{e} ^{\Bigl(W \Bigl(-\frac{x (t)}{c}\Bigr) + \operatorname{ln} (c)\Bigr)} + x (t)} $$ $x(0)=0$. Maple finds two solutions for this

added

The one that seems to work for real $c > 0$ is: $$x(t) = -\left(W_{-1} \left(\frac{-(-t + 2 c) \operatorname{e} ^{-2}}{c}\right) + 2\right) \operatorname{e} ^{\Biggl(tW_{-1} \Biggl(\frac{-(-t + 2 c) \operatorname{e} ^{-2}}{c}\Biggr) + 2\Biggr)} c $$

Answer 2

I think it might be easier the other way round: solve $dx(y)/dy=x(y)/y-1$, $x(c)=0$ to get $x(y)=(\log c-\log y)y$. Then we are left with $$\frac{dy(t)}{dt}=-\frac1{\log c+1-\log y(t)}=\frac1{\log(y(t)/ec)}.$$ This is solved by $y(t)=\frac{t+d}{W((t+d)/ce^2)}$. Plugging in $y(0)=c$ gives $d=-2c$ for the lower branch of $W$, hence, barring mistakes, \begin{align} y(t)&=\frac{t-2c}{W_{-1}\left(\frac{t-2c}{ce^2}\right)},\\\\ x(t)&=2(c-y(t))-t, \end{align} for $-c(e-2)< t< 2c$. Note that $\lim_{t\to2c-}(x(t),y(t))=(0,0)$ and $\lim_{t\to(2-e)c+}(x(t),y(t))=(-ce,ce)$.

Answer 3

$$dx/dy=x/y-1.$$ From this, you find $$\frac{d}{dy}(\frac{x}{y})=-\frac{1}{y}.$$