MathGroup Archive 2004

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

Search the Archive

Re: populating a list from system of odes


Hi, you mean


eqn = {a'[t] == k7 - k5 a[t] + k2 ab[t] - k1 a[t] b[t], 
    ab'[t] == -k2 ab[t] + k1 a[t] b[t], 
    b'[t] == k8 + k2 ab[t] - k6 b[t] - k1 a[t] b[t], 
    c'[t] == k3 ab[t] - k4 c[t]};

Union[Cases[eqn, _[t], Infinity]] /. a_Derivative[_][_] :> Sequence[]


Regards
  Jens

sean kim wrote:
> 
> hello group.
> 
> thanks all for the replies. they have been tremendiously helpful.
> 
> (about rossler vz lorenz. you are right, jens. what i did was to copy
> the first system that had attractor from the help menu.  lol)
> 
> let's say you have a system of ode's like..
> 
> {
> a'[t] == k7 - k5 a[t] + k2 ab[t] - k1 a[t] b[t],
> ab'[t] == -k2 ab[t] + k1 a[t] b[t],
> b'[t] == k8 + k2 ab[t] - k6 b[t] - k1 a[t] b[t],
> c'[t] == k3 ab[t] - k4 c[t]
> }
> 
> now what I want to do is get a list of the  variables in the system
> without having to type them out.
> 
> so that I would get {a[t],ab[t], b[t], c[t]}
> 
> as you can see from my previous posts, the problems i'm having are
> very similar.  I'm sur eit has to do with no tunderstanding the
> fundamentals of mathematica... but meanwhile, how would you guys
> approach that problem above?
> 
> once i have that, i can change the system to steady state easily with
> 
> % /. _'[t] -> 0
> 
> which gives
> 
> {
> 0 == k7 - k5 a[t] + k2 ab[t] - k1 a[t] b[t],
> 0 == -k2 ab[t] + k1 a[t] b[t],
> 0 == k8 + k2 ab[t] - k6 b[t] - k1 a[t] b[t],
> 0 == k3 ab[t] - k4 c[t]
> }
> 
> which i can use with the list i'm tryign to populate in Solve[]
> routine which then can be used in NDSolve[] with the codes you guys
> explained.
> 
> thank you all very much in advance.
> 
> sean


  • Prev by Date: one liner for a function?
  • Next by Date: Re: Mathematica and XML
  • Previous by thread: Re: populating a list from system of odes
  • Next by thread: Mathematica and XML