x=[0 27.2656 54.5311 81.7765 109.022 136.287 149.91 163.533 177.155 190.798 211.333 224.956 245.511 259.134 279.669 300.224 320.759 341.315 368.762 396.23 423.698 451.165 485.545 513.013 554.326 595.638 643.884 692.109 740.374 795.552 857.683 919.793];
y=[0.4 0.377767 0.355534 0.332132 0.308729 0.286497 0.274795 0.263094 0.251393 0.240861 0.229157 0.217456 0.206922 0.195221 0.183517 0.172982 0.161278 0.150744 0.139038 0.1285 0.117963 0.107426 0.0968867 0.0863497 0.0769769 0.067604 0.0593981 0.0500226 0.0429862 0.0359472 0.0312446 0.0253724];
f=@(c,x)c(1)*log(1+exp((c(2)-x)/c(3)));
c0=[1 -10 1];
[c,resnorm]=lsqcurvefit(f,c0,x,y);