diff --git a/AerialNavigation/rocket_powered_landing/rocket_powered_landing.py b/AerialNavigation/rocket_powered_landing/rocket_powered_landing.py
index 239f3629c1..bff35b2c27 100644
--- a/AerialNavigation/rocket_powered_landing/rocket_powered_landing.py
+++ b/AerialNavigation/rocket_powered_landing/rocket_powered_landing.py
@@ -462,11 +462,11 @@ def __init__(self, m, K):
         # x_t+1 = A_*x_t+B_*U_t+C_*U_T+1*S_*sigma+zbar+nu
         constraints += [
             self.var['X'][:, k + 1] ==
-            cvxpy.reshape(self.par['A_bar'][:, k], (m.n_x, m.n_x)) @
+            cvxpy.reshape(self.par['A_bar'][:, k], (m.n_x, m.n_x), order='F') @
             self.var['X'][:, k] +
-            cvxpy.reshape(self.par['B_bar'][:, k], (m.n_x, m.n_u)) @
+            cvxpy.reshape(self.par['B_bar'][:, k], (m.n_x, m.n_u), order='F') @
             self.var['U'][:, k] +
-            cvxpy.reshape(self.par['C_bar'][:, k], (m.n_x, m.n_u)) @
+            cvxpy.reshape(self.par['C_bar'][:, k], (m.n_x, m.n_u), order='F') @
             self.var['U'][:, k + 1] +
             self.par['S_bar'][:, k] * self.var['sigma'] +
             self.par['z_bar'][:, k] +
@@ -536,20 +536,21 @@ def solve(self, **kwargs):
             with warnings.catch_warnings():  # For User warning from solver
                 warnings.simplefilter('ignore')
                 self.prob.solve(verbose=verbose_solver,
-                            solver=solver)
+                            solver=cvxpy.ECOS)
         except cvxpy.SolverError:
             error = True
 
         stats = self.prob.solver_stats
 
-        info = {
-            'setup_time': stats.setup_time,
-            'solver_time': stats.solve_time,
-            'iterations': stats.num_iters,
-            'solver_error': error
-        }
+        if stats is not None:
+            info = {
+                'setup_time': stats.setup_time,
+                'solver_time': stats.solve_time,
+                'iterations': stats.num_iters,
+                'solver_error': error
+            }
 
-        return info
+            return info
 
 
 def axis3d_equal(X, Y, Z, ax):
diff --git a/requirements/requirements.txt b/requirements/requirements.txt
index cd4a01f88f..c3b0fa8511 100644
--- a/requirements/requirements.txt
+++ b/requirements/requirements.txt
@@ -1,7 +1,7 @@
 numpy == 2.2.0
 scipy == 1.14.1
 matplotlib == 3.10.0
-cvxpy == 1.5.3
+cvxpy == 1.6.0
 pytest == 8.3.4 # For unit test
 pytest-xdist == 3.6.1 # For unit test
 mypy == 1.13.0 # For unit test