function a=r3d(array,nrow,i,j,k);
%function a=r3d(array,nrow,i,j,k);

[ncol,mcol]=size(array);
ncol=ncol/nrow;
if i==0
     i=1:nrow;
end
if j==0
     j=1:ncol;
end
if k==0
     k=1:mcol;
end
ind=i'*ones(j/j(1))+nrow*ones(i'/i(1))*(j-1);
a=array(ind(:),k);
if (length(k)==1)&(length(j)>1)
    a=remat(a(:),length(i));
end

