A Eulerian finite difference technique is presented for solving the conservation equations for homogeneous compressible fluids. The method uses a Newton-Raphson-type relaxation scheme formed to enable the computation of sharp density gradients at all flow speeds. A derivation of the numerical algorithm and two example problems are presented.