n=10;T=5;k=T/n;
t=linspace(0,T,n+1);
u=zeros(1,n+1);
u(1)=1;
for i=1:n
    u(i+1)=u(i)-4.2*k*u(i);
end
plot(t,u)