This work considers the derivation procedure and evaluation of dynamic equations for isotropic micropolar circular cylinders by adopting a power series expansion method in the radial coordinate. Variationally consistent equations of motion together with pertinent sets of boundary conditions are expressed in a systematic fashion up to arbitrary order. The numerical results cover eigenfrequencies, mode shapes and field distributions over cross sections for axisymmetric and flexural motion adopting different sets of end boundary conditions for equations of different truncation orders of the present method. The results illustrate that the present approach may render benchmark solutions provided that higher order equations are used, and act as accurate approximate engineering solution for lower order equations.