diff --git a/Six_1g_spherical_benchmarks/Problem_1.py b/Six_1g_spherical_benchmarks/Problem_1.py
index 1f7bcf1..f0fe599 100644
--- a/Six_1g_spherical_benchmarks/Problem_1.py
+++ b/Six_1g_spherical_benchmarks/Problem_1.py
@@ -77,7 +77,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=sigt, c=c)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
diff --git a/Six_1g_spherical_benchmarks/Problem_2.py b/Six_1g_spherical_benchmarks/Problem_2.py
index 1f7bcf1..f0fe599 100644
--- a/Six_1g_spherical_benchmarks/Problem_2.py
+++ b/Six_1g_spherical_benchmarks/Problem_2.py
@@ -77,7 +77,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=sigt, c=c)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
diff --git a/Six_1g_spherical_benchmarks/Problem_4.py b/Six_1g_spherical_benchmarks/Problem_4.py
index 286f82a..c84eb19 100644
--- a/Six_1g_spherical_benchmarks/Problem_4.py
+++ b/Six_1g_spherical_benchmarks/Problem_4.py
@@ -40,7 +40,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=1.0, c=0.0)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
diff --git a/Six_1g_spherical_benchmarks/Problem_5.py b/Six_1g_spherical_benchmarks/Problem_5.py
index 520b3b8..12a2f5e 100644
--- a/Six_1g_spherical_benchmarks/Problem_5.py
+++ b/Six_1g_spherical_benchmarks/Problem_5.py
@@ -40,7 +40,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=1.0, c=0.0)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)