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