Z=0.5 m=1. b=0.2 beta=80. f(x)=(1.-1./Z)*x+b g(x)=(1.-m)*x+b set xrange[0:30./beta] p "sig.inp" u 1:7 w lp lw 2 fit[0:15./beta] f(x) "sig.inp" u 1:7 via Z,b fit[0:15./beta] g(x) "sig.inp" u 1:7 via m,b rep f(x)