A simple method for obtaining a stable reduced-order (

th-order) approximant

for a stable hlgh-order (

th-order) transfer function

is described. The method is an improvement over some previous methods [1], [2] in that it not only ensures matching of the first

time moments of

and

, but it also ensures that the (

)th time moment of

be as close as possible to the corresponding value in

.