% Intrinsic carrier concentration, T is in Kelvin
function result = ni(T)

E = 8.617e-5 * T;
result = 3.1e16 * T^(3/2) * exp(-0.603/E);