diff --git a/RELEASES.md b/RELEASES.md
index 2b2edcd9d..745a7de67 100644
--- a/RELEASES.md
+++ b/RELEASES.md
@@ -7,6 +7,7 @@
 - Added feature `grad=last_step` for `ot.solvers.solve` (PR #693)
 - Automatic PR labeling and release file update check (PR #704)
 - Reorganize sub-module `ot/lp/__init__.py` into separate files (PR #714)
+- Fix documentation in the module `ot.gaussian` (PR #718)
 
 #### Closed issues
 - Fixed `ot.mapping` solvers which depended on deprecated `cvxpy` `ECOS` solver (PR #692, Issue #668)
diff --git a/ot/gaussian.py b/ot/gaussian.py
index 4645d5fa4..002e69fb4 100644
--- a/ot/gaussian.py
+++ b/ot/gaussian.py
@@ -354,7 +354,7 @@ def bures_wasserstein_barycenter(
 
     The function estimates the optimal barycenter of the
     empirical distributions. This is equivalent to resolving the fixed point
-     algorithm for multiple Gaussian distributions :math:`\left{\mathcal{N}(\mu,\Sigma)\right}_{i=1}^n`
+    algorithm for multiple Gaussian distributions :math:`\left\{\mathcal{N}(\mu,\Sigma)\right\}_{i=1}^n`
     :ref:`[1] <references-OT-mapping-linear-barycenter>`.
 
     The barycenter still following a Gaussian distribution :math:`\mathcal{N}(\mu_b,\Sigma_b)`
@@ -452,7 +452,7 @@ def empirical_bures_wasserstein_barycenter(
 
     The function estimates the optimal barycenter of the
     empirical distributions. This is equivalent to resolving the fixed point
-     algorithm for multiple Gaussian distributions :math:`\left{\mathcal{N}(\mu,\Sigma)\right}_{i=1}^n`
+    algorithm for multiple Gaussian distributions :math:`\left\{\mathcal{N}(\mu,\Sigma)\right\}_{i=1}^n`
     :ref:`[1] <references-OT-mapping-linear-barycenter>`.
 
     The barycenter still following a Gaussian distribution :math:`\mathcal{N}(\mu_b,\Sigma_b)`