Operator technique has been applied for calculation of the quadrupole moment matrix-elements. Master formulas for the quadrupole moment matrix elements for pairs of Slater type, orbitals are derived, one using the expansion method for spherical harmonics and the other the transformed of the quadrupole moment matrix elements into overlap integrals for Mulliken. The numerical values of the quadrupole moment matrix elements evaluated by two methods are in agreement with each other and the calculated quadrupole moment for the ground state of HCl molecule is also in agreement with that of Nesbet.