B:=[ [ [ 1, 0 ], [ 1, 0 ], [ 1, 0 ], [ 1, 0 ] ],
     [ [ 2, 2, 1 ], [ 2, 3, 1 ], [ 2, 4, 2 ], [ 2, 4, 3 ] ],
     [ [ 4, 4, 2, 1, 3 ] ], [ ] ];;
Y:=RegularCWComplex(B);;
OrientRegularCWComplex(Y);
Print(Y!.orientation,"\n");
C:=ChainComplexOfRegularCWComplex(Y);
Print(C,"\n");
Print(List([0..3],n->C!.dimension(n)),"\n");
Print(BoundaryMatrix(C,1)*BoundaryMatrix(C,2),"\n");
