A computer code to calculate two-dimensional or axisymmetric compressible internal flows is developed based on LU-SGS algorithm. Flows of a quasi-one dimensional convergent-divergent nozzle (Denton nozzle), a hyperbolic slender nozzle, a ramjet nozzle, and the JPL nozzle are computed using the code. Computational results for various cases are compared with the analytic or experimental data available. The comparison shows that the code works properly.