A Wigner Distribution Approximation Method for Ordinary Differential Equations