???
Watch this Author's files
Was just playing around: L34: If you make it an mpower operation this function can support vector inputs.
fib = round( alpha.^n / sqrt(5) );
Jos's randpermfull better vectorized.
http://www.mathworks.com/matlabcentral/fileexchange/23257-randpermfull
Contact us at files@mathworks.com