A numerical algorithm is presented for solving laminar, steady, supersonic interacting boundary-layer flows for quasi-three-dimensional configurations. The interaction problem is treated as a boundary-value problem and a salient feature of the scheme is the direct implementation of the downstream boundary condition. Solutions are presented for axisymmetric and swept (yawed) compression ramps for both adiabatic and heat transfer conditions over a Mach number range of 2–6. The results are in good agreement with experimental data and existing theories for axisymmetric cases. For the swept (yawed) configurations, lack of experimental data makes a direct comparison impossible, but the present solutions are found to be in qualitative agreement with earlier studies. In addition, it is shown that the trends obtained for the sweep effects are well predicted by a simple extension of a two-dimensional asymptotic theory.