Solve coupled PDE

Hello. How can I input the following problem in FREEFEM?


It is actually a coupled PDE, and the coefficients in the weak form are matrices instead of scalars.

Thanks in advance

You have no problem to use freefem++ to solve numerical this problem if you have amesh of \mathcal D not to fine.

I takes 1/3 a hours to write a first solution.