function X = euler (f, x0, dt, T)
% Euler-Verfahren
% Input: <spaeter>

N = ceil(T/dt);
dt = T/N;

X = zeros(size(x0,1), N+1);

X(:,1) = x0;

for k = 1:N
    uk = f(X(:,k));
    X(:,k+1) = X(:,k) + dt*uk;
end