The virtual element method (VEM) is a relatively new discretization technique that can be seen as a generalization of the finite element method. In VEM there is no need to know explicitly the finite element space in each element—thus the virtual name. The method allows for arbitrary order approximations on polygonal meshes, and the polygons need to even be convex! In this project we looked at the method from an engineering perspective, comparing it with standard FEM. We also looked at creating hybrid VEM and FEM meshes. We showed the method converges optimally and can provide tremendous flexibility in analyzing complex problems. In the figure below, the method is applied to the Flamant problem, a half-space loaded with a point force.