m=1e-3;
t=-2*m:0.01*m:2*m;

a=0;

for i=1:4


a=a+250;

xt=(20*sin((2*pi*1000.*t)-(pi/3))).*(exp(-1.*a.*t));

subplot(4,1,i)
plot(t,xt)

end