We present our preliminary results from applying the Newton-Krylov-Schwarz method for the simulation of the electrical activity of the heart in two dimensions. We use the bidomain nonlinear equations, using a fully implicit time discretization scheme, and solving the resulting large system of equations with a Newton based algorithm at each step. We incorporate anisotropy into our model, and compare the results of using various conductivity ratios between the intra-cellular and extra-cellular domains. We also compare our results with previous work that has been done using less precise techniques, such as explicit and semi-implicit schemes, and less detailed models, such as the Monodomain model, and isotropic as well as quasi-isotropic Bidomain models. The results are obtained using PETSc of the Argonne National Laboratory.