A symmetric‐matrix time integration scheme for the efficient solution of advection‐dispersion problems