This paper introduces a procedure for solving a non-linear system of a partial differential equations of the parabolic type with random coefficients. An adaptive method transform the system into a set of ordinary differential equations which is solved after a suitable expansion of the probability density of the unknown random field in terms of orthogonal polynomials.