Numerical determination of flux components in vapor flow considers split boundary conditions in both directions to form grid lines for flow continuity equations. Matrix of coefficients obtains direct solution for potential flow problem without iteration.