The method of Bishop and Randič for evaluating force constants by direct differentiation and using only the ground-state wavefunction for the molecule in its equilibrium configuration is extended so as to handle wavefunctions built from a nonorthogonal basis set. The results are applied to a series of Gaussian wavefunctions for H2+.