MathGroup Archive 2001

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

Search the Archive

Re: Modified shooting method

  • To: mathgroup at
  • Subject: [mg30511] Re: Modified shooting method
  • From: Jens-Peer Kuska <kuska at>
  • Date: Fri, 24 Aug 2001 20:58:07 -0400 (EDT)
  • Organization: Universitaet Leipzig
  • References: <9m522e$qdc$>
  • Sender: owner-wri-mathgroup at


eigenvalue problems work not very well with simple shooting
The problem ist, that for a fixed value w[i] you have 
 - no solution if w[i] is not a eigenvalue
 - an infinite number of solutions if w[i] *is* an eigenvalue

The bes solution is to use a double or multiple
shooting method.  You should have a look at

file    mus.doc  mus.doc plus dependencies
by      Mattheij and Staarink
for     ordinary differential equation boundary-value problem solver
alg     multiple shooting
,       whence you may also fine DEQNS
prec    double
gams    I1b1 I1b2

An better way may be to use a matrix approach and solve 
the matrix eigenvalue problem. The simplest idea is to
us finite differences 
y''[x] -> (y[i-1]-2 y[i]+y[i+1])/h^2

for the inner points, and solve
the problem


Or to expand the solution in a orthogonal basis you like,
and solve the eigenvalueproblem
Hope that helps

Christophe Le Poncin-Lafitte wrote:
> Hello,
> Sorry for my english, which is poor.
> I'm working about a modelization of cyclonic and anticyclonic structures in
> geophysic ; so this is a study of the nonlinear adjustment of a density
> front.
> During our experiment, we observe the adjustment of a cylindrical structure,
> which becomes a little stable lens (during 40s), before the development of
> instabilities. But during the adjustment, we observe the development of
> inertial-gravity waves.
> Now we search to modelize our system.
> Our goal is to quantify the energy of this waves, to evaluate the importance
> of these for the comprehension of the adjustment first, then to understand
> the relative stability of the system, and finally to find why we observe
> instabilities, which destroy the system.
> I have a little problem to solve : a boundary value problem.
> My ODE is quite difficult :
> y"[x]+q[w,x]*y[x]==0
> where
> q[w,x]=w^2/(1-BesselI[0,x]/BesselI[0,rap])-1-3/4/r^2-BesselI[1,r]/Bessel[0,r

> rap is a parameter, which is a constant. w is the eigenmodes of
> inertial-gravity waves.
> My boundary conditions are  : y[0]=0 and y[rap]=0
> I take for initial condition on y' : y'[0]=1.
> In fact, I have to done my shooting on w, to isolate the different
> eigenmodes of inertial-gravity waves. So, this is a modified shooting
> method.
> I want to know if anybody has already treat this type of problem ?
> Regards,
> Christophe

  • Prev by Date: Re: Plot Joined Intelligently
  • Next by Date: RE:Controlling evaluation in Symbolize
  • Previous by thread: Modified shooting method
  • Next by thread: Solve system of ODEs