Coupled ODEs: not solving [REPOST]
- To: mathgroup at smc.vnet.net
- Subject: [mg6519] Coupled ODEs: not solving [REPOST]
- From: Patrick Jemmer <padz at joule.pcl.ox.ac.uk>
- Date: Sun, 30 Mar 1997 22:23:51 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
-- Hello all:
I was wondering if anyone might be able to help
me with my attempts to get MMA to solve the attacked set
of ODEs...
There are 8 equations, 1 for each unknown...
Am I simply hitting the ceiling in the complexity of such
systems MMA can handle ?
Cheers for any assistance.
Patrick.
DSolve[
{Derivative[1][En][t] ==
km4*EI[t] - k3*En[t] + km1*ES[t] + km3*Inh[t] - k4*En[t]*Inh[t] -
km2*En[t]*P[t] - k1*En[t]*S[t],
Derivative[1][S][t] == km1*ES[t] + k6*PI[t] - k1*En[t]*S[t] -
km6*Inh[t]*S[t],
Derivative[1][ES][t] ==
-(k2*ES[t]) - km1*ES[t] + km2*En[t]*P[t] + k1*En[t]*S[t],
Derivative[1][P][t] == k2*ES[t] - k7*P[t] - km2*En[t]*P[t] -
k5*Inh[t]*P[t] + km5*PI[t] + km7*T[t],
Derivative[1][Inh][t] == km4*EI[t] + k3*En[t] - km3*Inh[t] -
k4*En[t]*Inh[t] - k5*Inh[t]*P[t] + k6*PI[t] + km5*PI[t] -
km6*Inh[t]*S[t],
Derivative[1][EI][t] ==
-(km4*EI[t]) + k4*En[t]*Inh[t],
Derivative[1][PI][t] == k5*Inh[t]*P[t] - k6*PI[t] - km5*PI[t] +
km6*Inh[t]*S[t],
Derivative[1][T][t] == k7*P[t] - km7*T[t],
En[0] == 1, S[0] == 4, ES[0] == 0, P[0] == 0, Inh[0] == 0, EI[0] == 0,
PI[0] == 0, T[0] == 1},
{En[t], S[t], ES[t], P[t], Inh[t], EI[t], PI[t], T[t]},t]
_
_ __ __ _ __| |_ __ ___ __ _ _ __
| '_ \ / _` |/ _` | '_ ` _ \ / _` | '_ \
| |_) | (_| | (_| | | | | | | (_| | | | |
| .__/ \__,_|\__,_|_| |_| |_|\__,_|_| |_|
|_|
--------------------------------------------------------------
| Patrick Jemmer, | http://berne.pcl.ox.ac.uk/~padz |
| Physical & Theoretical | e-mail: padz at joule.pcl.ox.ac.uk |
| Chemistry Laboratory, | tel: +44-1865-2-75161 (work) |
| South Parks Road, | fax: +44-1865-2-75410 |
| Oxford OX1 3QZ | http://joule.pcl.ox.ac.uk/ (PTCL)|
--------------------------------------------------------------