A type of shock fitting method is used to solve some two and three dimensional flowproblems with interactions of various discontinuities. The numerical results show that highaccuracy is achieved for the flow field, especially at the discontinuities. Comparisons withthe Lax-Friedrichs scheme and the ENO scheme confirm the accuracy of the method.