new;
clear;
clf;

UENERGY_IN_JOULE = 1.602192E-19;
ULENGTH_IN_METER = 1.E-10;
UMASS_IN_KG = 1.66053E-27;
UTIME_IN_SECOND = ULENGTH_IN_METER*sqrt(UMASS_IN_KG/UENERGY_IN_JOULE);
UFREQ_IN_HZ = 1./UTIME_IN_SECOND/2/pi;
UFREQ_IN_THZ = UFREQ_IN_HZ / 1E12;

load a.out;
item = a(1,1);
n3 = max(a(:,2));
B(n3,n3) = 0;

for k = 2:item+1,
 B(a(k,1),a(k,2)) = a(k,3)+a(k,4)*i;
end;

w = sqrt(eig(B))*UFREQ_IN_THZ;

title ('Vibrational Density of States of SiC ')
hist(w,200);
xlabel ('THz');

PRINT;
