A new crossover dynamics mathematical model of monkeypox disease based on fractional differential equations and the <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" altimg="si553.svg" display="inline" id="d1e2699"><mml:mi>Ψ</mml:mi></mml:math

A novel crossover model for monkeypox disease that incorporates Ψ-Caputo fractional derivatives is presented here, where we use a simple nonstandard kernel function Ψ(t). We can be obtained the Caputo and Caputo–Katugampola derivatives as special cases from the proposed derivative. The crossover dynamics model defines four alternative models: fractal fractional order, fractional order, variable order, and fractional stochastic derivatives driven by fractional Brownian motion over four time intervals. The Ψ-nonstandard finite difference method is designed to solve fractal fractional order, fractional order, and variable order mathematical models. Also, the nonstandard modified Euler Maruyama method is used to study the fractional stochastic model. A comparison between Ψ-nonstandard finite difference method and Ψ-standard finite difference method is presented. Moreover, numerous numerical tests and comparisons with real data were conducted to validate the methods’ efficacy and support the theoretical conclusions.