clf;
load freq.out;

plot (freq(:,1), freq(:,2), '.');
xlabel ('frequency number n');
ylabel ('|J(\omega)|^2/2T_f [W/m/K]');
axis ([min(freq(:,1)) max(freq(:,1)) 0 max(freq(:,2))]);

smooth = filtfilt(ones(1,16)/16, [1], freq(:,2));
hold on;
plot (freq(:,1), smooth, 'r');

CHGX; CHGY;
CPRINT; 
