% Skript som implementerar Eulers metode
% Input: Startkravet y(x0)=y0 og F(x,y),
% som definerar differensiallikninga ved
% y'=F(x,y)

% Startkrav og maksimal x-verdi
x0=0;
y0=1;
xMax=4;

% Differensiallikning
F=@(x,y) sqrt(y)- x^2*cos(x/2);

% Steglengda
N=input('Gi talet på steg: ');
h=(xMax-x0)/N;

% Vektor med x-verdiar
xVektor=x0:h:xMax;

% Initerar x og y
x=x0; y=y0;

% for-løkke som implementerar Eulers metode
yVektor(1)=y0;
for n=1:N
% Eulers framover-metode
  y=y+F(x,y)*h;
% Oppdaterar vektor med y-verdiar
  yVektor(n+1)=y;
% Oppdaterar x
  x=x+h;
end  

% Plottar resultatet
plot(xVektor,yVektor,'k-','linewidth',1)