Code to solve stochastic elliptic PDE with standard Monte Carlo