Ask Your Question

Epidemiology equilibria solved, but what happened??

asked 2016-11-20 09:59:00 +0200

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 13:05:18 +0200

tmonteil gravatar image

updated 2016-11-20 13:14:23 +0200

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 09:59:00 +0200

Seen: 259 times

Last updated: Nov 20 '16