function [tout, yout]= odEul(FunFcn , t0, tfinal, y0, Nstp);
%
% Eulers metod (den är av ordning 1).
%
% FunFcn skall vara en funktionsfil som beräknar förstaderivatorna av y,
% dvs differentialekvationssystemet skrivet på standardform.
%
% t0 och tfinal är start- respektive sluttid.
% y0 är funktionens startvärde.
% Nstp är antalet steg intervallet skall delas upp i.
%
h = (tfinal - t0)/Nstp;        %steglängden beräknas
t = t0;
y = y0(:);                     %y blir en kolumnvektor
tout = t; 
yout = y.';                    %yout blir en radvektor med
                               % samma komponenter som y
                               % .' transponerar utan att komplex-
                               % konjugera(viktigt om y är komplex)        
for i = 1:Nstp,
   s1   = feval(FunFcn,t,y);   %motsvarar s1=FunFun(t,y);
   y    = y + h*s1;            %FunFun måste leverera en kolumnvektor;     
   t    = t + h;
   tout = [tout; t];           %lägg till ett värde i tout
   yout = [yout; y.'];         %lägg till en rad i yout
end;