function KriegJh

    globals;

    ux = (u(2:(N+1))-u(1:N))/h;
    Jh = (h/2)*sum((u-v).^2) + (mue*h)*sum(sqrt(ux.^2 + ep^2));

end
