Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2009

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

Search the Archive

Re: perturbation methods example from stephen lynch's book?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg100644] Re: perturbation methods example from stephen lynch's book?
  • From: sean k <sean_incali at yahoo.com>
  • Date: Wed, 10 Jun 2009 05:33:36 -0400 (EDT)
  • Reply-to: sean_incali at yahoo.com

I just realized that the codes I posted isn't complete.

Not sure what Lynch is doing in his notebook, but his Collect[] line doesn't seem to work. It's probably because he's using Dt[x, {t, 2}]  on just x not x[t]


epsilon="epsilon";
SetAttributes[{w1,epsilon},Constant]

x=x0+epsilon*x1+epsilon^2*x2;

Collect[(1+2 epsilon w1+epsilon^2 (w1^2+2 w2))Dt[x,{t,2}]+x-epsilon x^3,epsilon]

DSolve[{x0''[t]+x0[t]==0,x0[0]==1,x0'[0]==0},x0[t],t]

DSolve[{x1''[t]+x1[t]==Cos[t]^3+2 w1 Cos[t],x1[0]==0,x1'[0]==0},x1[t],t]
Simplify[%]




--- On Tue, 6/9/09, sean_incali at yahoo.com <sean_incali at yahoo.com> wrote:

> From: sean_incali at yahoo.com <sean_incali at yahoo.com>
> Subject: perturbation methods example from stephen lynch's book?
> To: sean_incali at yahoo.com
> Date: Tuesday, June 9, 2009, 4:57 PM
> Hello group,
> 
> This message is a it long.
> 
> I was reading Stephen Lynch's "dynamical systems with
> applications
> using mathematica" and noticed he's using asymptotics and
> perturbation
> methods in chapter 4, section 4 perturbation methods.
> 
> Except he's only showing the code for linstedt-poincare
> methods which
> fails for the example given in the book. (van der pol
> equation)
> 
> (*See Example 8:The Lindstedt-Poincare technique.*)
> SetAttributes[{w1, epsilon}, Constant]
> 
> x = x0 + epsilon*x1 + epsilonË?2*x2;
> 
> Collect[(1 + 2 epsilon w1 + epsilonË?2 (w1Ë?2 + 2 w2))
> Dt[x, {t, 2}] + x
> - epsilon xË?3, epsilon];
> 
> (*The O (1) equation.*)
> DSolve[{x0''[t] + x0[t] == 0, x0[0] == 1, x0'[0] == 0},
> x0[t], t]
> 
> (*The O (epsilon) equation.*)
> DSolve[{x1''[t] + x1[t] == Cos[t] Ë?3 + 2 w1 Cos[t], x1[0]
> == 0, x1'[0]
> == 0}, x1[t], t];
> Simplify[%]
> 
> 
> 
> He then discusses method of multiple scales to approximate
> the
> solutions to van der pol equation (x'' + x =  epsilon
> (1 - x^2) x',
> given x=a, and x'=0)
> 
> I was wondering if anyone worked out that example. (example
> 10 in
> chapter 4). If so, can someone kindly share the codes for
> it?
> ie.
> 1. solution to PDEs, O(1) and O(epsilon), that results from
> changing
> the time scales,
> 2. using TrigReduce to simplify,
> 3. remove the secular terms,
> 4. impose ICs and approximate the one term O(epsilon)
> solution, x_ms
> 
> I might be asking a lot, but I was hoping someone has the
> codes for
> it.
> 
> Thanks in advance
> 
> Sean
> 
> 
> 
> 


  • Prev by Date: Slow/jerky animations inside manipulate (more details)
  • Next by Date: Re: mathematica tutor for NYC high school student
  • Previous by thread: perturbation methods example from stephen lynch's book?
  • Next by thread: Re: perturbation methods example from stephen lynch's book?