%A standard Publication quality matlab plot file 
%Please change the vectors loaded
%Please change the axis limits according to your data
%Please comment out unnecessary lines

% load id_vg_cl;
% load id_vg_qm;
% load id_vd_cl;
% load id_vd_qm;

a=linspace(0,.6,100);
b=a.*a;
c=a.^3;
d=a.^4;


[hh1,h1,h2]=plotyy(a(:,1),a(:,2)*1e-6,c(:,1),c(:,2),'semilogy','plot');
hold on
[hh2,h3,h4]=plotyy(b(:,1),b(:,2)*1e-6,d(:,1),d(:,2),'semilogy','plot');
 
set(h1,'linewidth',[2]);
set(h2,'linewidth',[2]);
set(h3,'linestyle','--','linewidth',[2]);
set(h4,'linestyle','--','linewidth',[2]);

set(hh1(1),'xlim',[0 0.6],'ylim',[1e-9 1e-2]);
set(hh1(1),'xtick',[0:0.1:0.6],'box','off');
set(hh1(1),'ytick',[1e-9 1e-8 1e-7 1e-6 1e-5 1e-4 1e-3 1e-2]);
set(hh1(2),'xlim',[0 0.6],'ylim',[0 2800],'box','on');
set(hh1(2),'xtick',[0:0.1:0.6],'ytick',[0:400:2800]);

set(hh2(1),'xlim',[0 0.6],'ylim',[1e-9 1e-2]);
set(hh2(1),'xtick',[0:0.1:0.6]);
%set(hh2(1),'ytick',[1e-9 1e-8 1e-7 1e-6 1e-5 1e-4 1e-3 1e-2]);
set(hh2(2),'xlim',[0 0.6],'ylim',[0 2800]);
set(hh2(2),'xtick',[0:0.1:0.6],'ytick',[0:400:2800]);

h_xlabel_l=get(hh1(1),'xlabel');
h_ylabel_l=get(hh1(1),'ylabel');
h_ylabel_r=get(hh1(2),'ylabel');
set(h_xlabel_l,'string','V_{GS} or V_{DS} [V]','fontsize',[20]);
set(h_ylabel_l,'string','logI_{DS} [A/\mum]','fontsize',[20]);
set(h_ylabel_r,'string','I_{DS} [\muA/\mum]','fontsize',[20]);
grid off
set(hh1(1),'fontsize',[20],'linewidth',[2])
set(hh1(2),'fontsize',[20],'linewidth',[2])
set(hh2(1),'fontsize',[20],'linewidth',[2])
set(hh2(2),'fontsize',[20],'linewidth',[2])
set(hh1(1),'position',[0.15 0.15 0.70 0.74])
set(hh1(2),'position',[0.15 0.15 0.70 0.74])
set(hh2(1),'position',[0.15 0.15 0.70 0.74])
set(hh2(2),'position',[0.15 0.15 0.70 0.74])
set(h_xlabel_l,'position',[0.299233 1.2e-10 17.3205])
set(h_ylabel_l,'position',[-0.0928645 3.00206e-06 17.3205])
set(h_ylabel_r,'position',[0.70 1390.97 17.3205])

print -djpeg90 fig4.jpg