diff --git a/src/pysagas/cfd/oblique_prandtl_meyer.py b/src/pysagas/cfd/oblique_prandtl_meyer.py
index fc3f5c5..51f2535 100644
--- a/src/pysagas/cfd/oblique_prandtl_meyer.py
+++ b/src/pysagas/cfd/oblique_prandtl_meyer.py
@@ -507,7 +507,7 @@ def oblique_rho2_rho1(M1: float, beta: float, gamma: float = 1.4):
         Peter Jacobs
         """
         M1n = M1 * abs(np.sin(beta))
-        rho2_rho1 = (gamma + 1) * M1n**2 / 2 + (gamma - 1) * M1n**2
+        rho2_rho1 = (gamma + 1) * M1n**2 / (2 + (gamma - 1) * M1n**2)
         return rho2_rho1
 
     @staticmethod