Multi-rendezvous low-thrust trajectory optimization using costate transforming and homotopic approach