pub fn cmul(a: (f64, f64), b: (f64, f64)) -> (f64, f64)
Complex multiply: (a+ib)(c+id) = (ac-bd) + i(ad+bc)