A numerical scheme to study the propagation of action potentials in cardiac tissue using the bidomain model is presented. Previous methods have used explicit or implicit time stepping schemes. Explicit methods are easily implemented, but have severe time step restrictions. With implicit methods, there are no time step restrictions; however, a nonlinear system of equations must be solved at each time step. The scheme discussed here is a semi-implicit time step scheme which eliminates the severe time step restriction and requires that only linear systems be solved at each time step. These linear systems are solved using a multigrid method. |