#closeall;
nonorm;

f = getdata('monitor1','f');
Py=getdata("monitor1","Py");
x=getdata("monitor1","x");
T2 = 0.5*integrate(real(Py),1,x);
plot(f/1e12,log10(T2),'frequency','transmission');