function[v]=shwspeed(E,Ep,stps,mu1, mu2)
[v, imax]=espeed(E, Ep, stps, mu1, mu2);
v1=E*mu1;
v2=E*mu2;
hold off;
plot(E,v1,'b-');hold on; %axis([0 1300 0 5000000])
plot(E,v2,'r-');
plot(E,v,'g-');
title('Ve versus E');
xlabel('E field V/cm'); ylabel('Speed cm/s');
grid on;