#calculate the Vibrational degrees of freedom #initialization of variables atoms=5 #calculations f=3*atoms fvib=f-3-3 #results print '%s %d' %("Vibrational degrees of freedom = ",fvib)
Vibrational degrees of freedom = 9