montecarlo sampling of the Neufeld solution