# fit relative unit conductance to a Hill type relationship
k_half = 0.002
H = 1
f(x) = 1 / (1 + ((x/k_half)**H))
set logscale x
set xrange [1e-6:1e-2]
plot f(x),"IK1-pH_i-fit.data" using 2:3:4 w yerrorbars
FIT_LIMIT = 1e-6
fit f(x) "IK1-pH_i-fit.data" using 2:3 via k_half,H