Este estudio introduce el método de función auxiliar óptima (OAFM) para resolver un sistema no lineal de ecuaciones de Belousov-Zhabotinsky. El sistema se caracteriza por su dinámica compleja y se trata utilizando el operador de Caputo y conceptos del cálculo fraccional. El OAFM proporciona un enfoque sistemático para obtener soluciones analíticas aproximadas mediante la construcción de una función auxiliar. Al optimizar los parámetros de la función auxiliar, se obtiene una solución aproximada que coincide estrechamente con el comportamiento del sistema original. La efectividad y precisión del OAFM se demuestran a través de simulaciones numéricas y comparaciones con métodos existentes. El cálculo fraccional mejora la comprensión y modelado de la dinámica no lineal en el sistema de Belousov-Zhabotinsky. Este estudio contribuye al cálculo fraccional y a la dinámica no lineal, ofreciendo una herramienta poderosa para analizar y resolver sistemas complejos como la ecuación de Belousov-Zhabotinsky.