This paper deals with the construction of stable discrete numerical solutions of strongly coupled diffusion mixed partial differential systems. After discretization, using a matrix difference scheme, the resulting coupled mixed partial difference problem is treated using a discrete separation of variables method which avoids solving algebraic systems. Existence, stability of solutions, its construction and illustrative examples are considered.