Current electrophysiological models vary greatly in both complexity and accuracy. These models often involve a set of partial differential equations coupled with a large number of stiff ordinary differential equations (ODEs). The system of ODEs attempt to simulate the flow of ionic currents present in the cellular level of the heart and has been continually developed to provide an increasingly detailed description of cellular physiology. However, the stiffness has the effect of decreasing the speed of the solving process since the stability is often the limiting factor for the solution. An efficient method is needed to significantly decrease computational time in the numerical solution of two and three dimensional models of the electrical activity of the myocardium. In this presentation, we will be exploring the usage and efficiency of Gauss-type nested implicit Runge-Kutta technique to solve cardiac cell models. The method is of order 4 and has only explicit internal stage that leads to practical implementations. Comparison with other numerical methods employed in the context of electrocardiology will be presented.