Gauss has articulated some notion of recursion in his work on what we today call the Cooley-Tukey Fast Fourier Transform (FFT) algorithm.
where one finds the following translated passage of Gauss's original:
And so for this case, where most of the proposed values of the
function X, an Integral period of the arrangement, the number is
composite and $=\pi=\mu\nu$, in articles 25, 26, we learned that
through the division of that period into $\nu$ periods of $\mu$ terms, it
produces, when all values are given, the same satisfactory function,
which by the immediate application of the general theory applies to
the whole period; truly, that method greatly reduces the tediousness
of mechanical calculations, success will teach the one who tries it.
Gauss essentially describes recursion here, though the language is not explicit. That you can use the solution of a sub-problem of the same form to compute the whole was known to Gauss because indeed he would compute examples in the paper. So while his language to describe the phenomenon is indirect (i.e. there isn't a clearly formed concept of "recursion") it is understood to be part of the computation. Gauss notes that in this case this setup "greatly reduces the tediousness of mechanical calculations" which is an imprecise statement about the computational complexity of the FFT.
This work of Gauss was not published during his life-time and is generally dated around 1805.