function[v, maxe]=espeed(E, Ep, stps, mu1, mu2)
n1_=1./(70*exp(-Ep./E)+1);
n2_=70*n1_.*exp(-Ep./E);
v=E.*(n1_*mu1+n2_*mu2);
maxe=0