In this paper, a method for numerically solving the equations of a three-dimensional viscosityless transonic flow has been developed, based on the finite volume method. Stability and rapid convergence is ensured by the introduction of artificial viscosity terms of the second and fourth orders. The system of differential equations is solved with the help of approximative "LU" implicit factorization with Jameson-Yune flow separation.