What works is the following:
I assume that the function lam is a heaviside function, representing lam1 in region 1 and lam2 in region 2.Then I have to consider this in the PDE and it works. But I am sure this is not the most elegant approach. Any ideas to solve this problem in an easier way?
Attachments:
