I'm using this transfer function within a flute waveguide. Here's the excerpt about the transfer function.
"where f is the value of the volume flow advancing towards the
embouchure hole, k is a gain term, and h describes the deviation
of the center line of the air jet from the wedge. By changing the value of h it is possible to produce asymmetrical wavefom or- in other words - sounds that also have even harmonics in their spectrum. When h is zero, the spectrum of the output of the model has very weak even harmonics"
Admittedly the following excerpt is about another slightly different function (see screenshot), but the idea is still the same.
f(x) = -1 if x<= -1.325
x^3-x if 1.325<x<-1.325
1 if x>=1.325
"If the input signal does not go beyond (-1.325,1.325) range the waveshaper inverts the phase of the input signal and introduces a third harmonic. Once the range is exceeded, hard limiting occurs and inserts odd harmonics from clipping."
http://www.pdpatchrepo.info/hurleur/Screen_shot_2013-07-12_at_10.31.44.png