Because the Coronavirus epidemic spread in Iraq, the COVID-19 epidemic of people quarantined due to infection is our application in this work. The numerical simulation methods used in this research are more suitable than other analytical and numerical methods because they solve random systems. Since the Covid-19 epidemic system has random variables coefficients, these methods are used. Suitable numerical simulation methods have been applied to solve the COVID-19 epidemic model in Iraq. The analytical results of the Variation iteration method (VIM) are executed to compare the results. One numerical method which is the Finite difference method (FD) has been used to solve the Coronavirus model and for comparison purposes. The numerical simulation methods which are Mean Monte Carlo Finite difference (MMC_FD) and Mean Latin Hypercube Finite difference (MLH_FD), are also used to solve the proposed epidemic model under study. The obtained results are discussed, tabulated, and represented graphically. Finally, the absolute error is the tool used to compare the numerical simulation solutions from 2020 to 2024 years. The behavior of the Coronavirus in Iraq has been expected for 4 years from 2020 to 2024 using the proposed numerical simulation methods.
