function KriegDF

    globals;
    KriegD;

% Berechne Gradient von F
    D3 = zeros(N+1);
    D3(1,1:2) = D(1,1:2).^3; % die negative Wurzel^3 wird wieder negativ
    D3(N+1,N:(N+1)) = D(N+1,N:(N+1)).^3;
    for i = 2:N
        D3(i,i+1) = D(i,i+1).^3;
        D3(i,i-1) = D(i,i-1).^3;
        D3(i,i) = -(D3(i,i-1) + D3(i,i+1));
    end
    DF = eye(N+1) + mue*ep^2/h^2*D3;

% ohne Schleifen:
%   ux = (u(2:(N+1))-u(1:N))/h;
%   D3 = 1./sqrt(ux.^2 + ep^2);
%   D3 = D3.^3;
%   DF = diag([D3;0]) ...
%      + diag([0;D3]) ...
%      - diag(D3,-1)  ...
%      - diag(D3,+1);
%   DF = eye(N+1) + (mue*ep^2/h^2)*DF;

end