function Figure1()
%
%   This script reproduces approximants in Figure 1 of the manuscript
%
% 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
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

figure;

lambda2=8; nu2=10; eps=1; gamma=1;

Get_ApproximantByName_ZeroRestingCa('ExpExp',   lambda2, nu2, eps, gamma, 'k--');
Get_ApproximantByName_ZeroRestingCa('ExpPadeB', lambda2, nu2, eps, gamma, 'm-');
Get_ApproximantByName_ZeroRestingCa('ExpPadeA', lambda2, nu2, eps, gamma, 'm--');
legend('ExpExp', 'ExpPadeB', 'ExpPadeA');

figure;

lambda2=8; nu2=10; eps=0.1; gamma=0.1;

Get_ApproximantByName_ZeroRestingCa('ExpExp',   lambda2, nu2, eps, gamma, 'k--');
Get_ApproximantByName_ZeroRestingCa('ExpPadeB', lambda2, nu2, eps, gamma, 'm-');
Get_ApproximantByName_ZeroRestingCa('ExpPadeA', lambda2, nu2, eps, gamma, 'm--');
legend('ExpExp', 'ExpPadeB', 'ExpPadeA');

figure;

lambda2=0.003; nu2=100; eps=1; gamma=1;

Get_ApproximantByName_ZeroRestingCa('ExpExp',   lambda2, nu2, eps, gamma, 'k--');
Get_ApproximantByName_ZeroRestingCa('ExpPadeB', lambda2, nu2, eps, gamma, 'm-');
Get_ApproximantByName_ZeroRestingCa('ExpPadeA', lambda2, nu2, eps, gamma, 'm--');
Get_RBA_AnyRestingCa(nu2, eps, gamma, 0, 'r-');
legend('ExpExp', 'ExpPadeB', 'ExpPadeA', 'RBA');

figure;

lambda2=0.003; nu2=100; eps=0.1; gamma=0.1;

Get_ApproximantByName_ZeroRestingCa('ExpExp',   lambda2, nu2, eps, gamma, 'k--');
Get_ApproximantByName_ZeroRestingCa('ExpPadeB', lambda2, nu2, eps, gamma, 'm-');
Get_ApproximantByName_ZeroRestingCa('ExpPadeA', lambda2, nu2, eps, gamma, 'm--');
Get_RBA_AnyRestingCa(nu2, eps, gamma, 0, 'r-');
legend('ExpExp', 'ExpPadeB', 'ExpPadeA', 'RBA');
