Timeline for (In)Stability of equilibrium points in a vector plot (continued)?
Current License: CC BY-SA 4.0
19 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Sep 28, 2021 at 11:18 | comment | added | Math | When I tried performing Routh-Hurwitz criteria manually, the third condition fails, can you help? I will post a new question. | |
| Sep 22, 2021 at 11:23 | comment | added | Math | Thats what i used, i got the error "BoxRatios::nlist3: 1 is not a list of three numbers." | |
| Sep 21, 2021 at 15:27 | comment | added | Chris K |
The Fixed option in PlotEcoPhasePlane fixes one of the variables, taking a 2D slice through the 3D phase space.
|
|
| Sep 21, 2021 at 14:34 | comment | added | Math | When you plotted your 2d slices earlier, I didn't see the red highlight albeit we used the same code, why? | |
| Sep 21, 2021 at 14:22 | comment | added | Math | Having just skimmed through, this is fantastic! I will have a deeper look when I get some time. Thank you, this is greatly appreciated dear friend! I may have some additional questions aha | |
| Sep 21, 2021 at 14:15 | comment | added | Chris K |
1) Use Color in SetModel. 2) EcoPhasePlane plots 2D isoclines, which is not appropriate for your 3D model without fixing one variable. Use RuleListPlot to plot dynamics in this 2D slice instead. See my addition for details. BTW, a tip: try to use the same symbols in your different codes as in the paper - makes debugging easier!
|
|
| Sep 21, 2021 at 14:13 | history | edited | Chris K | CC BY-SA 4.0 |
added 1604 characters in body
|
| Sep 21, 2021 at 13:39 | comment | added | Math | No problem, I have figured it! two last questions: 1) is there a way to change the colour for $e$ and $i$ in the phase plots? I wish $e$ to be orange and $i$ to be red to match my code. 2) why is the red highlight around the phase plot? | |
| Sep 21, 2021 at 12:18 | comment | added | Chris K | I have a feeling that "help me recreate this paper" would be seen as outside the scope of this site. The usual emphasis is on more narrow, technical questions related to Mathematica, particularly ones that might help others. | |
| Sep 21, 2021 at 11:46 | comment | added | Math | To put it simply, I want to recreate what the paper has done for my model(using your eco package, isoclines are shown). Maybe I should post a separate question? | |
| Sep 20, 2021 at 16:51 | comment | added | Chris K |
Your system is 4D, reduced to 3D by the constraint $S+E+I+R=1$. Therefore the isoclines will be surfaces in 3D space as I've shown above. You could also take 2D slices as I've shown above, but it's not showing the full picture. To plot the dynamics of one variable vs another, use ParametericPlot. It's not clear what else needs to be answered here.
|
|
| Sep 20, 2021 at 15:56 | comment | added | Math | I want the isoclines as all models I have considered, I have used your package. I tried something similar and I can get the solutions to the ODE system but I cant get the phase plots. I will attach the solution plot soon | |
| Sep 20, 2021 at 15:50 | comment | added | Chris K |
Based on the paper, I assume you don't want isoclines or streams, just plotting a single trajectory. Since r is only defined algebraically, the easiest thing to do is ParametricPlot[ Evaluate[{s[t], 1 - s[t] - e[t] - i[t]} /. sol], {t, 0, 100}, AspectRatio -> 1]. For variables that are listed as Components you could use RuleListPlot[sol, {s, e}], so another approach would be to explicitly model r as well.
|
|
| Sep 20, 2021 at 15:32 | comment | added | Chris K | BTW you need to use $$ to mark LaTeX code | |
| Sep 20, 2021 at 15:07 | comment | added | Math | It will give me coding errors again for some reason but I have attached it. The main interest is the solution to the ODE system and the phase plots. | |
| Sep 20, 2021 at 14:56 | comment | added | Chris K | Sure, attach the paper and I'll take a look when I have time. | |
| Sep 20, 2021 at 14:48 | comment | added | Math | Nice answer yet again! Is there a way to plot say $s$ vs $r$? I am reading a paper(their model is a little more complex) but very similar. I can attach it in the question as an edit if you'd like to see. They have managed to produce results that I am unable to produce(albeit quite close), maybe you can try? | |
| Sep 19, 2021 at 14:37 | history | edited | Chris K | CC BY-SA 4.0 |
edited body
|
| Sep 19, 2021 at 14:30 | history | answered | Chris K | CC BY-SA 4.0 |