function E = FonctionEnergiePhi(VV,DeltaI,I,lambda); 
% 
% function E = FonctionEnergiePhi(VV, DeltaI,I,lambda); 
%
% Calcul la valeur de la fonctionnelle Phi en un point 
% VV (fonctionnelle a minimiser lors de la recherche du flot optique optimal 
% avec la methode du gradient conjugue et un parametre de relaxation
% lambda). 

[NablaI1,NablaI2] = Gradient(I); 
NablaI=cat(3,NablaI1,NablaI2); 
UU = DeltaI + sum(VV.*NablaI,3); 
[nabla11,nabla12] = Gradient(VV(:,:,1)); 
[nabla21,nabla22] = Gradient(VV(:,:,2));
EnergieGradVV = (norm(nabla11))^2 + (norm(nabla12))^2 + (norm(nabla21))^2 + (norm(nabla22))^2; 
E = (norm(UU))^2 + lambda*EnergieGradVV; 
end

