  We present the first formulation of the optimal polynomial approximation of
the solution of linear non-autonomous systems of ODEs in the framework of the
so-called $\star$-product. This product is the basis of new approaches for the
solution of such ODEs, both in the analytical and the numerical sense. The
paper shows how to formally state the problem and derives upper bounds for its
error.
