function Inskjutning(N) %N ar antal delintervall
%Intervallhalvering för lösning av F(s)=0:
TOL=1e-8;   %Tolerans för intervallhalveringsmetoden

a=0.5; b=1;     %Andpunkter pa intervallet
while b-a>TOL
    m=(a+b)/2;  %Mittpunkten av intervallet
    if F(a)*F(m)>0  %På samma sida om x-axeln
        a=m;
    else
        b=m;
    end
end
sfinal=(a+b)/2; %Lösningen av F(s)=0 med given tolerans.
%Plottning av lösningen
u=[1; sfinal]; h=1/N;
yvec=1;
for i=1:N
    u=u+h*[u(2); sqrt((i*h)^2+u(1)^2+u(2)^2)];
    yvec=[yvec;u(1)];
end
plot(1:h:2,yvec)

    function utdiff=F(s)
        %Eulers metod för lösning av differentialekvationen
        h=1/N;  %Steglängd
        u=[1; s];   %Startposition vid x=1.
        for i=1:N
            u=u+h*[u(2); sqrt((i*h)^2+u(1)^2+u(2)^2)];
        end
        utdiff=u(1)-3;
    end
end