In the treatment of Heisenberg Hamiltonian containing single-ion anisotropy by means of many-body Green's function method, the decoupling of the anisotropy term in higher order Green's functions usually takes an Anderson-Callen ͑AC͒ form. In this paper, possible improvement schemes of AC decoupling are discussed by comparison of the results with those of quantum Monte Carlo and frame rotation methods. We choose one scheme with a factor concerning the effect of the external field. Next, we discuss a possible difficulty of the frame rotation method in treating systems with single-ion anisotropies in more than one direction. Then we extended our method ͓Phys. Rev. B 70, 134424 ͑2004͔͒ to treat magnetic films. Some magnetic properties of ultrathin ferromagnetic films with thicknesses up to 16 monolayers are studied. The properties investigated include transition point, effective anisotropy coefficient, field-induced magnetization reorientation, and hysteresis loop. Several cases are investigated for uniaxial anisotropy and external field along different directions. The transition point, the effective anisotropy coefficient, the coercivity, and the loop area increase with increasing film thickness. The coercivity decreases and the loop area reduces with increasing temperature. The hysteresis loop along field direction is different from that along easy-axis direction. The coercivity and the area of the loop obtained in the former case are larger than the latter case. The reasons are analyzed by investigation of the trajectory of magnetization in detail. The influence of dipole interaction on the magnetic properties is discussed.