Random walk-on-spheres algorithms for solving mixed and Neumann boundary-value problems

We propose a new approach to constructing Monte Carlo methods for solving mixed boundary value problems for elliptic equations with constant coefficients. We derived a mean-value relation for point values of the solution. As a consequence, the walk-on-spheres algorithm can still be used even after a trajectory hits the reflecting boundary. Such an approach is significantly more efficient than the standard one.