function vt = threshold(vsb, type, vt0, gamma)
% analysis of threshold in the presence of body effect
phi = -0.6 * type;
vt = vt0 + gamma * (sqrt(abs(-phi + vsb)) - sqrt(abs(phi)));