abstract: We present an adaptation of Stein's method of multivariate normal approximation that is suitable for dynamical systems. Our method is direct in that it avoids "auxiliary randomization" common in the implementation of the method in stochastics. We also discuss some non-autonomous systems modeling non-equilibrium processes in physics.