clear, set(0,'DefaultFigureWindowStyle','docked')
N = 10; h = 1/N; lambda = -2; %modify these parameters to experiment
expEul = nan(1, N+1); impEul = nan(1, N+1); 
y0 = 1; expEul(1) = y0; impEul(1) = y0; 
%impMpr(1) = y0;

for ii = 1:N
    expEul(ii+1) = expEul(ii)*(1+h*lambda);
    impEul(ii+1) = impEul(ii)/(1-h*lambda);
end

t = linspace(0, 1, N+1);figure(1); 
plot(t, exp(lambda*t), t, impEul, '*-', t, expEul, 'k*-')
legend({'exact', 'impEul','expEul'},'location','northwestoutside','FontSize',24)