% Lösung der zusätzlichen AWA, die zur Durchführung 
% des Newtonverfahrens benötigt wird

function[dG] = dGG(a,b,N,us)

u0 = [0;1];
[x,u] = euler_system_dg(a,b,u0,N,us);
dG = - u(N+1);