k=[-18:.01:-3];
x=10. k;
a=((exp(x)-1)./x);plot(a);
b=((exp(x)-1)./(log(exp(x))));
plot(k,a,r');hold on
plot(k,b,'b');hold on