function F=taub(v)

global Tb

F=Tb*( 1.24 + 2.678 ./(1+exp((v+50)./16.027)) )  ;