def runme(writefilefunc, n=4, y0=[0.75, 0.75, 1.25], t0=0, t1=10, dt=0.01):
writefilefunc()
times, funcvals = ihf.simulateHillModel(n, y0, t0, t1, dt)
ihf.plotResults(times, funcvals, {
'linewidth': 4,
'label': [r'$x_1$', r'$x_2$', r'$x_3$']
}, {'ncol': 3})
def runme(RHS=RHS_5D_2015_10_21,
t1=200,
hillexp=10,
plotoptions={},
legendoptions={},
figuresize=(15, 10)):
y0 = np.array([1, 1, 1, 1, 1])
t, y = simulate(RHS, y0, t1=t1, hillexp=hillexp)
ihf.plotResults(t, y, plotoptions, legendoptions, figuresize)
def bistablecheck(n=7, y0=[0.75, 0.75, 1.25], t0=0, t1=10, dt=0.01):
eqns = [
lambda X: -X[0] + chf.negHillFunction(2, 1, 6, n, X[
1]) * chf.negHillFunction(3, 1, 2, n, X[2]),
lambda X: -X[1] + chf.negHillFunction(5, 1, 4, n, X[0]),
lambda X: -X[2] + chf.negHillFunction(4, 1, 3, n, X[1])
]
times, funcvals = ihf.plainIntegrate(eqns, y0, t0, t1, dt)
ihf.plotResults(times, funcvals, {
'linewidth': 4,
'label': [r'$x_1$', r'$x_2$', r'$x_3$']
}, {'ncol': 3})
def compare(test, n, y0, t0=0, t1=10, dt=0.01):
times, funcvals = ihf.simulateHillModel(n, y0, t0, t1, dt)
ihf.plotResults(times, funcvals)
times, funcvals = test.evalAnalytic(n, y0, t0, t1, dt)
ihf.plotResults(times, funcvals)
def run5pt1():
y0=np.array([198,162,321,503,638])
# y0=np.array([1,1,1,1,1])
t,y=simulate(RHS5pt1,y0)
ihf.plotResults(t,y)