diff --git a/desc/compute/_fast_ion.py b/desc/compute/_fast_ion.py
index 43f388c88..f86c30bd3 100644
--- a/desc/compute/_fast_ion.py
+++ b/desc/compute/_fast_ion.py
@@ -143,7 +143,7 @@ def _Gamma_c(params, transforms, profiles, data, **kwargs):
     )
 
     def Gamma_c(data):
-        """∫ dλ ∑ⱼ [v τ γ_c²]ⱼ π²/4."""
+        """∫ dλ 〈 ∑ⱼ [v τ γ_c²]ⱼ 〉 π²/4."""
         bounce = Bounce2D(
             grid,
             data,
@@ -297,7 +297,7 @@ def _Gamma_c_Velasco(params, transforms, profiles, data, **kwargs):
     )
 
     def Gamma_c(data):
-        """∫ dλ ∑ⱼ [v τ γ_c²]ⱼ π²/4."""
+        """∫ dλ 〈 ∑ⱼ [v τ γ_c²]ⱼ 〉 π²/4."""
         bounce = Bounce2D(
             grid,
             data,
diff --git a/desc/compute/_neoclassical.py b/desc/compute/_neoclassical.py
index bbc3fe629..6eb70d547 100644
--- a/desc/compute/_neoclassical.py
+++ b/desc/compute/_neoclassical.py
@@ -199,7 +199,7 @@ def _epsilon_32(params, transforms, profiles, data, **kwargs):
     )
 
     def eps_32(data):
-        """(∂ψ/∂ρ)⁻² B₀⁻³ ∫ dλ λ⁻² ∑ⱼ Hⱼ²/Iⱼ."""
+        """(∂ψ/∂ρ)⁻² B₀⁻³ ∫ dλ λ⁻² 〈 ∑ⱼ Hⱼ²/Iⱼ 〉."""
         # B₀ has units of λ⁻¹.
         # Nemov's ∑ⱼ Hⱼ²/Iⱼ = (∂ψ/∂ρ)² (λB₀)³ ``(H**2 / I).sum(axis=-1)``.
         # (λB₀)³ d(λB₀)⁻¹ = B₀² λ³ d(λ⁻¹) = -B₀² λ dλ.