real T = 300;
real p = 3.7e-10;
real WLC(real x, real L, real p=p, real T=T, real kB=kB,
- real xLmax=0.99, real fmax=fnmax) {
+ real xLmax=0.99, real fmax=fnmax) {
real f = infinity;
if (x < 0)
return infinity;
}
void graphSawtooth(string file="datafile", real k,
- int xcol=0, int fcol=1,
- real xscale=1, real fscale=1, real dx=0, real df=0,
- pen p=red) {
+ int xcol=0, int fcol=1,
+ real xscale=1, real fscale=1, real dx=0, real df=0,
+ pen p=red) {
file fin=input(file).line();
real[][] a=fin.dimension(0,0);
a=transpose(a);
real[] f=a[fcol];
x = x - f/k; /* Remove cantilever extension */
graphData(x=x, y=f, xscale=xscale, yscale=fscale,
- dx=dx, dy=df, p=p, dots=false);
+ dx=dx, dy=df, p=p, dots=false);
}
/* expt.dat already scaled to nm and pN */