Ask Your Question

using logistic differential equatin (reproducing Ti Nspire CAS code)

asked 2016-07-29 02:43:11 -0500

etb gravatar image

updated 2016-07-29 02:55:38 -0500

I’m trying to reproduce some results from Ti Nspire CAS to use a differential equation. Only secondarily am I interested to actually solve the differential equation.

My Ti Nspire codes is as follows

deSolve(y’ = 2*10^(-5)*y(1500-y) and y(0) = 50, t, y)

and yields the following result

y = 1500*e^((3*t)/100) / (e^((3*t)/100) + 29)

image description

Can anyone give me a hit as to how I can do this in Sagemath?

edit retag flag offensive close merge delete

1 answer

Sort by » oldest newest most voted

answered 2016-07-29 03:56:16 -0500

mforets gravatar image

updated 2016-07-29 03:57:04 -0500

t = var('t')
y = function('y')(t)
ye = desolve(diff(y,t) == 2*10^(-5)*y*(1500-y), y, ics=[0,50])

ye = ye*3/100
yt = solve(ye.simplify_log(), y)

$$ \newcommand{\Bold}[1]{\mathbf{#1}}\left[y\left(t\right) = \frac{1500 e^{\left(\frac{3}{100} t\right)}}{e^{\left(\frac{3}{100} t\right)} + 29}\right] $$

Remark. Here I multiplied by $3/100$, since otherwise the simplify_log() was not factoring the terms containing $y(t)$. For instance,

f = 2*log(x) - 2*log(x-1)

gives $\log(x^2/(x-1)^2)$, but

f = 1/3*log(x) - 1/3*log(x-1)

gives $-\frac{1}{3}\log(x-1) + \frac{1}{3}\log(x)$, instead of $\frac{1}{3}\log(x/(x-1))$ (and in the first case I would prefer to get rather $2\log(x/(x-1))$; neither factor() nor simplify_full() do it -- why?).

edit flag offensive delete link more


Thanks a lot. I cannot answer your question, but I'll get back to you here if I find an answer to that.

etb gravatar imageetb ( 2016-07-29 04:01:44 -0500 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower


Asked: 2016-07-29 02:43:11 -0500

Seen: 53 times

Last updated: Jul 29 '16