pub fn interval_find_zeros(
    c1: Float,
    c2: Float,
    c3: Float,
    c4: Float,
    c5: Float,
    theta: Float,
    t_interval: Interval,
    zeros: &mut [Float; 8],
    zero_count: &mut u8,
    depth: usize
)