Numerical solution of systems of quasilinear hyperbolic differential equations by means of the method of nebencharacteristics in combination with extrapolation methods