Ask Your Question

Epidemiology equilibria solved, but what happened??

asked 2016-11-20 02:59:00 -0500

kaylavb23 gravatar image

var('beta gamma S I R N')

dS = - beta * S * I / N dI = beta * S * I / N - gamma * I dR = gamma * I

solve([dS == 0, dI == 0, dR == 0], [S, I, R]) [[S == r7, I == 0, R == r8]]

Why am i getting different values for S and R which correlate to how many times I run the code? At first I thought that it was because I did an interpolate above, but when I went and modified the code a bit at the beginning I was still left with increasing equilibria values.

the next run will be [[S == r9, I == 0, R == r10]] Does anyone know why this is happening??

edit retag flag offensive close merge delete

1 answer

Sort by ยป oldest newest most voted

answered 2016-11-20 06:05:18 -0500

tmonteil gravatar image

updated 2016-11-20 06:14:23 -0500

The solution [[S == r7, I == 0, R == r8]] should be read as follows: the set of solutions is an infinite two-parameter space, which can be parametrized by two anonymous parameters r7 and r8 that Sage created for the purpose of describing the solution set. Each time you solve a new equation whose set of solutions is huge, Sage will create new parameters (symbolic variables) to serve as a parameters to describe the solutions set.

edit flag offensive delete link more

Your Answer

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

Add Answer

Question Tools


Asked: 2016-11-20 02:59:00 -0500

Seen: 41 times

Last updated: Nov 20 '16