Matlab - Productoria

Buen día. Agradezco su ayuda:

Tengo una matriz que representa el tiempo ( de tamaño m x n) ordenados asi:

t1 t4 t7

t2 t5 t8

t3 t6 t9 ...
.
.
Para cada uno de esos tiempos t, se va a evaluar una matriz cuadrada y se quiere hacer el producto de ellas asi: donde f1, f2 f3 son cierta función f, evaluada cada elemento de la 1er comunda de t

M1=[ f1(t1), f2(t1) ; f3(t1), f4(t1) ] ,
M2=[ f1(t2), f2(t2) ; f3(t2), f4(t2) ]
M2=[ f1(t3), f2(t3) ; f3(t3), f4(t3) ]
Mn=[ f1(tn), f2(tn) ; f3(tn), f4(tn) ]

Hacer el producto de M1.*M2.*M3 .* Mn y almacenar esa matriz resultado

luego similar pero para la otra comumna de t4 t5 y t6
: donde f1, f2 f3 son cierta función f, evaluada cada elemento de la 2er comunda de t


M1=[ f1(t4), f2(t4) ; f3(t4), f4(t4) ] ,
M2=[ f1(t5), f2(t5) ; f3(t5), f4(t5) ]
M2=[ f1(t6), f2(t6) ; f3(t6), f4(t6) ]
Mn=[ f1(tn), f2(tn) ; f3(tn), f4(tn) ]


De nuevo hacer el producto de M1.*M2.*M3.*Mn y almacenar esa matriz resultado y así hasta terminar con todas las comunas de la matriz de tiempo t.