Las ondas gravitacionales producidas por fusiones de estrellas de neutrones binarias ofrecen una ventana única para entender el comportamiento de la materia en condiciones extremas. En este contexto, modelamos analíticamente el efecto de la materia en las ondas gravitacionales de las fusiones de estrellas de neutrones binarias. Comenzamos con un sistema de agujeros negros binarios, aprovechando el formalismo post-newtoniano para la fase de inspiral y el modelo de un cuerpo hacia atrás para la fusión. Combinamos los dos métodos para generar una forma de onda base y validamos nuestros resultados contra simulaciones de relatividad numérica. A continuación, integramos los efectos de marea en fase y amplitud para tener en cuenta la interacción entre la materia y el espacio-tiempo utilizando el modelo NRTidal y probamos su precisión contra las predicciones de relatividad numérica para dos ecuaciones de estado, encontrando una discrepancia alrededor de la fusión. Posteriormente, levantamos la restricción sobre los coeficientes para que sean independientes de la deformabilidad de marea y los recalibramos utilizando las predicciones de relatividad numérica. Obtenemos mejores ajustes para la fase y la amplitud alrededor de la fusión y somos capaces de extender el modelado de fase más allá de la fusión. Implementamos nuestro método en un nuevo código Python de código abierto y fácil de usar, guiado por un Jupyter Notebook, llamado . Nuestra investigación ofrece nuevas perspectivas sobre el modelado analítico del efecto de las mareas en las ondas gravitacionales de las fusiones de estrellas de neutrones binarias.