add_noise_disc_gaussian function

Add discrete normal noise with mean 0 to predicted values with constant variance