ON THE NUMERICAL SOLUTION OF SYSTEM OF NONLINEAR PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS