% generate a postsript figure of a saturation curve

global a c k s f g ;

# acks.txt has only numbers (in text, one per line, no comments)
a = 3.1661413854234736753  # "a"
c = 0.19814705422872894158         # "c"
k = 2.0000000000000004441  # "k"            # k not used
s=  1.0307746920239901911  # "s"

f = [0:.0001:1];
g = (s*f.*k^(a*c))./(f.^(1/c)*(k^(a)-1)+s^(1/c)).^c;
gset xrange[0:1]
gset yrange[0:1]
gset size square
%gset key off
gset nokey
xlabel('f(q)');
ylabel('g(f(q)) = f(kq)');
plot(f,g)
%replot

%gset terminal latex  
%gset output 'pork_fig_acs.tex'
%replot

gset terminal postscript eps 25
gset output 'fig_acs.eps'
replot

