vertices = 4*[...
    0 0 0
    0 0 1
    0 1 0
    0 1 1
    1 0 0
    1 0 1
    1 1 0
    1 1 1];
    
stickers = [...
    0 1 1
    0 1 3
    0 3 1
    0 3 3
    1 0 1
    1 0 3
    3 0 1
    3 0 3
    1 1 0
    1 3 0
    3 1 0
    3 3 0
    4 1 1
    4 1 3
    4 3 1
    4 3 3
    1 4 1
    1 4 3
    3 4 1
    3 4 3
    1 1 4
    1 3 4
    3 1 4
    3 3 4];
faces = [...
    1 2 6 5 1
    1 3 4 2 1
    8 4 3 7 8
    6 8 7 5 1
    2 4 8 6 2
    1 5 7 3 1];

SquareTemp1 = [1 2; 3 4]; SquareTemp2 = [1 4; 2 3];
curFactStick1 = [2 3]; curFaceStic2 = [1 4];
x=vertices(:,1); y=vertices(:,2); z=vertices(:,3);
stickNum =1;

hold on;
% for i=1:8
% text(x(i),y(i),z(i),num2str(i));
% end

for faceIt=1:6
    curFace = faces(faceIt,:);
    curSquare1 = curFace(SquareTemp1);
    curSquare2 = curFace(SquareTemp2);
    
    pl1 = plot3(x(faces(faceIt,:)), y(faces(faceIt,:)), ...
        z(faces(faceIt,:)));
    midx1 = (x(curSquare1(:,1)) + x(curSquare1(:,2)))/2;
    midy1 = (y(curSquare1(:,1)) + y(curSquare1(:,2)))/2;
    midz1 = (z(curSquare1(:,1)) + z(curSquare1(:,2)))/2;
    
    midx2 = (x(curSquare2(:,1)) + x(curSquare2(:,2)))/2;
    midy2 = (y(curSquare2(:,1)) + y(curSquare2(:,2)))/2;
    midz2 = (z(curSquare2(:,1)) + z(curSquare2(:,2)))/2;
    
    plot3(midx1,midy1,midz1);
    plot3(midx2,midy2,midz2);

    for i=1:size(stickers,1)
        text(stickers(i,1), stickers(i,2), stickers(i,3), num2str(i));
    end
end

permMat = zeros(6, 24);
for i=1:6
    permMat(i,:) = 1:24;
end

permMat(1, 1:4) = [3 1 4 2];
permMat(2, 5:8) = [6 8 5 7];
permMat(3, 9:12) = [11 9 12 10];
permMat(4, 13:16) = [15 13 16 15];
permMat(5, 17:20) = [18 20 17 19];
permMat(6, 21:24) = [23 21 24 23];

diffPerms(1,:) = [22 21 6 5 9 10 17 18];
diffPerms(2,:) = [21 23 14 13 11 9 1 2];
diffPerms(3,:) = [3 1 5 7 13 15 19 17];
diffPerms(4,:) = [24 23 8 7 11 12 19 20];
diffPerms(5,:) = [22 24 16 15 12 10 3 4];
diffPerms(6,:) = [14 16 20 18 4 2 6 8];

for it = 1:6
    imagePerm = [diffPerms(it,3:8) diffPerms(it,1:2)];
    permMat(it, diffPerms(it,:)) = imagePerm;
end

myperms = {permMat(1,:), permMat(2,:), permMat(3,:), permMat(4,:), ...
    permMat(5,:), permMat(6,:)};