function[x,u] = schiess(a,b,ua,ub,N)

% Startwert der Anfangssteigung
s = 1;

% Newtoniteration zur Nullstellenbestimmung
for  n = 1:10
    [G,us] = GG(s,a,b,ua,ub,N);
    [dG] = dGG(a,b,N,us);
    s = s - G/dG;
end

% Lösung der AWA mit berechneter Anfangssteigung
u0 = [ua;s];
[x,u] = euler_system(a,b,u0,N);