% Explizite Eulersche Methode

clear

y0=1;
T=0.7;
n=500;
tau=T/n;

t=[0:tau:T];
[a b]=size(t);

Y(1)=y0;
for k=2:b
    Y(k)=Y(k-1)+tau*f(Y(k-1));
end

% plot
plot(t,Y,'k')

