clear;
close all;

len1=31;
len2=31;

v_xip3=linspace(0.0,2.0,len1)*1.0e-3;
v_se=linspace(0.0,4.0,len1)*1.0e-3;
v_sr=0.8e-3; 
v_es=1.799e-3;
delay2=0.06;
% v_p3xi=2.5e-3;  % bian 
v_ep3=0.3*1.0e-3;
% v_sr=1*1.0e-3;
v_p1xi=3.01e-4;
open1=1; 
open2=1;
KK=1;


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





State=zeros(len1,len2);
FD=zeros(len1,len2);

for i=1:len1
    tic
    parfor j=1:len2
        [p1,p2,p3,p4,p5,p6]=BGCT_subfun2(delay2,v_sr,v_p1xi,open1,open2,KK,v_xip3(i),open3,open4,AA,v_es,v_ep3,v_se(j));
        State(i,j)=p5;
        FD(i,j)=p6;
    end
    toc
end

save data

figure(1)
subplot(221),imshow(State,[]),set(gca,'ydir','normal'),colorbar;
subplot(222),imshow(FD,[]),set(gca,'ydir','normal'),colorbar;
colormap(jet);