
%Differensiallinkningen er y' = F(x,y)
F =@(x,y)   -x/y ;
%sqrt(1-y^2) ;
%y*(1-y);
%- sqrt(abs(y));
%-x/y ;  

a=0;  %startverdi for x
b=1;   %stoppverdi for x
N=100; %antall delintervaller
d= (b-a)/N ;  %steglengde 

y0 =0.5 ; % Initialbetingelsen er y(a) = y0


X=a:d:b ; %radvektor med N+1 verdier fra a til b med steglengde d
Y=zeros(1,N+1); %Radvektor med 0-elementer
Y(1) =y0; %setter Y(a) = y0

for n= 1:N
    Y(n+1) = Y(n) + d * F( X(n) , Y(n) ) ;
    %Eulers metode implementeres 
end
   plot(X, Y)  %plotter resultatet
   
   hold on
 [z,w]= ode45(F,[a b],y0); %benytter en av matlabs innebygde metoder 
 %for å estimere løsnigner til differensiallikninger
 plot(z,w, 'r') % Plotter resultatet med rød kurve.
 hold off