#---------------------------------------------------------
# Code for Macaulay2

Q = ZZ/2[M,B,A,L,P] / ideal(M^2-M,B^2-B,A^2-A,L^2-L,P^2-P)

RingElement | RingElement :=(x,y)->x+y+x*y;
RingElement & RingElement :=(x,y)->x*y;

fM = A;
fB = M;
fA = A | (L & B);
fL = P | (L & (1+B));
fP = M;

(fM, fB, fA, fL, fP)

I = ideal(fM+M, fB+B, fA+A, fL+L, fM+M);
G = gens gb I

#---------------------------------------------------------
# Code for Cyclone in AlgoRun

NUMBER OF VARIABLES: 5
NUMBER OF STATES :   2

M = A
B = M
A = B*A*L+B*L+A
L = B*L*P+B*L+L*P+L+P
P = M

#---------------------------------------------------------
# Code for BoolNet in RStudio

targets, factors
M, A
B, M
A, A | (L & B)
L, P | (L & !B)
P, M


library(Boolnet)
lacModel <- loadNetwork("lac-operon")
lacAttractors <- getAttractors(lacModel)
plotStateGraph(lacAttractors)