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

# In[1]:


import actonet
import biolqm


# In[2]:


model = biolqm.load("https://raw.githubusercontent.com/jgtz/StableMotifs/master/TLGLNetwork.txt", "booleannet")


# In[3]:


bn = biolqm.to_minibn(model, ensure_boolean=True)
inputs = bn.constants()
inputs


# In[4]:


for n in inputs:
    bn[n] = n


# In[5]:


act = actonet.ActoNet(bn, inputs)


# In[6]:


r = act.reprogramming_fixpoints(Apoptosis=1, maxsize=2)
r


# In[7]:


r.as_table()


# In[8]:


r = act.reprogramming_fixpoints(Apoptosis=0, maxsize=2)
r


# In[9]:


r.as_table()


# In[ ]: