tex/src/figures/fit-space/king_vs_best.asy

   1 import wtk_graph;
   2
   3 size(15cm,10cm,IgnoreAspect);
   4
   5 scale(Log, Linear);
   6 real xscale=1;
   7 real yscale=1;
   8
   9 /* f(x) = A + log10(x) + B */
  10 real fn_logxliny(real x, real[] params) {
  11   return params[0] * log10(x) + params[1];
  12 }
  13
  14 real xmin = 1.01457e-05;
  15 real xmax = 0.10364;
  16
  17 graphFile("Best_2002_detailed_unfolding_pathway/fig3a.dat", xcol=1, ycol=0,
  18           xscale=xscale, yscale=yscale, p=psoft, mpath=m30,
  19           t="Best valley", dots=true);
  20 graphFile("fit-valley.dat", xscale=xscale, yscale=yscale, p=phard, mpath=m30,
  21           t="King valley", dots=true);
  22 fitFile("Best_2002_detailed_unfolding_pathway/fig3a.dat.fit.dat",f=fn_logxliny,
  23         xmin=xmin, xmax=xmax, xscale=xscale, yscale=yscale, p=psoft);
  24 fitFile("fit-valley.dat.fit.dat", f=fn_logxliny,
  25         xmin=xmin, xmax=xmax, xscale=xscale, yscale=yscale, p=phard);
  26
  27 label(sLabel(""), point(N), N);
  28 xaxis(sLabel("$k_{u0}$ ($1/$s)"), BottomTop, LeftTicks);
  29 yaxis(sLabel("$x_u$ (m)"), LeftRight, RightTicks);
  30
  31 add(legend(), point(E), 20E, UnFill);