The system for producing magnetic field constitutes an important component of magnetic refrigerator. Many researchers have directed significant effort to increase the magnetic field intensity, because the magnetocaloric effect at the Curie temperature increases with the power of 2/3 of the magnetic field. In this study, we report the simulation of the magnetic field intensity at polar axis of a Halbach cylinder (HC) by i) changing the length and thickness of the HC, ii) having with or without gap of the HC, and iii) surrounding the HC with a soft magnet shell, acting as a shielding. We simulated the field distribution of a HC with a finite size. Furthermore, the detailed numerical results of the magnetic field distribution and its dependence on shielding are presented in this study.