A conceptually new approximation method to study the time-frequency properties of dynamical systems characterized by linear ordinary differential equations is presented. We bypass solving the differential equation governing the motion by writing the exact Wigner distribution corresponding to the solution of the differential equation. The resulting equation is a partial differential equation in time and frequency. We then show how it lends itself to effective approximation methods because in the time frequency plane there is a high degree of localization of the signal. Numerical examples are given and compared to exact solutions.