Simulations of cardiac defibrillation are associated with considerable numerical challenges. The cell models have traditionally been discretized by first order explicit schemes, which are associated with severe stability issues. The sharp transition layers in the solution call for stable and efficient solvers. We propose a second order accurate numerical method for the Luo-Rudy phase 1 model of electrical activity in a cardiac cell, which provides sequential update of each governing ODE. An \textit{a priori} estimate for the scheme is given, showing that the bounds of the variables typically observed during electric shocks constitute an invariant region for the system, regardless of the time step chosen. Thus the choice of time step is left as a matter of accuracy. Conclusively, we demonstrate the theoretical result by some numerical examples, illustrating second order convergence for the Luo-Rudy 1 model.