A systematic method to fit calculated to observed vibrational frequencies has been developed and implemented in a computer program. The procedure consists of the refinement of a scaled quantum mechanical force field (SQMFF) previously obtained according to Pulay's method. The key step in the process is the generation of an intermediate matrix, C LambdaC(T), which is then refined. The above step produces only small corrections to the scaled force constants, yielding a considerable improvement of the fitted frequencies. This scheme of refinement can be carried out using any kind of coordinates. To show the reliability and performance of the proposed method, the force fields of two very different systems, as benzene and tetranitromethane, have been chosen as example tests.