A numerical code, which utilizes the stellarator expansion, is developed and tested. It is used to investigate the magnetohydrodynamic stability properties of several stellarator configurations, including Heliotron E, Wendelstein VII‐A, a modular‐coil device, and ATF‐1.