function divI = Divergence(II); 

% function divI = Divergence(II); 
% 
% Calcule la divergence d'un champ de vecteurs 

[N1,N2,N3] = size(II); 
T1 = [2:N1 1]; 
T2 = [2:N2 1]; 
divI = (II(:,T2,1)-II(:,:,1)) + (II(T1,:,2)-II(:,:,2)); 
end

