The purpose of this study is to perform the numerical results of the global stability of infectious disease mathematical model by using the stochastic computing scheme. The design of proposed solver is presented by one of the efficient and reliable schemes named as Bayesian regularization neural network (BRNN). The global stability of infectious disease mathematical nonlinear model is categorized into susceptible, infected, recovered and vaccinated. The construction of dataset is performed through the Runge-Kutta scheme in order to lessen the mean square error (MSE) by dividing the statics as training 74 %, while 13 % for both testing and endorsement. The proposed stochastic process contains log-sigmoid merit function, twenty neurons and optimization through RBNN for the numerical solutions of the global stability of infectious disease mathematical system. The best training values for each model's case are performed around 10–11. The scheme's correctness is performed by the matching of the results and the minor calculated absolute error performances. Moreover, the regression, state transmission, error histogram and MSE indicate the trustworthiness of the designed solver.