diff --git a/demos/demo_references.bib b/demos/demo_references.bib
index 31fe279..0277537 100644
--- a/demos/demo_references.bib
+++ b/demos/demo_references.bib
@@ -1,4 +1,4 @@
-@Article{LeVeque:1996,
+@article{LeVeque:1996,
   author={LeVeque, Randall J.},
   title={{High-Resolution Conservative Algorithms for Advection in Incompressible Flow}},
   journal={SIAM Journal on Numerical Analysis},
@@ -9,7 +9,7 @@ @Article{LeVeque:1996
   doi={10.1137/0733033}
 }
 
-@Book{Hundsdorfer:2003,
+@book{Hundsdorfer:2003,
   title={Numerical solution of time-dependent advection-diffusion-reaction equations},
   author={Hundsdorfer, Willem H and Verwer, Jan G and Hundsdorfer, WH},
   volume={33},
diff --git a/demos/gray_scott.py b/demos/gray_scott.py
index 8ef350d..c059405 100644
--- a/demos/gray_scott.py
+++ b/demos/gray_scott.py
@@ -160,5 +160,5 @@ def qoi():
 #
 # .. rubric:: References
 #
-# .. bibliography:: demo_references.bib
+# .. bibliography::
 #    :filter: docname in docnames
diff --git a/demos/point_discharge2d.py b/demos/point_discharge2d.py
index b1eee69..5db8e0d 100644
--- a/demos/point_discharge2d.py
+++ b/demos/point_discharge2d.py
@@ -9,7 +9,7 @@
 #
 # Consider the same steady-state advection-diffusion test case as in the
 # motivation for the Goalie manual: the "point discharge with diffusion"
-# test case from :cite:`RCJ:14`. In this test case, we solve
+# test case from :cite:`Riadh:2014`. In this test case, we solve
 #
 # .. math::
 #   \left\{\begin{array}{rl}
@@ -39,11 +39,11 @@ def get_function_spaces(mesh):
 
 
 # Point sources are difficult to represent in numerical models. Here we
-# follow :cite:`WBHP:22` in using a Gaussian approximation. Let
+# follow :cite:`Wallwork:2022` in using a Gaussian approximation. Let
 # :math:`(x_0,y_0)=(2,5)` denote the point source location and
 # :math:`r=0.05606388` be a radius parameter, which has been calibrated
 # so that the finite element approximation is as close as possible to the
-# analytical solution, in some sense (see :cite:`WBHP:22` for details). ::
+# analytical solution, in some sense (see :cite:`Wallwork:2022` for details). ::
 
 
 def source(mesh):
@@ -210,3 +210,9 @@ def qoi():
 # metric-based mesh adaptation to the point discharge test case considered here.
 #
 # This tutorial can be dowloaded as a `Python script <point_discharge2d.py>`__.
+#
+#
+# .. rubric:: References
+#
+# .. bibliography::
+#    :filter: docname in docnames
diff --git a/demos/solid_body_rotation.py b/demos/solid_body_rotation.py
index af807b5..1c52a74 100644
--- a/demos/solid_body_rotation.py
+++ b/demos/solid_body_rotation.py
@@ -234,5 +234,5 @@ def qoi():
 #
 # .. rubric:: References
 #
-# .. bibliography:: demo_references.bib
+# .. bibliography::
 #    :filter: docname in docnames
diff --git a/goalie/math.py b/goalie/math.py
index bd8971d..a361ff0 100644
--- a/goalie/math.py
+++ b/goalie/math.py
@@ -9,7 +9,7 @@ def bessi0(x):
     """
     Modified Bessel function of the first kind.
 
-    Code taken from :cite:`VVP+:92`.
+    Code taken from :cite:`Vetterling:1992`.
     """
     if isinstance(x, np.ndarray):
         if np.isclose(x, 0).any():