1 2 3 4 5 6 7 8
clc //initialization of variables atoms=5 //calculations f=3*atoms fvib=f-3-3 //results printf("Vibrational degrees of freedom = %d",fvib)