diff --git a/notebooks/AU-Mic.ipynb b/notebooks/AU-Mic.ipynb
index 754f376..106d73e 100644
--- a/notebooks/AU-Mic.ipynb
+++ b/notebooks/AU-Mic.ipynb
@@ -114,13 +114,28 @@
     "T_phot = 3700  # Plavchan (2009)\n",
     "T_spot1 = T_spot2 = T_spot3 = 2300\n",
     "\n",
-    "blackbody = lambda *args: ActiveStar()._blackbody(*args)\n",
+    "# blackbody = lambda *args: ActiveStar()._blackbody(*args)\n",
     "\n",
-    "phot = jnp.array([float(blackbody(tess_mean_wavelength, T_phot))])\n",
+    "# Download and bin PHOENIX model spectra to compute contrast:\n",
+    "kwargs = dict(\n",
+    "    bins=1, # Turn the TESS bandpass into a single bin\n",
+    "    min = 0.7 * u.micron,\n",
+    "    max = 0.9 * u.micron,\n",
+    "    log=False\n",
+    ")\n",
+    "\n",
+    "photosphere, spot1, spot2, spot3 = [\n",
+    "    bin_spectrum(\n",
+    "        get_spectrum(T_eff=T_eff, log_g=4.52, Z=0.12, cache=True), **kwargs\n",
+    "    )\n",
+    "    for T_eff in [T_phot, T_spot1, T_spot2, T_spot3]\n",
+    "]\n",
+    "\n",
+    "phot = jnp.array([float(photosphere.flux[0].value)])\n",
     "spectrum = jnp.array(\n",
-    "    [[blackbody(tess_mean_wavelength, T_spot1)],\n",
-    "     [blackbody(tess_mean_wavelength, T_spot2)],\n",
-    "     [blackbody(tess_mean_wavelength, T_spot3)]]\n",
+    "    [[spot1.flux[0].value],\n",
+    "     [spot2.flux[0].value],\n",
+    "     [spot3.flux[0].value]]\n",
     ")"
    ]
   },
@@ -605,7 +620,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.5"
+   "version": "3.10.13"
   }
  },
  "nbformat": 4,