fn dst1(input: &[f64], sin_table: &[Vec<f64>], output: &mut [f64])
1D DST-I: X[k] = Σ x[n] * sin(π*(n+1)*(k+1)/(N+1))