Surface Integral Equations are often used to model electromagnetic scattering phenomena. Large-scale problems can efficiently be solved with Boundary Element Methods (BEM) of which the Method of Moments (MoM) has found widespread use in the computational electromagnetics community. The framework of integral equations allows for the design of many different formulations for a wide range of scattering problems. Most of them are a clever combination of the basic electric and magnetic field integral operators. In this paper, the open-source library BEM++ will be used as a powerful tool to build different integral equation formulations and preconditioners.