This paper is concerned with the dynamics of a system of nonlinear finite difference equations which arise from a class of parabolic boundary-value problems. The main purpose is the determination of a critical physical parameter σ* such that for σ * the solution blows up at some finite n*. Also discussed is the convergence and blowing up property of un, for various initial vector uo whenσ< σ*.