We investigate mathematical models of the Hepatitis B and C viruses in the study, considering vaccination effects into account. By utilising fractional and ordinary differential equations, we prove the existence of equilibrium and the well-posedness of the solution. We prove worldwide stability with respect to the fundamental reproduction number. Our numerical techniques highlight the biological relevance and highlight the effect of fractional derivatives on temporal behaviour. We illustrate the relationships among susceptible, immunised, and infected populations in our epidemiological model. Using comprehensive numerical simulations, we analyse the effects of fractional derivatives and highlight solution behaviours. Subsequent investigations will examine the impact of regional heterogeneity, providing significant perspectives for epidemiological research.
