I am having some troubles to get the logistic equation from its
differential form. For example:
sage: tau = var('tau') # proportionality constant
sage: k = var('k') # max population size
sage: t = var('t') # independent variable
sage: p = function('p', t) # dependent variable
sage: myode =
The problem is that the solution that desolve returns:
sage: myode = tau*diff(p,t) == p*(1-p/k)
sage: sol(t) = desolve(de=myode, ivar=t, dvar=p)
sage returns:
-tau*log(-k + p(t)) + tau*log(p(t)) == c + t
isn't in an easy enough form for plot to display. You're setting
sol(t) not to an