In this study, a nonlinear deterministic mathematical model that evaluates two important therapeutic measures of the COVID-19 pandemic: vaccination of susceptible and treatment for infected people who are in quarantine, is formulated and rigorously analyzed. Some of the fundamental properties of the model system including existence and uniqueness, positivity, and invariant region of solutions are proved under a certain meaningful set. The model exhibits two equilibrium points: disease-free and endemic equilibrium points under certain conditions. The basic reproduction number,
R
0
, is derived via the next-generation matrix approach, and the dynamical behavior of the model is explored in detail. The analytical analysis reveals that the disease-free equilibrium solution is locally as well as globally asymptotically stable when the associated basic reproduction number is less than unity which indicates that COVID-19 dies out in the population. Also, the endemic equilibrium point is globally asymptotically stable whenever the associated basic reproduction number exceeds a unity which implies that COVID-19 establishes itself in the population. The sensitivity analysis of the basic reproduction number is computed to identify the most dominant parameters for the spreading out as well as control of infection and should be targeted by intervention strategies. Furthermore, we extended the considered model to optimal control problem system by introducing two time-dependent variables that represent the educational campaign to susceptibles and continuous treatment for quarantined individuals. Finally, some numerical results are illustrated to supplement the analytical results of the model using MATLAB ode45.