clear
Le=0.2e-6;
Lg=1.8e-6;
S1=400e-12;
E=12.8;
E0=8.85e-12;
V=25;
P=25e-3;
e0=1.6e-19;
tau=0.27e-12;
mu=200e-4;
hv=2.2784e-019;
R=72;
w=1e11:0.01e12:1e13;
k=(1-(tan(pi*Le/(4*(Le+Lg))))^4)^0.5;
C1=pi*(1+E)*E0*S1
C2=2*(Le+Lg)*log(2*(1+k^0.5)/(1-k^0.5));
C=C1/C2;
S2=400e-12;
G0=(P*e0*tau*mu)/(hv*S2);
P1=(V*G0)^2*R;
P2=2*(1+(w.*tau).^2).*(1+(w.*R*C).^2);
y=P1./P2;
y1=1e9*y;
x=w./(2*pi*1e12);
Le=0.5e-6;
Lg=1.8e-6;
S1=400e-12;
E=12.8;
E0=8.85e-12;
V=25;
P=25e-3;
e0=1.6e-19;
tau=0.27e-12;
mu=200e-4;
hv=2.2784e-019;
R=72;
w=1e11:0.01e12:1e13;
k=(1-(tan(pi*Le/(4*(Le+Lg))))^4)^0.5;
C1=pi*(1+E)*E0*S1
C2=2*(Le+Lg)*log(2*(1+k^0.5)/(1-k^0.5));
C=C1/C2;
S2=400e-12;
G0=(P*e0*tau*mu)/(hv*S2);
P1=(V*G0)^2*R;
P2=2*(1+(w.*tau).^2).*(1+(w.*R*C).^2);
y=P1./P2;
y2=1e9*y;
x=w./(2*pi*1e12);
Le=1.0e-6;
Lg=1.8e-6;
S1=400e-12;
E=12.8;
E0=8.85e-12;
V=25;
P=25e-3;
e0=1.6e-19;
tau=0.27e-12;
mu=200e-4;
hv=2.2784e-019;
R=72;
w=1e11:0.01e12:1e13;
k=(1-(tan(pi*Le/(4*(Le+Lg))))^4)^0.5;
C1=pi*(1+E)*E0*S1
C2=2*(Le+Lg)*log(2*(1+k^0.5)/(1-k^0.5));
C=C1/C2;
S2=400e-12;
G0=(P*e0*tau*mu)/(hv*S2);
P1=(V*G0)^2*R;
P2=2*(1+(w.*tau).^2).*(1+(w.*R*C).^2);
y=P1./P2;
y3=1e9*y;
x=w./(2*pi*1e12);
loglog(x,y1,'r',x,y2,'g',x,y3,'k')
set(gca,'linewidth',1);
axis([0.06 6 0.01 5])
legend('\fontsize{12}\rm le=0.2','\fontsize{12}\rm le=0.5','\fontsize{12}\rm le=1.0')
xlabel('\fontsize{12}\rm frequentcy (THz) ')
ylabel('\fontsize{12}\rm output power(nw)')
