Higher-order basis functions are widely used to model currents and fields in numerical simulations of electromagnetics problems because of the greater accuracy and computational efficiency they can provide. Different problem formulations, such as method of moments (MoM) and the finite element method (FEM) require different constraints on basis functions for optimal performance, such as normal or tangential continuity between cells. In this thesis, a method of automatically generating bases that satisfy the desired basis constraints is applied to a MoM formulation for scattering problems using surface integral equations. Numerical results demonstrate the accuracy of this approach, and show good system matrix conditioning when compared to other higher-order bases.