This function computes ideal data and perturbs them with random gaussian noise with standard deviation X_std and Y_std. All these perturbations are done using X and Y data saved in 'original' arrays so the function may be run multiple times and will use original source datasets for simulations. To keep X or Y from being perturbed, - simply enter 0 in place of X_std or Y_std.
For details see source code in: Dataset.m