n := 3; V := combinat[powerset](n); read(`lowerbound-data.maple`); read(`lowerbound-programs.maple`): lList := [seq( seq(seq([converttonumber(switch(o,v)),i],v = V), o=OO2[i]), i=1..nops(OO2))]: List := map2(op,1,lList): # List contains now a list of ALL orientations Classnumber := map2(op,2,lList): Sinks := seq(sinks(o)[1],o=List): dualprob := # probabilities for the orientations in each orientation class [39691/16433616, 3845/4108404, 3845/4108404, 2977/16433616, 3427/2738936, 1945/2054202, 19129/16433616, 19129/8216808, 1609/1027101, 0, 39691/8216808, 1029/1369468, 0, 3845/4108404, 0, 3427/1369468, 0,0,1810/342367]; add(dualprob[Classnumber[i]],i=1..nops(List)) = 1; erg := bestforSearcher({$1..nops(List)}); evalf(erg);