% Explizites Eulerverfahren für die im Newtonverfahren
% zusätzlich zu lösende AWA

function[x,u] = euler_system_dG(a,b,u0,N,us)

h = (b-a)/N;
v(:,1) = u0;   

for i = 1:N
    x(i) = a+(i-1)*h;
    v(:,i+1) = v(:,i) + h*fw_dg(x(i),v(:,i),us(i));    
end    

u = v(1,:);
x(N+1) = b;