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