FEM(Finite Element Method) is carried out to investigate the eigen-properties of cylindrical HAPCF(Hole-Assisted Photonic Crystal Fiber). HAPCF is made of a conventional structure with air holes which are arranged symmetrically with respect to the core in the core-cladding interface. The eigen-mode consistent with the boundary condition of HAPCF is calculated numerically by applying FEM to the vector Helmholtz equation. In this process, the high order matrix equation is resulted and the Krylov-Schur iteration method is applied to obtain the eigen-mode. As a result, the several eigen-modes which reveal the prominent physical properties of HAPCF are obtained and by them the characteristics of the eigen-properties can be understood. In order to promote the physical understanding of the eigen-modes, the spacial distribution of transverse electric field, their corresponding magnitude of the longitudinal electric field and 2 dimensional projection of them are represented by pictorially.