Code covered by the BSD License
function partition(ii) global mark mark(ii)=mark(ii)-1; for jj=1:3 mark(ii+jj*2^(2*mark(ii)))=mark(ii); end
Contact us at files@mathworks.com