Calculating program for compressible flows with area changes