The computation of slowly moving Shockwaves is considered. A numerical error is introduced into the solution during an explicit computation, which results in the generation of spurious entropy. The generation of this spurious entropy and its effect on the pressure wave in Euler flows is examined, in one dimension by considering a moving shock in a constant area duct, and in two dimensions by considering an unsteady aerofoil flow. An efficient temporal discretisation is presented, whose advantages over the explicit scheme are two-fold. First, the numerical entropy generation due to shock motion is eliminated, and second, the CPU requirements for a two-dimensional unsteady computation are reduced by a factor of five from that required for the explicit scheme, for no loss of accuracy.