Tehran University of Medical Sciences Journal of Biostatistics and Epidemiology 2383-4196 7 1 2021 05 15 25 35 Lakshmi Sridhar Department of Chemical Engineering , University of Puerto Rico, Mayaguez, Puerto Rico. 2021 01 05 2021 02 02 Introduction: A rigorous multiobjective optimal control strategy (that does not require the use of weighting functions) of the epidemic models that consider vaccination and treatment strategies is presented. Modifications of the standard susceptible-infectious-removed, susceptible-exposed-infectious-removed, and the modified susceptible-infectious-removed models are dynamically optimized to minimize the number of infected individuals while, controlling the rate at which the individuals are vaccinated and treated. Method: The optimization program, Pyomo , where the differential equations are automatically converted to a Nonlinear Program using the orthogonal collocation method is used for performing the dynamic optimization calculations. The Lagrange-Radau quadrature with three collocation points and 10 finite elements are chosen. The resulting nonlinear optimization problem was solved using the solver BARON 19.3, accessed through the Pyomo-GAMS27.2 interface.  Results: The computational results how that the multiobjective optimal control profiles generated by this strategy are very similar to those produced when weighting functions are used.  Conclusion: The main conclusion of this work is to demonstrate that one can perform a rigorous dynamic optimization of epidemic models without the use of weighting functions that have the potential to produce some uncertainty and doubt in the results, in addition to dealing with unnecessary additional variables. https://jbe.tums.ac.ir/index.php/jbe/article/view/423