clear
dt=0.001
t=0:dt:10;
x=zeros(size(t));
x(1)=1;
alpha=0.2;
for i=2:length(t)
    x(i) = x(i-1) + alpha*dt*x(i-1);
end
subplot(2,1,1);
plot(t,x);
y=exp(alpha*t);
subplot(2,1,2);
plot(t,y);