% Fig.2

clear;
close all;

V=0*1.0e-3;
delay2=0.062;
v_sr1=0.3*1.0e-3;
v_sr2=0.75*1.0e-3;
v_sr3=1.17*1.0e-3;
v_sr4=1.55*1.0e-3;







% v_sr=1.0e-3; 
v_es=1.799e-3;
% delay2=0.06;
v_p3xi=0.101e-3;
v_p3e=0.8e-3;
% v_sr=1*1.0e-3;
v_p1xi=3.01e-4;
open1=1;
open2=1;
KK=1;


open3=1;open4=1;
AA=1;



[tt,xx1]=BGCT_subfun(delay2,v_sr1,v_p1xi,open1,open2,KK,v_p3xi,open3,open4,AA,v_es,v_p3e,V);
[tt,xx2]=BGCT_subfun(delay2,v_sr2,v_p1xi,open1,open2,KK,v_p3xi,open3,open4,AA,v_es,v_p3e,V);
[tt,xx3]=BGCT_subfun(delay2,v_sr3,v_p1xi,open1,open2,KK,v_p3xi,open3,open4,AA,v_es,v_p3e,V);
[tt,xx4]=BGCT_subfun(delay2,v_sr4,v_p1xi,open1,open2,KK,v_p3xi,open3,open4,AA,v_es,v_p3e,V);

figure(1),
subplot(221),plot(tt,xx1),set(gca,'YDir');
subplot(222),plot(tt,xx2),set(gca,'YDir');
subplot(223),plot(tt,xx3),set(gca,'YDir');
subplot(224),plot(tt,xx4),set(gca,'YDir');