sim.add(" x'=r2*x*(1-x/K2) ", 1, plot=1) sim.params(r=2, K=50, r2=3, K2=50) sim.run(0, 20) # ## adding some data # In[31]: t_data = array([7, 14, 21, 28, 35, 42, 49, 56, 63, 70, 77, 84], float) h_data = array([ 17.93, 36.36, 67.76, 98.10, 131, 169.5, 205.5, 228.3, 247.1, 250.5, 253.8, 254.5 ]) sim = Simulation() sim.figsize = (6, 4) sim.add(" y'=r*y*(1-y/K) ", 1, plot=True) sim.params(r=.2, K=270) sim.add_data(t=t_data, y=h_data, plot=True) sim.run(0, 100) # In[35]: model = MCMCModel(sim, r=Uniform(0, 1), K=Uniform(100, 400)) # In[36]: model.run_mcmc(500) # In[37]:
# <codecell> x,y=sim.t,sim.mice plot(x,y,'r--') xlabel('Days') ylabel('Number of Mice') # <markdowncell> # ## Predator-Prey Dynamics # <codecell> sim=Simulation() sim.figsize=(20,5) sim.add("deer' = r*deer*(1-deer/K)-c*deer*wolf", initial_value=350, plot=True) sim.add("wolf' = -Wd*wolf+D*deer*wolf", initial_value=50, plot=True) sim.params(r=0.25,D=0.001,c=0.005,Wd=0.3,K=1e500) print sim.equations() sim.run(0,500) # fig=sim.figures[0]
from pylab import plot, xlabel, ylabel # In[6]: x, y = sim.t, sim.mice plot(x, y, 'r--') xlabel('Days') ylabel('Number of Mice') # ## Predator-Prey Dynamics # In[8]: sim = Simulation() sim.figsize = (20, 5) sim.add("deer' = r*deer*(1-deer/K)-c*deer*wolf", initial_value=350, plot=True) sim.add("wolf' = -Wd*wolf+D*deer*wolf", initial_value=50, plot=True) sim.params(r=0.25, D=0.001, c=0.005, Wd=0.3, K=1e500) print(sim.equations()) sim.run(0, 500) # fig=sim.figures[0] # fig.savefig('predprey_dynamics1.pdf') # fig.savefig('predprey_dynamics1.png') # fig=sim.figures[1]