THREE-DIMENSIONAL LOW-THRUST INTERPLANETARY TRAJECTORY OPTIMIZATION