The solution of systems of non-autonomous linear ordinary differential equations is crucial in various applications, such as nuclear magnetic resonance spectroscopy. We introduced a new solution expression in terms of a generalization of the Volterra composition. Such an expression is linear in a particular algebraic structure of distributions, which can be mapped onto a subalgebra of infinite matrices. It is possible to exploit the new expression to devise fast numerical methods for linear non-autonomous ODEs. As a first example, we present a new method for the operator solution of the generalized Rosen-Zener model, a system of linear non-autonomous ODEs from quantum mechanics. The new method’s computing time scales linearly with the model’s size in the numerical experiments. A second example is the analysis of temporal network, where the new expression might lead to novel extension of subgraph centrality indexes. Streaming link (MS Teams): https://teams.microsoft.com/l/meetup-join/19%3ad48916dee7364e0bbb08216ea621de47%40thread.tacv2/1738144836095?context=%7b%22Tid%22%3a%2224c5be2a-d764-40c5-9975-82d08ae47d0e%22%2c%22Oid%22%3a%226f4c2a62-3ec4-4b20-9af3-43b578f2b18e%22%7d