  We study the impedance spectra of woodwind instruments with arbitrary
axisymmetric geometry. We perform piecewise interpolations of the instruments'
profile, using interpolating functions amenable to analytic solutions of the
Webster equation. Our algorithm optimizes on the choice of such functions,
while ensuring compatibility of wavefronts at the joining points. Employing a
standard mathematical model of a single-reed mouthpiece as well as the
time-domain reflection function, which we derive from our impedance results, we
solve the Schumacher equation for the pressure evolution in time. We make
analytic checks that, despite the nonlinearity in the reed model and in the
evolution equation, solutions are unique and singularity-free.
