f(x)=a+b*x
fit f(x) "<(head -n 20 data.dat && tail -n 21 data.dat)" via a,b
plot "data.dat" w p, a+b*x w l