Currently most of the numerical resolution of optimal control problems is done using direct methods where the increase accuracy of indirect methods is overshadowed by the necessary analytical derivation required beforehand. With recent developments from the control-toolbox ecosystem team the application of indirect methods as become more streamline enabling a wider range of problems to be solved, like, for example, optimal control problems applied to the transmission of infectious diseases. In this work, we aim to extend the application of indirect methods to optimal control problems applied to epidemiological models, using the control-toolbox [25].