package("polmat")

math

 

P:=matrix([[85/128, 5/16],[5/16, 85/128]]) + matrix([[17/64, 1/8],[1/8, 17/64]])*(z+1/z)

math

 

polmat::toeplitz(P, z, 2)

math