We investigate a set of nonlinear matrix equations with nonnegative matrix coefficients which has arisen in applied sciences. There are papers where the minimal nonnegative solution of the set of nonlinear matrix equations is computed applying the different procedures. The alternate linear implicit method and its modifications have intensively investigated because they have simple computational scheme. We construct a new decoupled modification of the alternate linear implicit procedure to compute the minimal nonnegative solution of the considered set of equations. The convergence properties of the proposed iteration are derived and a sufficient condition for convergence is derived. The performance of the proposed algorithm is illustrated on several numerical examples. On the basis of the experiments we derive conclusions for applicability of the computational schemes.