Abstract: In this paper we consider the singular perturbation boundary-value problem of the following coupling type system of convection-diffusion equations We advance two methods: the first one is the initial value solving method, by which the original boundary-value problem is changed into a series of unperturbed initial-value problems of the first order ordinary differential equation or system so that an asymptotic expansion is obtained;the second one is the boundary-value solving method, by which the original problem is changed into a few boundary-value problems having no phenomenon of boundary-layer so that the exact solution can be obtained and any classical numerical methods can be used to obtain the numerical solution ofconsismethods can be used to obtain the numerical solution of consistant high accuracy with respect to the perturbation parameter ε.