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)