This study aimed to minimize the tumor cell population by using minimal medicine for chemotherapy treatment while maintaining the effector-immune cell population at a healthy threshold. Therefore, a mathematical model was developed in the form of ordinary differential equations (ODE), and the solution to Multi-Objective Optimal Control Problem (MOOCP) was obtained using Multi-Objective Optimization algorithms. In this study, the interaction of the tumor-effector cell population with chemotherapy was investigated using Pure MOOCP and Hybrid MOOCP methods. The handling of constraints and the Pontryagin Maximum Principle (PMP) differ among these methods. Swarm Intelligence (SI) and Evolutionary Algorithms (EA) were used to process the results of these methods. The numerical outcomes of SI and EA are displayed via Pareto Optimal Front. In addition, the solutions from these algorithms were further analyzed using the Hypervolume Indicator. The findings of this study demonstrate that the hybrid method outperforms Pure MOOCP via Multi-Objective Differential Evolution (MODE). The MODE produces a point on the Pareto Optimal Front with minimal distance to the origin, where the distance represents the best solution.