Function quartiles 

pub fn quartiles(data: &[f64]) -> (f64, f64, f64)

Returns the lower quartile (Q1), median (Q2), and upper quartile (Q3).

Uses linear interpolation (same method as percentile). Returns (0.0, 0.0, 0.0) for empty input.