This research presents a method of modeling and numerical simulation of a reluctance stepper motor using reduced finite-element time-stepping technique. In presented model, the circuit equations are reduced to non-stationary differential equations, i.e. the inductance mapping technique is used to find relationship between coil inductance and rotor position. A strongly coupled field-circuit model of the stepper motor is presented. In analyzed model the magnetostatic field partial differential equations are coupled with rotor motion equation and solved simultaneously in each iterative step. The nonlinearity problem is solved using Newton-Raphson method with spline approximation of the B-H curve.