Numerical solution of hyperbolic systems with different time scales using asymptotic expansions