Analytic regularity for a singularly perturbed system of reaction-diffusion equations with multiple scales

We consider a coupled system of two singularly perturbed reaction-diffusion equations, with two small parameters 0 < ε ≤ μ ≤ 1, each multiplying the highest derivative in the equations. The presence of these parameters causes the solution(s) to have boundary layers which overlap and interact, based on the relative size of ε and μ. We show how one can construct full asymptotic expansions together with error bounds that cover the complete range 0 < ε ≤ μ ≤ 1. For the present case of analytic input data, we present derivative growth estimates for the terms of the asymptotic expansion that are explicit in the perturbation parameters and the expansion order.