3 coupled elliptic PDE's.How to implement? Finite difference derivative?

I made a beginning