Script for exporting Force constants matrix

[fly]

I am eager to get the script for exporting Force constants matrix after geometry optimization. who can help me?
 

[zh]

The force on each atom can be printed out easily by calling "calculateAtomicForces(self_consistent_calculation = scf)". You can manually collect all of them and construct the force constant matrix using the free code named "PHON"(http://www.homepages.ucl.ac.uk/~ucfbdxa/phon/)