corrected blackbody luminosity density equation

.mailmap

@@ -48,6 +48,8 @@ Barnabás Barna <[email protected]>
 
 Caroline Sofiatti <[email protected]>
 
+Cecelia Powers <[email protected]>
+
 Chaitanya Kolliboyina <[email protected]>
 
 Chinmay Talegaonkar <[email protected]>

docs/physics/montecarlo/initialization.ipynb

@@ -216,7 +216,7 @@
    "metadata": {},
    "source": [
     "We define important constants, and for comparison's sake, we code the Planck distribution function\n",
-    "$$L_\\nu (\\nu)=\\frac{8\\pi r_\\mathrm{boundary\\_inner}^2 h\\nu^3}{c^2}\\frac{1}{\\exp\\left(\\frac{h\\nu}{k_BT_\\mathrm{inner}}\\right)-1}$$\n",
+    "$$L_\\nu (\\nu)=\\frac{8\\pi^2 r_\\mathrm{boundary\\_inner}^2 h\\nu^3}{c^2}\\frac{1}{\\exp\\left(\\frac{h\\nu}{k_BT_\\mathrm{inner}}\\right)-1}$$\n",
     "where $L_\\nu$ is the luminosity density (see [Basic Spectrum Generation](../spectrum/basic.ipynb)) with respect to frequency, $\\nu$ is frequency, $h$ is Planck's constant, $c$ is the speed of light, and $k_B$ is Boltzmann's constant:\n"
    ]
   },
@@ -237,7 +237,7 @@
     "def planck_function(nu):\n",
     "    return (\n",
     "        8\n",
-    "        * np.pi\n",
+    "        * np.pi**2\n",
     "        * r_boundary_inner**2\n",
     "        * h\n",
     "        * nu**3\n",