clear
reset
set term qt
set title "P(p), error-bars and gaussian-fit"
plot [-2.5:6.0] [0:] "AP0-0-80000.8_2000_15" lw 1.5 lc 7 w l
replot "AP0-0-80000.8_2000_15" every 30 lc 7lw 1.5  w e
f(x)=exp(-a*(x-1)**2)
fit [-2.5:6] f(x) "AP0-0-80000.8_2000_15" via a
replot f(x) lc "black"
