This manuscript is related to undertaking a mathematical model (susceptible, vaccinated, infected, and recovered) of rotavirus. Some qualitative results are established for the mentioned challenging childhood disease epidemic model of rotavirus as it spreads across a population with a heterogeneous rate. The proposed model is investigated using a novel approach of fractal calculus. We compute the boundedness positivity of the solution of the proposed model. Additionally, the basic reproduction ratio and its sensitivity analysis are also performed. The global stability of the endemic equilibrium point is also confirmed graphically using some available values of initial conditions and parameters. Sufficient conditions are deduced for the existence theory, the Ulam–Hyers (UH) stability. Specifically, the numerical approximate solution of the rotavirus model is investigated using efficient numerical methods. Graphical presentations are presented corresponding to a different fractional order to understand the transmission dynamics of the mentioned disease. Furthermore, researchers have examined the impact of lowering the risk of infection on populations that are susceptible and have received vaccinations, producing some intriguing results. We also present a numerical illustration taking the stochastic derivative of the proposed model graphically. Researchers may find this research helpful as it offers insightful information about using numerical techniques to model infectious diseases.