We have developed a special computing code for calculation of nuclear shape changes and quadrupole moments ( Q) of Promethium Isotopes. It has been shown from these calculations that by increasing neutron number, deformation parameter also increase for more heavier isotopes which means more deformation from spherical shape. By comparison with Nilsson level energy diagrams we can infer quadrupole deformation parameter (β₂) and hence calculate quadrupole moments of these isotopes.