diff --git a/MITRotor/BEMSolver.py b/MITRotor/BEMSolver.py
index daea9bc..559e657 100644
--- a/MITRotor/BEMSolver.py
+++ b/MITRotor/BEMSolver.py
@@ -130,7 +130,7 @@ def Ctprime(self, grid: Literal["sector", "annulus", "rotor"] = "rotor"):
         return average(self.geom, Ctprime, grid=grid)
 
 
-@adaptivefixedpointiteration(max_iter=500, relaxations=[0.25, 0.5, 0.96])
+@adaptivefixedpointiteration(max_iter=1000, relaxations=[0.25, 0.5, 0.96])
 class BEM:
     """
     A generic BEM class which facilitates dependency injection for various models.
diff --git a/MITRotor/TangentialInduction.py b/MITRotor/TangentialInduction.py
index cd493ec..d44369c 100644
--- a/MITRotor/TangentialInduction.py
+++ b/MITRotor/TangentialInduction.py
@@ -57,7 +57,7 @@ def __call__(
         eff_yaw = calc_eff_yaw(yaw, tilt)
         aprime = (
             np.clip(aero_props.C_tau_corr, -2, 2)
-            / (4 * np.maximum(geom.mu_mesh, 0.1) ** 2 * tsr * (1 - aero_props.an) * np.cos(eff_yaw))
+            / (4 * np.maximum(geom.mu_mesh, 0.1) * tsr * (1 - aero_props.an) * np.cos(eff_yaw))
         )
 
         return aprime