Abstract
In this manuscript, an HIV/AIDS epidemic model including five compartments of S,I,A,T and R is proposed under fractal-fractional-order derivative. The existence theory utilizing Schaefer- and Banach-type fixed point theorems for the solution of considered model is constructed. Additionally, Ulam-Hyers (UH) and generalized UH stability conditions via nonlinear functional analysis are established. A fractional type of two-step Lagrange polynomial known as fractional Adams-Bashforth (AB) method is developed for numerical simulation of the considered model. The simulated results for various fractal-fractional orders are tested on some existing real data of disease spread in South Africa and show that the values of S,I,A,T are decreased after treatment was started. In addition, R, the population of infected people who have changed their sexual habits sufficiently after starting the treatment and changing their sexual behavior, increased gradually. Finally, it is shown that for all five compartments of the proposed model of HIV/AIDS, the smaller values of fractal-fractional order have better performance than larger values.