A modification of Boris Algorithm in cylindrical coordinate used to solve Lorentz equation for axisymmetrical plasma simulation is proposed. The velocity updates are calculated directly using matrix operation. The inertial forces are evaluated at n-0.5 time steps. The algorithm is tested for several scenarios: electric field only, electric field dominated, |E| = |B|, ExB drift, gyration test, and a test devised by Delzanno. The result is compared to the standard Boris solver in cylindrical coordinate. The results show that even though generally the solver is not that accurate compared to the standard Boris solver, especially for a larger time step width, it is comparatively faster. Thus, the solver might be useful for plasma simulations that do not require detailed trajectory for each simulated particle.