Sorry for spamming,
but I figured my initial code didn't make sense in the way mu was
computed; so I changed it to the following. The error message stays
the same, though.
====================
t = var ('t')
S_0 = 1.5
X_0 = 0.05
Y_XS = 0.5
K_S = 0.007
mu_max = 0.8
X = function ('X', t)
S = function ('S', X)
def mu (S):
return (mu_max * S) / (K_S + S)
dXdt = diff (X, t) == mu(S) * X
dSdt = diff (S, t) == -mu(S) * X / Y_XS
desolve_system ([dXdt, dSdt], [X, S], ics = [0, X_0, S_0])
====================
Richard
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
