A new numerical method based on the method of auxiliary sources (MAS) is presented for analyzing the propagation characteristics of the EH1 mode of a half-mode substrate integrated waveguide (HMSIW). A set of filament auxiliary sources are located on a cylinder surface inside the shorting-via of the