ellipticModel with general b.c.
We need a test for mixed type b.c. - specially partly defined: for example square domain with Gamma_0 = Robin Gamma_1 = Dirichlet Gamma_2 = Neuman Gamma_3 = natural possibly non linear. We need an approach to define the different parts of the boundary for Robin/Neuman