compute FreeFem formal expression of the rotational
v = rot(u)
a column vector, ffeq data structure (tlist)
a ffeq data structure (tlist)
u = defvar(['u1';'u2']) rot(u)