El documento está dedicado al análisis teórico e implementación numérica de un modelo matemático de un sistema de reacción-difusión no lineal en la plataforma COMSOL Multiphysics. Se considera el problema aplicado de la simulación por computadora del cambio de polarización en películas ferroeléctricas delgadas. El modelo se basa en el enfoque termodinámico de Landau-Ginzburg-Devonshire-Khalatnikov y se formaliza como un problema de valor inicial-límite para una ecuación diferencial parabólica semilineal. Se explicaron los fundamentos teóricos del modelo. La aplicación del diseño de la interfaz de usuario se desarrolló con COMSOL Multiphysics. Se realizaron una serie de experimentos computacionales para estudiar la histéresis ferroeléctrica y las dependencias de la temperatura de la polarización en el ejemplo de una película de titanato de bario ferroeléctrico.