#!/usr/bin/env python
# coding: utf-8

# # Tyson1991

# In[1]:


get_ipython().run_line_magic('matplotlib', 'inline')
from ecell4.prelude import *


# This model is described in the article:

# In[2]:


ref = citation(1831270)
ref


# In[3]:


from IPython.display import display, Markdown
display(Markdown("Abstract: {}".format(ref.abstract)))


# In[4]:


with reaction_rules():
    YT = Y + YP + M + pM
    CT = C2 + CP + M + pM

    ~Y > Y | 0.015 / CT
    Y > ~Y | 0.0 * Y
    CP + Y > pM | 200.0 * CP * Y / CT
    pM > M | pM * (0.018 + 180 * ((M / CT) ** 2))
    M > pM | 0.0 * M
    M > C2 + YP | 1.0 * M
    YP > ~YP | 0.6 * YP
    C2 > CP | 1000000.0 * C2
    CP > C2 | 1000.0 * CP

m = get_model()


# In[5]:


show(m)


# In[6]:


run_simulation(100.0, model=m, y0={'CP': 0.75, 'pM': 0.25}).plot(legend=True)