How do you multiply two transfer functions numerically?

I have two transfer function estimates for two blocks in my feedforward path that are in series. I would like to numerically multiply both of these together to form a single block equivalent. To estimate the transfer function for each block, I have experimental data and I have used the tfestimate function.

The values for each transfer function estimate looks like this. Of course, I've only put the first several values for each estimate here for brevity.

Txy1=
  -0.234430873110630 + 0.00000000000000i
  -0.235894779161570 - 0.00319325918105480i
  -0.240404763278006 - 0.00644705375023638i
  -0.248330648404303 - 0.00989933822199415i
  -0.260341562132277 - 0.0138713441364568i
  -0.277494712873904 - 0.0190444192662199i
Txy2=
5.64813402210633 + 0.00000000000000i
5.73106608167536 - 0.276024783164775i
5.98705210834285 - 0.574385252351786i
6.43786839588527 - 0.921841502572328i
7.11988058673111 - 1.35448893694931i
8.08153862775710 - 1.92307018423326i

I want to compute what Txy1*Txy2 is. Does anyone know how I should do this? I'm thinking that I would "foil" element by element?