lat = 41.02054 
lon = -72.15608
date = Date.new(2024, 11, 10)

observer = Astronoby::Observer.new(
  latitude: Astronoby::Angle.from_degrees(lat),
  longitude: Astronoby::Angle.from_degrees(lon)
)

moon = Astronoby::Moon.new(time: date.to_time)
moon_events = moon.observation_events(observer: observer) # <-- infinite loop hit

def iterate(initial_rising, initial_transit, initial_setting)
        delta = 1
        corrected_rising = initial_rising
        corrected_transit = initial_transit
        corrected_setting = initial_setting
        until delta < ITERATION_PRECISION <--- delta oscillates between two values > `ITERATION_PRECISION`
          iterate = RiseTransitSetIteration.new(
            observer: @observer,
            date: @date,
            coordinates_of_the_next_day: @coordinates_of_the_next_day,
            coordinates_of_the_day: @coordinates_of_the_day,
            coordinates_of_the_previous_day: @coordinates_of_the_previous_day,
            shift: shift,
            initial_rising: corrected_rising,
            initial_transit: corrected_transit,
            initial_setting: corrected_setting
          ).iterate
          delta = iterate.sum
          puts delta
          corrected_rising = rationalize_decimal_time corrected_rising + iterate[0]
          corrected_transit = rationalize_decimal_time corrected_transit + iterate[1]
          corrected_setting = rationalize_decimal_time corrected_setting + iterate[2]
        end
        [corrected_rising, corrected_transit, corrected_setting]
      end