This work develops a local theory of the inhomogeneous coupled Schrödinger equations [Formula: see text]. Here, one treats the critical Sobolev regime [Formula: see text], where [Formula: see text] is the index of the invariant Sobolev norm under the dilatation [Formula: see text]. To the authors’ knowledge, the technique used in order to prove the existence of an energy local solution to the above-mentioned problem in the sub-critical regime s < s c, which consists of dividing the integrals on the unit ball of [Formula: see text] and its complementary, is no more applicable for s = s c. In order to overcome this difficulty, one uses two different methods. The first one consists of using Lorentz spaces with the fact that [Formula: see text], which allows us to handle the inhomogeneous term. In the second method, one uses some weighted Lebesgue spaces, which seem to be suitable to deal with the inhomogeneous term | x|− γ. In order to avoid a singularity of the source term, one considers the case p ≥ 2, which restricts the space dimensions to N ≤ 3.