MathGroup Archive 2009

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

Search the Archive

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


Hi Simon,

Thanks for the reply.

Actually I've found some resources on wolfram website.

http://library.wolfram.com/infocenter/MathSource/5269/

It seems like an excellent resource for those who wants to learn and
use asymptotics in mathematica. Only problem is that it uses
mathematica 3...

And Lynch's example 7 and typos come from his duffing equation.
Instead of x''+x+eps x^3 =0, he uses x''+x =eps x^3

Differences show in O(e) eqn.

o1 = DSolve[{x0[t] + x0''[t] == 0, x0'[0] == 0, x0[0] == 1 }, x0[t], t]
[[1, 1]]
DSolve[{x0[t]^3 + x1[t] + x1''[t] == 0 /. {x0[t] -> Cos[t]}, x1'[0] ==
0, x1[0] == 0 }, x1[t], t][[1, 1]]
oe = % // Simplify // Expand
"x~xp=" o1[[2]] + \[Epsilon] (oe)[[2]]
% // TraditionalForm


vs.

o1 = DSolve[{x0[t] + x0''[t] == 0, x0'[0] == 0, x0[0] == 1 }, x0[t], t]
[[1, 1]]
DSolve[{-x0[t]^3 + x1[t] + x1''[t] == 0 /. {x0[t] -> Cos[t]}, x1'[0]
== 0, x1[0] == 0 }, x1[t], t][[1, 1]]
oe = % // Simplify // Expand
"x~xp=" o1[[2]] + \[Epsilon] (oe)[[2]]
% // TraditionalForm


  • Prev by Date: Re: Creating Matrix from vectors specific issue.
  • Next by Date: Re: Correction to "Fundamental Theorem of Calculus and Mathematica"
  • Previous by thread: Re: perturbation methods example from stephen lynch's book?
  • Next by thread: Graphics3D not visible under X11 (was: Compositing, 3D graphics, and KDE 4.2)