A method to replace a continuous univariate distribution with a discrete distribution that takes MN different values is analysed. Both distributions share the same rth moments for r = 0,..., 2N - 1 ...