The modified Newton-Sabatier method is extended to the inverse scattering problem of coupled channels at fixed energy. The problem is solved for a system of coupled radial Schrödinger equations with a potential matrix independent of the angular momentum of the relative motion. The new method is applied to the example of two channels coupled by complex-valued square well potentials.