script to calculate vdos