simulating Gaussian domains

mainly some helper functions (creating distributions etc..), the
simulation itself is not complete yet.