Trying to get exact solution for a nonlinear DE

  • Thread starter Thread starter Paalfaal
  • Start date Start date
  • Tags Tags
    Nonlinear
Paalfaal
Messages
13
Reaction score
0
Ok, this should be quite simple. I've been looking at this problem for quite some time now, and I'm tired.. Please help me!

The equation to solve is r' = r(1-r^2)

The obvious thing to do is to do partial fraction expansion and integrate(from r_{0} to r):

∫ (\frac{1}{r} + \frac{1}{1-r^2})dr = t

After some trig substitution: ln[r\sqrt{\frac{r+1}{1-r}}] = t (evaluated from r_{0} to r). Then take the exponential on both sides.

r\sqrt{\frac{r+1}{1-r}} = r_{0}\sqrt{\frac{r_{0}+1}{1-r_{0}}}e^{t}

This leads to kind of a nasty expression which can't be solved explicitly for r(t) (at least it seems that it can't be solved explicitly).. So this is where I'm stuck.

The solution to this ploblem is: r(t) = \frac{r_{0}e^{t}}{\sqrt{1+r_{0}(e^{2t}-1)}}

Again, this should be fairly simple, but I'm stuck.. Please help me!
 
Physics news on Phys.org
∫ (\frac{1}{r} + \frac{1}{1-r^2})dr = t is false. r is missing on the second fraction.
 
Thank you! Got it now :)
 
Paalfaal said:
Ok, this should be quite simple. I've been looking at this problem for quite some time now, and I'm tired.. Please help me!

The equation to solve is r' = r(1-r^2)

The obvious thing to do is to do partial fraction expansion and integrate(from r_{0} to r):

∫ (\frac{1}{r} + \frac{1}{1-r^2})dr = t
That's an incorrect "partial fraction expansion". 1- x^2= (1- x)(1+ x) so
\frac{1}{r(1- r^2)}= \frac{1}{r}- \frac{1}{r- 1}- \frac{1}{r+ 1}
There is no "trig substitution" involved in integrating that. the integral is
ln(r)- ln(r-1)- ln(r+1)+ C= ln\left(\frac{r}{r^2-1}\right)+ C

After some trig substitution: ln[r\sqrt{\frac{r+1}{1-r}}] = t (evaluated from r_{0} to r). Then take the exponential on both sides.

r\sqrt{\frac{r+1}{1-r}} = r_{0}\sqrt{\frac{r_{0}+1}{1-r_{0}}}e^{t}

This leads to kind of a nasty expression which can't be solved explicitly for r(t) (at least it seems that it can't be solved explicitly).. So this is where I'm stuck.

The solution to this ploblem is: r(t) = \frac{r_{0}e^{t}}{\sqrt{1+r_{0}(e^{2t}-1)}}

Again, this should be fairly simple, but I'm stuck.. Please help me!
You have integrated incorrectly.
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...

Similar threads

Back
Top