% Av Erik Dalsryd 2017-10-11
clear all, close all, clc
y = 1;   yPrim = 1;   % Begynnelsev�rden
u = [y  yPrim]';     % Variabelbyte
h = 0.2;       % Stegl�ngd
xStart = 0;   % Start, oberoende variabeln
xSlut   = 1.6;   % Slut
xVek   = [xStart:h:xSlut]';   % x-vektor
U  = u';   % Spara startv�rdena

for  x = xVek(1:end-1)'   % I n�st sista steget ber�knas sista u
    u = u + h*funk(x, u);   % Fram�t Euler
    U = [U  ; u'];     % Lagrar u
end

plot (xVek, U(:,1))
xlabel ('x')
ylabel ('y')
