An efficient computational algorithm is developed for solving linear matrix equations whose characteristic interaction matrix H' differs by a matrix of rank r from an interaction matrix H of a system ...