x=a:.01:b;
x=x(:);
y=exp(-x).*x;
plot(x,y);
shg