This paper presents a novel numerical technique for the identification of effective and basic reproduction numbers,ReandR0, for long-term epidemics, using an inverse problem approach. The method is based on the direct integration of the SIR (Susceptible–Infectious–Removed) system of ordinary differential equations and the least-squares method. Simulations were conducted using official COVID-19 data for the United States and Canada, and for the states of Georgia, Texas, and Louisiana, for a period of two years and ten months. The results demonstrate the applicability of the method in simulating the dynamics of the epidemic and reveal an interesting relationship between the number of currently infectious individuals and the effective reproduction number, which is a useful tool for predicting the epidemic dynamics. For all conducted experiments, the results show that the local maximum (and minimum) values of the time-dependent effective reproduction number occur approximately three weeks before the local maximum (and minimum) values of the number of currently infectious individuals. This work provides a novel and efficient approach for the identification of time-dependent epidemics parameters.