function lfo(phase) (
  phase = (phase - floor(phase)) * 4; // phase is 0..4
  phase < 1 ?  1 - (1 - phase)*(1 - phase)*(1 - phase) : // 0..1
  phase < 2 ?  1 - (phase - 1)*(phase - 1)*(phase - 1) : // 1..0
  phase < 3 ? -1 + (3 - phase)*(3 - phase)*(3 - phase) : // 0..-1
              -1 + (phase - 3)*(phase - 3)*(phase - 3) ; // -1..0
);