This work presents a reliable algorithm to obtain approximate analytical solutions for a strongly coupled system of singularly perturbed convection–diffusion problems, which exhibit a boundary layer at one end. The proposed method involves constructing a zero-order asymptotic approximate solution for the original system. This approximation results in the formation of two systems: a boundary layer system with a known analytical solution and a reduced terminal value system, which is solved analytically using an improved residual power series approach. This approach combines the residual power series method with Padé approximation and Laplace transformation, resulting in an approximate analytical solution with higher accuracy compared to the conventional residual power series method. In addition, error estimates are extracted, and illustrative examples are provided to demonstrate the accuracy and effectiveness of the method.