Born-Oppenheimer breakdown effects in the determination of diatomic internuclear potentials: Application of a least-squares fitting procedure to the HCl molecule