Numerical Solutions for Optimal Control of Monodomain Equations in Cardiac Electrophysiology