pwSimulateNewRealization

 mode = 1 (default):

 Simulates new data points based on the
 original data points plus Gaussian noise
 with a standard deviation as defined in the observables.

 Repeated call of this function always adds the noise
 to the original values (which are backuped internally).

 mode = 2 (Monte-Carlo):

 The new data is based on the model values plus noise:

 yData = yModel + randn(size(yModel)) .* yDataStd;

 mode = 3 (Bootstrap):

 All data points are drawn with replacement from the set of all available data points.