Mathematical models are now an important tool for studying cardiac electrophysiology and arrhythmias. Utilizing these models to quantify behavior and make predictions requires solving the models computationally using numerical schemes. We discuss different solution methods and other computational considerations for simulating cardiac action potentials in single cells and multicellular preparations.