function [c4n,n4e,Db,Nb] = triang_cube(d)
if d == 1
    c4n = [0;1]; n4e = [1,2]; Db = 1; Nb = 2;
elseif d == 2
    c4n = [0,0;1,0;1,1;0,1]; n4e = [1,2,3;1,3,4];
    Db = [1,2]; Nb = [2,3;3,4;4,1];
elseif d == 3  
    c4n = [0,0,0;1,0,0;1,1,0;0,1,0;0,0,1;1,0,1;1,1,1;0,1,1];
    n4e = [1,2,3,7;1,6,2,7;1,5,6,7;1,8,5,7;1,4,8,7;1,3,4,7];
    Db = [1,2,3;1,4,3];
    Nb = [2,3,7;2,7,6;1,2,6;1,6,5;5,6,7;1,8,5;5,7,8;1,4,8;...
        4,7,8;3,4,7];
end
