MathGroup Archive 2003

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

Search the Archive

Re: Solving DEs by NDSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg41035] Re: Solving DEs by NDSolve
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Wed, 30 Apr 2003 04:19:09 -0400 (EDT)
  • Organization: Universitaet Leipzig
  • References: <b8lgh1$99s$1@smc.vnet.net>
  • Reply-to: kuska at informatik.uni-leipzig.de
  • Sender: owner-wri-mathgroup at wolfram.com

Hi,

and what is a[t]  in

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])),
                   ^^^^
                   here ??

>     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}]

Regards
  Jens

"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


  • Prev by Date: Re: Show[] applied to a Plot[] and a ListPlot[] ??
  • Next by Date: RE: eqiprobable intervals
  • Previous by thread: Re: Solving DEs by NDSolve
  • Next by thread: Writing Applications for Mac OS X