use: kb=1.38e-23 q=1.6e-19 iref=1e-6 C=1e-9 delta_u=5.0 delta_t=1e-7 t=273+27 #function for Ut, temperature is x! u(x)=kb*x/q #expo function i(x,t)=iref*exp(x/u(t)) #time proportional to C and delta u (low/high level) , inverse to i plus some delay t(x,t)=C/i(x,t)*delta_u+delta_t #frequency f(x,t)=1.0/t(x,t) m.c.