  We consider linear dynamical systems with quadratic output. We first define
the two transfer functions, a single-variable and a multivariate one, that
fully describe the dynamics of these special nonlinear systems. Then, using the
samples of these two transfer functions, we extend the AAA algorithm to model
linear systems with quadratic output from data.
