A viscous-inviscid interaction method to calculate the subsonic and transonic flow over nozzle afterbodies with supersonic jet exhausts was developed. The method iteratively combines a relaxation solution of the full potential equation for the inviscid external flow, a shock capturing-shock fitting inviscid jet solution, an integral boundary layer solution, a control volume method for treating separated flows, and an overlaid mixing layer solution. A computer program called RAXJET which incorporates the method, illustrates the predictive capabilities of the method by comparison with experimental data is described, a user's guide to the computer program is provided. The method accurately predicts afterbody pressures, drag, and flow field properties for attached and separated flows for which no shock induced separation occurs.