A procedure for the realization of a transfer function matrix using a mixture of Markov parameters and moments is given. The procedure results in the determination of

, and

, and does not involve the determination of A
-1and

, etc., as needed in existing techniques.