An advanced integration method for solving reaction-diffusion-type equation
s for cardiac conduction is suggested. Operator splitting and adaptive time
step methods were used in this method, which ran significantly speed up in
tegration while preserving accuracy.