% Cube surface data set [0,a]*[0,b]*[0,c]

function [Rx,Ry,Rz,NM]=Bcube(a,b,c,n)

[x,y,z]=cube_eq(a,b,c,n);

NM=0;
for i=1:n
    for j=1:n
        for k=1:n
            if(i==1|i==n|j==1|j==n|k==1|k==n)
                NM=NM+1;
                Rx(NM)=x(i);
                Ry(NM)=y(j);
                Rz(NM)=z(k);
            end
        end
    end
end