Central Difference Convexification Method for Soft-Landing Trajectory Optimization in Planetary Powered Descending Phase