MathGroup Archive 2009

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

Search the Archive

Please help debug this set of simple diff eqs- having trouble with

  • To: mathgroup at smc.vnet.net
  • Subject: [mg101944] Please help debug this set of simple diff eqs- having trouble with
  • From: angela <amb458 at cornell.edu>
  • Date: Fri, 24 Jul 2009 06:10:36 -0400 (EDT)

mechanism = NDSolve[{
    s'[t] == -k1f  s[t] b2[t] + k1r c[t],
    c'[t] == k1f [a] [s] - k1r c[t],
    p'[t] == k2[c] == bp[t] == b'[t],
    a'[t] == -k1f a[t] s[t] + k1r c[t],
    b2'[t] == -k1f s[t] b2[t] + k2f c[t],
    s[0] == 0.005;
    b[0] == 0.025;
    p[0] == 0;
    bp[0] == 0;
    b[0] == 0}
   /. {k1f -> 0.1, k1r -> 0.01, k2 -> 1},
  {b2, s, c, p, b, bp},
  {t, 0, 1500}]


NDSolve::deqn: Equation or list of equations expected instead of True \
in the first argument {s^\[Prime][t]==0.01 c[t]-0.1 b2[t] \
s[t],<<4>>,True}. >>


  • Prev by Date: Re: documentation bug?
  • Next by Date: Re: Select any
  • Previous by thread: Re: Forcing vertical scale on Histogram
  • Next by thread: Graph issue