Abstract In this work a multi-point constraint unfitted finite element method for the solution of Poisson equation is presented. Key features approach are strong enforcement essential boundary, and interface conditions. This, along with stability method, achieved through use constraints that applied to so-called ghost nodes lie outside physical domain. Another key benefit lies in fact that, as ...