load GVg

figure(1)
[hh1,h1,h2]=plotyy(field,GVg(:,1),field,GVg(:,1),'semilogy','plot'); hold on
hold on
set(h1,'linestyle','-','color','k','linewidth',[2]);
set(h2,'linestyle','-','color','k','linewidth',[2]);

set(hh1(2),'xlim',[-0.1 0.1],'ylim',[0 1.5]);
set(hh1(1),'xtick',[-0.1:0.05:0.1],'box','on');
set(hh1(2),'ytick',0:0.3:1.5,'box','on');
set(hh1(1),'xlim',[-0.1 0.1],'ylim',[1e-4 10]);
set(hh1(1),'ytick',[1e-4 1e-3 1e-2 0.1 1 10],'box','on');
set(hh1(1),'position',[0.2 0.2 0.65 0.65])
set(hh1(2),'position',[0.2 0.2 0.65 0.65])

h_xlabel_l=get(hh1(1),'xlabel');
h_ylabel_l=get(hh1(1),'ylabel');
h_ylabel_r=get(hh1(2),'ylabel');
set(hh1(1),'fontsize',[20],'linewidth',[2])
set(hh1(2),'fontsize',[20],'linewidth',[2])

[hh1,h1,h2]=plotyy(field,GVg(:,51),field,GVg(:,51),'semilogy','plot'); hold on

hold on
set(h1,'linestyle','--','color','m','linewidth',[2]);
set(h2,'linestyle','--','color','m','linewidth',[2]);

set(hh1(2),'xlim',[-0.1 0.1],'ylim',[0 1.5]);
set(hh1(1),'xtick',[-0.1:0.05:0.1],'box','on');
set(hh1(2),'ytick',0:0.3:1.5,'box','on');
set(hh1(1),'xlim',[-0.1 0.1],'ylim',[1e-4 10]);
set(hh1(1),'ytick',[1e-4 1e-3 1e-2 0.1 1 10],'box','on');
set(hh1(1),'position',[0.2 0.2 0.65 0.65])
set(hh1(2),'position',[0.2 0.2 0.65 0.65])

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','\epsilon_G [V/A]','fontsize',[24]);
set(h_ylabel_l,'string','G_{D} [e^2/h]','fontsize',[24]);
set(h_ylabel_r,'string','G_{D} [e^2/h]','fontsize',[24]);
set(hh1(1),'fontsize',[20],'linewidth',[2])
set(hh1(2),'fontsize',[20],'linewidth',[2])

%print -dtiff -r300 Fig4a