pub fn fourier_series_evaluate(
    coefficients: &[SeriesCoefficient],
    period: f64,
    t: f64
) -> f64