how do you code a probabilistic boolean gene network?

We are using the bioinformatics toolbox. It includes this information

wire = [4 2 3;5 3 2;3 6 1;5 4 6;6 1 2;3 5 6];

rule = {[231 230], 64, [5 7], 108, 61, [62 60]};

prob = {[1/2 1/2], 1, [1/2 1/2], 1, 1, [1/2 1/2]};

and it needs function [pnet, rnet] = ruletree(rule,prob)

and the ruletree needs this suggested output.

pnet = 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250 0.1250

rnet = 111111 111112 112111 112112 211111 211112 212111 212112

I just need a starting point, I am really lost here.