diff --git a/aerosandbox/optimization/opti.py b/aerosandbox/optimization/opti.py
index 355ccf5f..63688ca2 100644
--- a/aerosandbox/optimization/opti.py
+++ b/aerosandbox/optimization/opti.py
@@ -1168,6 +1168,7 @@ def constrain_derivative(self,
         )
         duals = self.subject_to(
             np.diff(variable) == integrals,
+            _stacklevel=_stacklevel + 1
         )
 
         return duals