Implement in Matlab the rotary loudspeaker effect discussed on page $23$ of the session $7$

Lecture notes $($also on page $93$ of Zoelzer et al$,$ and in Disch et al$).$ Write this as a Matlab

function that can be called as follows