He arreglado un pequeño error que no permitía la simulación para posiciones negativas del soliton. También he introducido la función "atan2" para calcular la fase de la función de onda, aunque inesperadamente sigo teniendo el mismo problema de saltos en "2pi"
He arreglado un pequeño error que no permitía la simulación para posiciones negativas del soliton. También he introducido la función "atan2" para calcular la fase de la función de onda, aunque inesperadamente sigo teniendo el mismo problema de saltos en "2pi"