Function clenshaw 

pub fn clenshaw(szeta: f64, czeta: f64, f: &[f64], k: i32) -> f64

Evaluate y = sum(F[k] * sin((2*k+2) * zeta), k, 0, K-1) by Clenshaw summation. zeta is specify by its sine and cosine, szeta and czeta.