  We propose a functional implementation of \emph{Multivariate Tower Automatic
Differentiation}. Our implementation is intended to be used in implementing
$C^\infty$-structure computation of an arbitrary Weil algebra, which we
discussed in the previous work.
