Euler in Matlab implementieren.

Hallo Forum, ich möchte gängige numerische Verfahren zur Lösung von 
Anfangswertproblemen in Matlab implementiern. Habe aber einige 
Verständnisslücken.
Nehmen wir als einfachen Fall das Euler-Verfahren mit der Vorschrift:
x_k+1 = x_k + T*f(x_k, u_k) x_k ist der aktuelle Wert zum k-ten 
Abtastzeitpunkt und T die Schrittweite der Abtastung. u ist die 
Eingangsgröße

Jetzt zu implementierung. Ich bräuchte doch bestimmt eine Funktion, die 
auf jeden Fall irgendeinen Anfangswert und einen Zeit-Vektor als 
Paramenter hat.

Also zB:
function x = euler(x0, t)
Aber da hörts bei mir auch schon auf und ich möchte auch nicht mit aus 
dem Netz eine Lösung suchen, sondern selbst bzw mit eurer Hilfe auf 
diese kommen ;). Kann mir jemand weiterhelfen?
Ich würde eine for-Schleife vermuten wo man dann immer das aktuelle x 
durch den vorher berechneten Wert für x ersetzt, oder eben am Anfang das 
x0 benutzt.

Vielen Dank!
LG, Franzi