Based upon Wu's general three-dimensional turbomachine flow theory, a set of computer programs was established to measure flow along the relative stream surfaces S1 and S2, employing nonorthogonal curvilinear coordinates and nonorthogonal velocity components. The streamline extension technique, matrix direct solution, and line relaxation method are used in these programs. With the employment of fictitious grid points, the internal storage capacity required in the matrix direct solution is reduced by nearly two-thirds. The matrix direct solution and the line relaxation method are found to converge nearly 10 times faster than the streamline extension technique. These programs constitute a design system of axial flow turbomachine blades based upon three-dimensional flow. Experience shows that iterative solutions considering three-dimensional flow obtained by the successive calculation of the Sj and S2 surface flows are rapid and practical. The major results of a number of typical applications are included in the paper.