A low dispersion, high accuracy finite element method for first order hyperbolic systems in several space variables