An algorithm is proposed for the accurate and efficient integration of the linear mixed hardening elastoplastic constitutive equations in finite element analysis based on the assumption of constant strain rate vector with respect to the plastic Lagrange multiplier over the time step. An exact sample problem and iso-error maps are given to test the accuracy of the present method as well as two other commonly used ones.