Add example of how to compute nullclines using bifurcation kit #1194
Conversation
TorkelE
changed the title
docs/src/steady_state_functionality/examples/nullcline_plotting.md
Outdated
Show resolved
Hide resolved
| using Catalyst | ||
| t = default_t() | ||
| @parameters v K n | ||
| @species X(t) Y(t) | ||
| rxs = [ | ||
| Reaction(hillr(Y, v, K, n), [], [X]), | ||
| Reaction(hillr(X, v, K, n), [], [Y]), | ||
| Reaction(1.0, [X], []), | ||
| Reaction(1.0, [Y], []), | ||
| ] | ||
| @named bs_switch = ReactionSystem(rxs, t) | ||
| bs_switch = complete(bs_switch) |
There was a problem hiding this comment.
Why create it programmatically? Seems like there is no advantage in this case (i.e. we aren't showing off any programmatic features)?
There was a problem hiding this comment.
I think because the second half of the tutorial was programmatic I figured it was nice/easier to do a full programmatic one. However, you might be right that it doesn't really add anything.
…g.md Co-authored-by: Sam Isaacson <[email protected]>
|
You might be right about extracting functions. I had some idea why this made sense (for more complicated curves, although here I used a simpler case because I am familiar with it). You might be right that we really can do this generally without involving BifurcationKit. Maybe I will go back to the drawingboard a bit and think if I can streamline things. |
WIP