How to increase the speed of this code?

Here, I want to generate a 3D lattice points in the form of 3 column tuples. Here is my code:

    % code
  pts = [];
for i = 0.01:0.01:1
for j = 0.01:0.01:1
for k = 0.01:0.01:1
a = [i j k];
pts = [pts;a];
end
end
end
  end

It is too slow when I run this script in the Matlab, is there any way to increase the speed of the operation? Or is there any way to decrease the number of the for loop?

Thanks

Answer 1

[Roger Stafford]

Try 'ndgrid':

   [X,Y,Z] = ndgrid(0.01:0.01:1);
   pts = [Z(:),Y(:),X(:)];   %(Corrected)

Answer 2

[Azzi Abdelmalek]

i = 0.01:0.01:0.1
j = 0.01:0.01:.1
k = 0.01:0.01:.1
[ii,jj,kk]=meshgrid(i,j,k)
out=[kk(:) ii(:) jj(:)]

[Roger Stafford]

@Azzi: I think these will not be in the same order requested by Zhongruo.