Re: populating a list from system of odes
- To: mathgroup at smc.vnet.net
- Subject: [mg45933] Re: populating a list from system of odes
- From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
- Date: Fri, 30 Jan 2004 04:16:00 -0500 (EST)
- Organization: Universitaet Leipzig
- References: <bvaq98$ah0$1@smc.vnet.net>
- Reply-to: kuska at informatik.uni-leipzig.de
- Sender: owner-wri-mathgroup at wolfram.com
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