clear, set(0,'DefaultFigureWindowStyle','docked')
N = 200000; h = 1/N;
Eul = nan(1, N+1);
y0 = 1; Eul(1) = y0;

for ii = 1:N
    Eul(ii+1) = Eul(ii) + h*(atan(Eul(ii)));
end

t = linspace(0, 1, N+1);figure(1);
plot(t, Eul, '*-')
legend({'Eul'},'location','northwestoutside','FontSize',24)