compute FreeFem formal expression of the trace
v = tr(u)
a square matrix, ffeq data structure (tlist)
a ffeq data structure (tlist)
u = defvar(['u1' 'u3';'u2' 'u4']) tr(u)