diff --git a/examples/gaussian_nondiagcov_legacy.py b/examples/gaussian_nondiagcov_legacy.py
index f36d6bb8..a411937f 100644
--- a/examples/gaussian_nondiagcov_legacy.py
+++ b/examples/gaussian_nondiagcov_legacy.py
@@ -279,7 +279,7 @@ def run_example(
             # 2D surface plot of posterior.
             # ==================================================================
             ln_posterior_func = partial(ln_posterior, inv_cov=inv_cov)
-            ln_posterior_grid, x_grid, y_grid = ex_utils.eval_func_on_grid(
+            ln_posterior_grid, x_grid, y_grid = hm.utils.eval_func_on_grid(
                 ln_posterior_func,
                 xmin=xmin,
                 xmax=xmax,