A Superconducting Rotary Machine (SRM) is characterized by an air-cored machine with its rotor iron and stator iron teeth removed. For this reason, the SRM is featured by 3D magnetic flux distribution, which decreases in the direction of axis. Therefore, 3D magnetic field analysis method is required to know about characteristic of magnetic field distribution of SRM. In this paper, 3D flux distribution of SRM is calculated by analytical method. The magnetic field distribution of the field coils is calculated by Biot-Savart equation. The magnetic core is represented by magnetic surface polarities. This paper describes the combined use of above methods for the total field computation, and compares results of analytical method and 3D FEM(Finite Element Method).