%% Explizites Euler-Verfahren zur Lösung von ODE (-Systeme) 1. Ordnung 
function [t,z]=ode_euler(odefunction,t,z0)
%% Anfangswerte, falls es sich um ein System handelt, 
% ist z0 ein Vektor, und wird in die erste Zeile von z geschrieben
z(1,:)=z0;
%% Schleife über das Integrationsintervall
% am Ende der Schleife enthält (im Falle eines Systems)
% jede Spalte von z eine Lösung: 1. Spalte z1, 2. Spalte z2 ...
for ii=2:(length(t));
    h=(t(ii)-t(ii-1));
     % F kann auch ein Vektor sein 
    F=odefunction(t(ii-1),z(ii-1,:)');
    z(ii,:)=z(ii-1,:)+h*F';
end
%%
end
