Matlab - graficar la funcion de distribución de fermi-dirac

Hola : la verdad soy en terrible en matlab y necesito su ayuda necesito graficar esta función en matlab,
f=1/(e^(E-u)/KT)
donde:
K=1.38*10^-23

debo graficar para T=0 , T=300, T=1000, T=5000
u=cte y E es la variable independiente.

Les agradeceria mucho su ayuda

% define constantes
K = 1.38e-23;
u = 1;

% T1 es 'eps' para que no haya división por cero
T1 = eps;
T2 = 300;
T3 = 1000;
T4 = 5000;

% define valores de variable independiente
E = -15 : .1 : -10;

% define tus funciones paramétricas
f1 = 1./(exp(E-u)/(K*T1));
f2 = 1./(exp(E-u)/(K*T2));
f3 = 1./(exp(E-u)/(K*T3));
f4 = 1./(exp(E-u)/(K*T4));

% grafica las líneas
plot(E,f1,'r', E,f2,'g', E,f3,'b', E,f4,'m')
legend('T=0', 'T=300', 'T=1000', 'T=5000')

Este código ocupe yo para T= 0TF ; T=0.1 TF ; T=0.2 TF y T=0.5TF
Es por eso que dejé K=1 (TF: temperatura de Fermin.
Y con valores de relación entre la energía y la energía de Fermin que van de 0 a 1.5 (valores de variable independiente) y con saltos de 0.001

hola buenas tardes, disculpa que te moleste pero no entiendo porque consideraste a k =1? podrias explicarme ??