The parameterized Cauer ladder network (CLN) method, which is a recently developed nonlinear model order reduction (MOR) method, was applied to the analysis of a variable reluctance (VR) stepper motor. To represent the magnetic saturation, we parameterized the stator network using the coil magnetic flux and input sequence. The rotor network adopted the CLN to investigate the effect of the eddy current. The CLN method was effective in analyzing the continuous rotation, stepwise rotation, and step-out motion.