function Figure1_Best_Approximant()
%
%  This script plots *only the best* approximant for each of the four
%  parameter sets in Figure 1
%
% Closed-form approximations for stationary single-channel Ca2+ nanodomains
%          in the presence of buffers with two binding sites
%
%         Victor Matveev and Yinbo Chen (NJIT) March 2021
% (GPL) Free to use, share, modify, with appropriate attribution
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Names = ["RBA", "ExpPadeA", "ExpPadeB", "ExpExp"];
% Method's legend string to be added to the last panel of each figure

Params = [8,      10,    1,    1;
          8,      10,  0.1,  0.1;
          0.003, 100,    1,    1;
          0.003, 100,  0.1,  0.1 ];

for k = 1 : size(Params, 1)

  lambda2 = Params(k, 1);
  nu2     = Params(k, 2);
  eps     = Params(k, 3);
  gamma   = Params(k, 4);
  
  method = Select_Best_Buffer_Approximant(lambda2, nu2, eps, gamma);
  figure;
  Get_ApproximantByName_ZeroRestingCa(method, lambda2, nu2, eps, gamma, 'b-');
  legend(Names(method));

end
