fn dst_sin_table(n: usize) -> Vec<Vec<f64>>
Build sin table for DST-I of given length. sin_table[n][k] = sin(π*(n+1)*(k+1)/(N+1))