Table 2 Matlab function expressing the boolean network.

   new = zeros(lnow,1);

   Lig1 = 1;

Lig2 = 2;

A1 = 3;

   A2 = 4;

Inputk3r = 5;

K1 = 6;

   K2 = 7;

K3p1 = 8;

K3p2 = 9;

   K3r = 10;

Ph1 = 11;

Ph2 = 12;

   R1 = 13;

R2 = 14;

Geneexp = 15;

   Tf1 = 16;

   for i = 1:lnow

inode = setnow(i) ;

switch inode

case A1

new(i) = nodes(R1);

case A2

new(i) = and((nodes(R1) == 0), nodes(R2));

case K1

new(i) = and((nodes(Ph1) == 0), nodes(A1));

case K2

new(i) = or(and(nodes(A2), nodes(Ph2)), and(nodes(K1), nodes(A2)));

case K3p1

new(i) = and(nodes(K1), nodes(K3r));

case K3p2

new(i) = and(nodes(K3r), nodes(K2));

case K3r

new(i) = nodes(Inputk3r);

case Ph1

new(i) = nodes(K3p1);

case Ph2

new(i) = nodes(A1);

case R1

new(i) = nodes(Lig1);

case R2

new(i) = nodes(Lig2);

case Geneexp

new(i) = nodes(Tf1);

case Tf1

new(i) = or(or(and(nodes(K3p1), nodes(K3p2)), nodes(K3p2)), nodes(K3p1));

   end

end