An unconditionally stable numerical method for the Luo-Rudy 1 model used in simulations of defibrillation.

Numerical simulations of defibrillation using the Bidomain model coupled to a model of membrane kinetics represent a serious numerical challenge. This is because very high voltages close to defibrillation electrodes demand that extreme time step restrictions be placed on standard numerical schemes, e.g.