Abstract
We are concerned with discussing the numerical approximation for the random heat model with boundary conditions of the third type in one dimension. The random (forward time, backward space) finite difference scheme is used to find the stochastic solution. The impression of the consistency and the stability technique under the mean square spectral radius is studied.