Isoparametric FEM in Matlab

A short Matlab implementation realizes a flexible isoparametric finite element method up to quadratic order for the approximation of elliptic problems in two-dimensional domains with curved boundaries. Triangles and quadrilaterals equipped with varying quadrature rules allow for mesh refinement. Numerical examples for the Laplace equation with mixed boundary conditions indicate the flexibility of isoparametric finite elements.