This notebook provides a quick test for bifurcation plots in MuMoT.¶

import mumot

mumot.__version__

%%model
$
U -> A : g_A
U -> B : g_B
A -> U : a_A
B -> U : a_B
A + U -> A + A : r_A
B + U -> B + B : r_B
A + B -> A + U : s
A + B -> B + U : s
$

model1 = mumot.parseModel(In[2])

model2 = model1.substitute('U = N - A - B')

model3 = model2.substitute('a_A = 1/v_A, a_B = 1/v_B, g_A = v_A, g_B = v_B, r_A = v_A, r_B = v_B')

model4 = model3.substitute('v_A = \mu + \Delta/2, v_B = \mu - \Delta/2')

bifurcation1 = model4.bifurcation('s','A-B', choose_xrange = [0,5], initWidgets={'initBifParam': [4.5,0,5,0.1]})
# to get started, here is a parameter combination that works
# set \mu = 3
# set \Delta = 0.1 for unfolding of pitchfork
# set \Delta = 0 for pitchfork

bifurcation1.showLogs()

bifurcation2 = model4.bifurcation('\\Delta','A-B', initBifParam = 2.0, params = [('\\mu', 4), ('s', 4)], 
                                  choose_xrange = [-2,2])

bifurcation2.showLogs()

multiController4 = mumot.MuMoTmultiController([model4.bifurcation('s','A', silent = True), 
                                             model4.bifurcation('s','B', silent = True)], 
                                            shareAxes = False)

multiController4.showLogs()