%  Euler framåt för differentialekvationer i flera variabler
% 
%
%
y=[1;1]

b=1;
a=0;
N=400;
h=(b-a)/N;

t=0;
yv1=[];
tv=[];
% Framåt Euler
for i=1:N
    % För ODE-system är det bäst att göra funktionsfiler (inte @-notation)
    y=y+h*mittsystem_f(t,y);   

    % Såhär får man inte göra. Dessa två rader är ett typiskt fel 
    % studenter brukar göra.
    % Problem: Man får inte uppdatera en ekvation i taget eftersom
    % man måsta använda det gamla y(1)-värdet
    %y(1)=y(1) + h*(3*y(1)+y(2)*y(2));
    %y(2)=y(2) + h*(y(1)-4*y(2));

    t=t+h;
    yv1=[yv1;y(1)];
    tv=[tv;t];
end
hold on;
plot(tv,yv1,'r--');
