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

function [Dx,Dy,Dz,ND]=Dcube(a,b,c,n)

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

Dx=zeros(n^3,1);
Dy=zeros(n^3,1);
Dz=zeros(n^3,1);

ND=0;
for i=1:n
    for j=1:n
        for k=1:n
                ND=ND+1;
                Dx(ND,1)=x(i);
                Dy(ND,1)=y(j);
                Dz(ND,1)=z(k);
        end
    end
end