diff --git a/examples/bulk_j_an.py b/examples/bulk_j_an.py
index 285854498106a10cc67c6b9d77d5b5e5a1f16a11..8c46fac29c7517e4a5e82876a8da70d72f5ad2b7 100644
--- a/examples/bulk_j_an.py
+++ b/examples/bulk_j_an.py
@@ -54,7 +54,7 @@ def j_anom(T, h, n=5, alpha_soc=0.01, perturbative=False):
 
     # Solver for evaluating current perturbatively in SOC
     eta = 1.0
-    sol = get_solver(soc_alpha=alpha_soc, eta=eta, h=h, D=0)
+    sol = get_solver(soc_alpha=alpha_soc, eta=eta, h=h, D=1.0)
 
     E_typical = 10 + 10 * T
     w, a = get_matsubara_sum(T, E_typical)