Discussion Forum

I have made a physics builder that solves an eigenvalue differential equation given some input geometry and input potential V(x,y,z). I would like to be able to set V to take different values in different parts of the domain. How can I set this? My script is attached. The geometry is two hexagons with common origin and a slab on the outer hexagon. I would like the function V which enters in my differential equation to take on three different values, one for each of the three geometries. How can I do that? Could anyone be so kind as to show how to do it in the attached script.
Edit: Right now I am trying with a geometry of two concentric circles with different radii. I am then using an expression like a + (b - a) * (x^2 + y^2 > R1^2) - b * (x^2 + y^2 > R2^2) for the stepwise potential. Is there any way to do this for a more complicated geometry like my hexagonal one?