Function sine_power 

fn sine_power(w: f64, w_low: f64, order: i32) -> Scaled

w^{μ/2} for 1 - y² = w + w_low (0 < w <= 1).