diff --git a/SimPEG/EM/Static/DC/ProblemDC.py b/SimPEG/EM/Static/DC/ProblemDC.py
index 40d4183f..04a07c57 100644
--- a/SimPEG/EM/Static/DC/ProblemDC.py
+++ b/SimPEG/EM/Static/DC/ProblemDC.py
@@ -128,6 +128,10 @@ class Problem3D_N(BaseDCProblem):
         MeSigma = self.MeSigma
         Grad = self.mesh.nodalGrad
         A = Grad.T * MeSigma * Grad
+
+        # Handling ... singularity
+        A[0,0] = A[0,0] + 1.
+
         # if self._makeASymmetric is True:
         #     return V.T * A
         return A
diff --git a/tests/em/static/test_DC.py b/tests/em/static/test_DC.py
index 1d83d32c..99b384aa 100644
--- a/tests/em/static/test_DC.py
+++ b/tests/em/static/test_DC.py
@@ -104,7 +104,7 @@ class DCProblemTestsN(unittest.TestCase):
 
     def test_misfit(self):
         derChk = lambda m: [self.survey.dpred(m), lambda mx: self.p.Jvec(self.m0, mx)]
-        passed = Tests.checkDerivative(derChk, self.m0, plotIt=False, num=3)
+        passed = Tests.checkDerivative(derChk, self.m0, plotIt=False)
         self.assertTrue(passed)
 
     def test_adjoint(self):
@@ -120,7 +120,7 @@ class DCProblemTestsN(unittest.TestCase):
 
     def test_dataObj(self):
         derChk = lambda m: [self.dmis.eval(m), self.dmis.evalDeriv(m)]
-        passed = Tests.checkDerivative(derChk, self.m0, plotIt=False, num=3)
+        passed = Tests.checkDerivative(derChk, self.m0, plotIt=False)
         self.assertTrue(passed)
 if __name__ == '__main__':
     unittest.main()