Abstract <p>A mathematical model is developed to analyze the dynamics of corruption by identifying significant parameters. The model allows us to simulate and observe changes in outcomes when these parameters are varied using Python. The model is proved to be both epidemiologically and mathematically well posed. We showed that all solutions of the model are positive and bounded with initial conditions in a certain meaningful set. The existence of unique corruption free and endemic equilibrium points are investigated and the basic reproduction number is computed. Then, we study the local asymptotic stability of these equilibrium points. The analysis shows that the system has a locally asymptotically stable corruption-free equilibrium point when the reproduction number is less than one and locally asymptotically stable endemic equilibrium point when the reproduction number is greater than one. The simulation result shows the agreement with the analytical results.</p>