Modelled COVID-19 pandemic with a system of 9 first order differential equations. The system was fitted to the values of the pandemic in Italy, UK, India, Brazil and Sweden, and numerically solved using MCMC statistical methods in python’s lmfit module. Estimates of the real number of infected people and predictions for the future were then made.