How to create a symbolic matrix where each entry is a function of time?

I am trying to create a matrix where each entry of the matrix is a function of time. For example, this is what I want the output to look like:

A =
[ A1_1(t), A1_2(t)]
[ A2_1(t), A2_2(t)]

I have tried:

syms t
A = sym('A(t)', [2 2])

but this does not work. Any suggestions?

Answer 1

[Walter Roberson]

   A = sym(zeros(2,2));
   for J = 1 : size(A,1)
     for K = 1 : size(A,2)
       A(J,K) = feval(symengine,sprintf('A%d_%d', J, K), 't');
     end
   end