The level set method is a powerful way of tracking surfaces by defining the surface as a zero level set of a continuous function that is usually a signed distance function. The level set method is one of the best methods for simulating multi-phase flow because it can easily handle fast topological changes, as well as splitting and merging of fluids. In this thesis, a standard level set method was implemented in C++, using the finite element method library deal.II, to simulate incompressible two-phase flow on some benchmark problems. The results show a significant change of mass in the simulations, something that should not be allowed to happen when simulating incompressible fluids. The mass changes mainly occur in the reinitialization phase, where the level set function is rebuilt to look more like a signed distance function.