This function computes Y data using current best-fit parameters to use as a source dataset for Monte-Carlo perturbation.
For details see source code in: Dataset.m