<<~/Programs/OctagonalPet/OctaPet/Mathematica/bound2.m

(*This runs the main routine from bound2.m*)
  A=Table[Table[Table[Table[{i1,i2,i3,i4},{i1,0,1}],{i2,0,1}],{i3,0,1}],{i4,0,1}];
  A=Flatten[A,3];

  B=Table[Table[Table[Table[{i1,i2,i3,i4},{i1,1,4}],{i2,1,4}],{i3,1,4}],{i4,1,4}];
  B=Flatten[B,3];

 LIST1=Table[Table[GetBounds[A[[i]],B[[j]]],{i,1,Length[A]},{j,1,Length[B]}]]
 LIST2=Flatten[LIST1]
 Max[LIST2]
