Hi, I'm glad to announce that the resample() function also works if the optional 'blocksize' parameter is set to 1. As naturally expected, this is then equivalent to omitting the parameter, giving a standard one-by-one resampling. This behavior is not documented, however; in the docs it says "If this argument is given it should be a positive integer greater than or equal to 2." thanks, sven -- Example: <hansl> matrix m = mnormal(10,2) print m matrix mr = resample(m, 1) print mr </hansl>