We have developed a special computing code for calculation of nuclear quadrupole moments versus deformation parameter β. For some even-even isotopes, it has been seen that by increasing neutron number, deformation parameter also increase, which means more deformation from spherical shape.