MathGroup Archive 1998

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Nonlinear Differential Equation


Hi David,
   Interesting problem; I evaluated the DSolve below and got an
expression with two constants c[1] and c[2]. Applying your boundary
conditions (manually) I think the constants are c[1]=1/4-I Pi/2 and
c[2]=0. Then by playing around with the integral, I'll guess that the
exact amplitude and period are given by the expressions A and P below.
(Technically, you have to take eps->0 in the expression for P). You'll
see these numbers are pretty close to the ones you quoted. Hope this
helps.

DSolve[w''[x]+2 w[x]/(1-w[x]^2)==0,w[x],x]

A=Sqrt[1-E^(-1/2)];
P[eps_]:= N[2 Sqrt[2]NIntegrate[(1/2+ Log[1-t^2])^(-1/2),{t,0,A-eps}]];


N[A]
0.627271

P[10^(-7)]
3.70814

Alan Lewis

David Djajaputra wrote:
> 
> Hello,
> 
> I'm just wondering if anyone knows how to solve this DE:
> 
> y''(x) = - 2 y(x) / [1 - y^2(x)],   y(0) = 0,  y'(0) = 1.
> 
> The function is periodic. Mathematica says that it's transcendental but
> it cannot solve it. I'm curious if the amplitude and/or the period can
> be obtained in closed forms. By the way, if you plot it the amplitude
> is
> 
> about 0.6 and the period is about 3.7.
> 
> Much thanks in advance,
> 
> David


  • Prev by Date: Re: Q: Setting up Mathematica
  • Next by Date: Outputting Strings to a file
  • Previous by thread: Nonlinear Differential Equation
  • Next by thread: Problem with WMF Part 2