Re: ODE system - shooting for a terminal point
- To: mathgroup at smc.vnet.net
- Subject: [mg121178] Re: ODE system - shooting for a terminal point
- From: Oliver Ruebenkoenig <ruebenko at wolfram.com>
- Date: Sat, 3 Sep 2011 08:03:26 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201108311004.GAA12094@smc.vnet.net>
On Wed, 31 Aug 2011, The Mighty Algernon wrote:
> Hi,
>
> I wonder if you could help me.
>
> I have a system of ODE in which I have to figure out the initial
> conditions of part of the differential equations. I managed to do this
> by the shooting method by writing the system as a function of the
> initial conditions.
>
> For example solution[yini_,zini_]:=NDSolve[......
>
> then I used FindRoot on the set of terminal restrictions to find yini
> and zini (these are the initial conditions). This procedure worked
> although it is terribly sensitive to the first guesses on the FindRoot
> and some results I get are obviously wrong. But this is something I
> can refine.
>
> Now I have a similar goal but the variable I'm trying to guess is not
> an initial condition of a differential equation but the terminal value
> of the range of the exogenous variable. For example,
>
> If I have a system of 2 ODE , y[x] and z[x], I want to know for {x,
> 0,xtarget} what is my xtarget. Think of x as time and I want to figure
> when is the terminal time of the system such that some restrictions at
> the terminal time are met.
Hi,
perhaps this could be done as an Event?
Oliver
>
> My first ideia was to adapt the above procedure, writing
>
> solution[xtarget_]:=NDSolve[......
> then giving the set of restrictions and then using FindRoot to get
> xtarget, but Mathematica gives me an error message the xtarget is not
> a real value.
>
>
> My question is :
>
> Is it possible to adapt the procedure outlined above, meaning that
> I'm just doing some kind of error here (that I'll have to correct), or
> Mathematica does not handle this kind of problems directly in NDSolve
> and I have to write some sort of loop with ever more refined guesses
> for xtarget?
>
> My best regards
>
> Joao Pereira
>
>
>
>
>
>
>
>