clim = 25; cA0 = 1/exp(1);
r1=0.80:0.01:1.00;
n=length(r1);
St = @(cAclim,r1) log((1/(cAclim*r1)-(1-r1)/r1))-1/(1-(1-r1)*cAclim);
caclim = arrayfun(@(i) fzero(@(cAclim) St(cAclim,r1(i)),cA0),1:n);
plot(r1,caclim*exp(1), 'k',1, 1,'k*')
v=axis; axis([0.95*v(1), 1.05*v(2), 0.99*v(3), v(4)])
xlabel('R1'); ylabel('c1')

set(gcf, 'PaperUnits', 'centimeters') 
set(gcf, 'PaperPosition', [3 10 11 8.25])
print -depsc 'fol99r1c.eps'

