We consider in this work hexagonal finite volume methods for two-dimensional applications. Ordinary and compact seven-point schemes with three-color iterative updates are investigated with application to elliptic equations and reaction-diffusion systems. In particular, electrical activity of human ventricular tissue is observed. Monodomain model for human ventricular cell is used to simulate action membrane potential, and electrical wave propogation acrossing two dimensional myocardiac tissue can be obtained. The results have shown good agreement with previous research works.