Subroutine for calculation of matrix padé approximants - CPC, Mendeley (Mendeley, 1976)