Solution of modified matrix equations

An efficient computational algorithm is developed for solving linear matrix equations whose characteristic interaction matrix ${\bf H}'$ differs by a matrix of rank r from an interaction matrix ${\bf H}$ of a system whose solution is known, or easily computed. The formulation clearly manifests the effects of physical parameters characterizing the interactions in the system. The exact inverse is given for a nonsymmetric Toeplitz matrix with arbitrary corner elements.