Limpié(simplifiqué) un poco el código
En la función actualizar metí los osciladores
Modifiqué la función tiempo_colision_particula_oscilador, aumentando n para aumentar la precisión en el cálculo del tiempo