Before I tried to build a boolean network, I think it is quite complicated. But it comes that building a boolean network is pretty simple. I even did it in MATLAB. IT does not need ode function to solve the problem. Just put the rate of change logic equations in a ordinary equation, to calculate several steps (step size always 1). You can get the answer. The sample code is listed in the following.
function test;
a = 1;
b = 1;
c = 1;
for i = 1:10
[a, b, c] = calabc(a, b, c);
[a b c]
end
function [a,b,c] = calabc(a0, b0, c0)
a = b0;
b = a0 | c0;
c = (a0 & b0) | (b0 & c0) | (a0 & c0);