Function phase_spectrum 

pub fn phase_spectrum(xs: &[f64]) -> Vec<f64>

Compute the phase spectrum of a real signal.

Returns arg(X[k]) for k = 0..N/2+1.