Merge pull request #1579 from othmar-humboldt/patch-1

Update change_ecc_inc.py to allow for incl. reduction

src/poliastro/core/thrust/change_ecc_inc.py

@@ -85,7 +85,7 @@ def a_d(t0, u_, k_):
             np.cos(nu)
         )  # The sign of ß reverses at minor axis crossings
 
-        w_ = cross(r_, v_) / norm(cross(r_, v_))
+        w_ = (cross(r_, v_) / norm(cross(r_, v_))) * np.sign(inc_f - inc_0)
         accel_v = f * (np.cos(beta_) * thrust_unit + np.sin(beta_) * w_)
         return accel_v