A number of new layer methods solving the Neumann problem for semilinear parabolic equations is constructed by using probabilistic representations of their solutions. The methods exploit the ideas of weak sense numerical integration of stochastic differential equations in bounded domain. In spite of the probabilistic nature these methods are nevertheless deterministic. Some convergence theorems are proved. Numerical tests are presented
