The idea is to approximate the solution of a variational problem using a finite-dimensional space of piecewise polynomials that interpolate homogeneous boundary conditions at boundary nodes. To the curved region,we use the interpolation method, uses the normal interpolation in the unit interior, but the unit on uses the Lobatto multinomial interpolation, then gives its error estimate. |