Re: Solving DEs by NDSolve
- To: mathgroup at smc.vnet.net
- Subject: [mg41052] Re: Solving DEs by NDSolve
- From: Selwyn Hollis <hollisse at mail.armstrong.edu>
- Date: Wed, 30 Apr 2003 04:23:04 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Payal,
Look at your third equation. There's a lower-case a that should be an
upper-case A.
(Don't feel bad... it was not at all easy to find!)
Cheers,
Selwyn Hollis
-----------------------------------------------------------
Sheth, Payal R. wrote:
Hi,
I have been having consistent problems applying the
NDSolve function to a set of differential equations.
I keep getting the error stating
NDSolve::"ndnum": "Encountered non-numerical value for a derivative at \
\!\(t\) == \!\(1.9000180679486`*^-311\)."
I do not have a training in differential equations as such, so it is
very
difficult for me to nail down the problem. Any help would be
extremely appreciated.
NDSolve[{A'[t] == -k1*e[t]*l[t] +
k11*c[t] + k2*e[t] - k22*A[t] - ((vmaxe*A[t])/(Kme + A[t])) +
kcat*c[t],
e'[t] == +((vmaxe*A[t])/(Kme + A[t])) - k2*e[t] + k22*A[t],
l'[t] == -k1*a[t]*l[t] + k11*c[t] + ((vmaxl*b[t])/(Kml + b[t])),
b'[t] == kcat*c[t] - ((vmaxl*b[t])/(Kml + b[t])),
c'[t] == k1*A[t]*l[t] - k11*c[t] - kcat*c[t],
A[0] == 0.00, e[
0] == 100, l[0] == 100, b[0] == 0, c[0] == 0}, {A, e, l,
b, c}, {t, 0, 100}]
k1 = 0.09;
k11 = 0.6;
kcat = 25;
k2 = 1;
k22 = 0.01;
vmaxe = 450;
Kme = 50;
vmaxl = 1;
Kml = 100;
Best Regards,
Payal Sheth