I came across an old post and notice an syntax different from what I am costume in WM. The original post is here
(*definition of the modified spherical Bessel function of the first \
kind from https://mathematica.stackexchange.com/questions/10767/fast-spherical-\
harmonics-radiative-transfer*)
SphericalBesselI[0, 0] := SphericalBesselI[0, 0] = 1;
SphericalBesselI[l_, z_] := Sqrt[\[Pi]/(2 z)] BesselI[l + 1/2, z];
(*expansion coefficients for the Henyey-Greenstein phase function*)
\[Sigma][l_][g_] := g^l
(*diagonal matrix elements*)
h[l_][g_, \[Omega]_] :=
h[l][g, \[Omega]] = (2 l + 1) (1 - \[Omega] \[Sigma][l][g]);
(*low-order determinants*)
Look at line 7 and 9. It has a double brackets. Not sure but will [Sigma][l_][g_] := g^l be same as [Sigma][l_,g_] := g^l?