Implementation of Romberg's variance reduction method for Monte Carlo, in R