annotate private/randnat.m @ 15:a6d5597bd922

Trying alternative fixed point computation in mc_global_info1>mc_fixpt2;
Also factored out randnat from some other functions.

function X=randnat(M)
% randnat - Draw random natural number from uniform distribution
%
% randnat :: 
%    M:natural   ~'maximum value', 
% -> [[D]->[M]] ~'size D array of numbers in 1..M'.
%
% The size is specified as in the RAND function, except that if
% only one dimension is given, we create a vector, not a square
% matrix, eg

X=1+floor(M*rand);