added new Monte Carlo procedures with variance reduction methods. Note that the best variance reduction for the run time seems to be with the antithetic variable technique, not moment matching or other more complex techniques.