In this paper, a five-compartment model is used to explore the dynamics of the COVID-19 pandemic, taking the vaccination campaign into account. The present model consists of five components that lead to a system of five ordinary differential equations. In this paper, we examined the disease from the perspective of a fractal fractional derivative in the Caputo sense with a power law type kernal. The model is also fitted with real data for Pakistan between June 1, 2020, and March 8, 2021. The fundamental mathematical characteristics of the model have been investigated thoroughly. We have calculated the equilibrium points and the reproduction number for the model and obtained the feasible region for the system. The existence and stability criteria of the model have been validated using the Banach fixed point theory and the Picard successive approximation technique. Furthermore, we have conducted stability analysis for both the disease-free and endemic equilibrium states. On the basis of sensitivity analysis and the dynamics of the threshold parameter, we have estimated the effectiveness of vaccination and identified potential control strategies for the disease using the proposed model outbreaks. The stability of the concerned solution in Ulam-Hyers and Ulam-Hyers-Rassias sense is also investigated. For the proposed problem, some results regarding basic reproduction numbers and stability analysis for various parameters are represented graphically. Matlab software is used for numerical illustrations. Graphical representations are given for different fractional orders and for various parametric values.

In this paper, we investigate a fractal-fractional-order mathematical model with the influence of hospitalized patients and the impact of vaccination with fractal-fractional operators. The respective derivatives are considered in the Caputo, Caputo Fabrizio, and Atangana–Baleanu senses of fractional order $$ \alpha $$α and fractal dimension $$ \tau $$τ . For the proposed problem, some results regarding basic reproduction number and stability are given. Using the next-generation matrix approach, we have investigated the global and local stability of several types of equilibrium points. We provide a detailed analysis of the existence and uniqueness of the solution. Moreover, we fit the model with the real data of Pakistan from June 01, 2020, till March 24, 2021. Then, we use the fractal-fractional derivative to find a numerical solution for the model. MATLAB software is used for numerical illustration. Graphical presentations corresponding to different parameteric values are given as well.