En el contexto del Modelo de Dos Dobletes de Higgs tipo III (2HDM-III), que incluye acoplamientos con cambio de sabor leptónico estudiamos la producción de los bosones de Higgs neutros en colisiones electrón-positrón y su decaimiento con cambio de sabor a muón antitau, σ(e+e− → Zφ → Zµτ¯), donde φ : (h0, H0), en el futuro colisionador llamado International Linear Collider (ILC). Se analizaron varios escenarios (Benchmark Points) encontrándose que el número de eventos para este proceso es muy grande teniendo la posibilidad de ser detectado en ILC. Los cálculos asociados a este proceso se realizaron por medio de los softwares LanHEP y CalcHEP, los cuales se eligieron porque ofrecen un manejo intuitivo de los lagrangianos, las reglas de Feynman y manejo de los parámetros para el cálculo de las secciones eficaces, los anchos de decaimiento y las razones de ramificación del proceso de interés, con lo que se calcula el número de eventos y con esto comprobar la predictibilidad del modelo.