diff --git a/ElectromagneticScattering/scattering_init.py b/ElectromagneticScattering/scattering_init.py
index 75999adcf8b99b5ca63b631b3b0e192fad87ea7d..51761f3fdbebc097872ebb94cca385f5e6920df5 100644
--- a/ElectromagneticScattering/scattering_init.py
+++ b/ElectromagneticScattering/scattering_init.py
@@ -38,7 +38,7 @@ def field_VSH_expansion(post_filename):
     FF_erCrossXnm_p = np.zeros((npts_theta,npts_phi,p_max),dtype=complex)
     #####################################
     ##### sn, pn ,un
-    ##### Brian's recurrence relations
+    ##### Brian Stout's recurrence relations (B_ arrays)
     B_Pnpms        = np.zeros((npts_theta,m_max+1,n_max+1))
     B_unpms        = np.zeros((npts_theta,m_max+1,n_max+1))
     B_snpms        = np.zeros((npts_theta,m_max+1,n_max+1))
@@ -88,7 +88,7 @@ def field_VSH_expansion(post_filename):
     
     #####################################
     ##### sn, pn ,un
-    ##### Brian's recurrence relations
+    ##### Brian Stout's recurrence relations (B_ arrays)
     m_max = n_max
     p_max = n_max*n_max+2*n_max
     aM_nm   = np.zeros(p_max,dtype=complex)