In [2]:

    

from dzdy import *
import matplotlib.pyplot as plt

da = Director()
da.load_pc('scripts/pBAD.txt')
da.load_dc('scripts/BAD.txt')


    

In [3]:

    

bad = da.new_mc('BAD', 'ABM', tar_pc='pBAD', tar_dc='BAD')
## Key command
bad.add_behaviour('BR', 'LifeRate', s_death='Dead', s_birth='Young', rate=0.2, dt=1)
bad.set_observations(states=['Young', 'Middle', 'Old', 'Dead'], transitions=['Die'])


    

In [4]:

    

mod, out =  da.simulate('BAD', y0={'Young': 100, 'Middle': 100, 'Old': 100}, to=20)


    

In [5]:

    

out.plot()
plt.show()