1. I'll post here a simplified version of my problem Say you have a matrix A, and you want all its components to be functions, for example: A11 = Sin(a_n) A12 = Cos(a_n) A21 = Sin(b_n) A22 = Cos(b_n) And I want to be able to do this in mathematica so as to have the matrix A a function of n. 2. Relevant equations How can I do this, such that when I ask MatrixForm[A] for example, it would give me the matrix with all components in terms of a_4 and b_4 ? 3. The attempt at a solution Rather long and wrong, so better not post it here.