This function uses original X and ideal Y data to produce randomly perturbed datasets for Monte-Carlo analysis.
For details see source code in: Dataset.m