static double cos_squareWave(double phi) {
  /*double phase = frac(phi/(pi()*2));
  ret phase < .5 ? 1 : -1;*/
  todo();
}