f = dlmread("modes.dat",',');
plot(f(:,1),f(:,2:end),'b-',f(:,1),f(:,1),'k-');
xlabel("wavevector kx (units of 2\pi/a)");
ylabel("frequency (units of 2\pic/a)");
axis([0 0.5 0 0.4]);
