Different countries used the growth Gompertz function at the beginning of the COVID-19 pandemic to model the number of cumulative infected cases since it provides reasonable results. Such a model allows only one mode, but the pandemic evolution has exhibited a multimodal behavior due to the different waves and variants of the COVID-19 virus. Thus, Gompertz’s classical growth model is not well suited to describe a long pandemic with different virus variants. This work presents generalizations of the Gompertz model that can reproduce a multimodal behavior to model the dynamics of infected cases. The models are applied to COVID-19 data from Nuevo Leon, Mexico.Diferentes países usaron la función de crecimiento Gompertz al principio de la pandemia por COVID-19 para modelar el numero acumulado de infectados dado que proporcionaba un ajuste razonable. Este modelo permite una única moda, pero la pandemia evoluciono exhibiendo un comportamiento multimodal debido a las diferentes olas y variantes del COVID-19. Por tanto, el modelo Gompertz clásico de crecimiento no ajusta bien para describir una pandemia larga con diferentes variantes del virus. Este trabajo presenta generalizaciones del modelo Gompertz donde se pueda capturar un comportamiento multimodal para modelar la dinámica de los casos infectados. Este modelo es aplicado a datos de COVID-19 de Nuevo León, México.