// Données du problème
Omega   = [0,1]; // domaine
f       = -10;   // chargement
psi     = -0.6;  // contrainte
// Paramètres de l'algorithme numérique
N       = 16;    // nombre de noeuds internes
lambda = 10;    // paramètre pour l'algorithme d'Uzawa
epsilon = 1e-6;  // tolérance pour l'algorithme d'Uzawa
kmax    = 1500;  // nombre d'itérations maximum
// Algorithme d'Uzawa
qh     = ones(N+2, 1);  // multiplicateurs de Lagrange
uh_old = zeros(N+2, 1);
for k=0:kmax
  [uh, xDOFs] = ef(N, Omega, XXXX, 1); // minimisation libre à qh fixé
  qh = max(0, XXXX);                   // mise à jour de qh
  if(norm(uh-uh_old)<=epsilon) then
    break;                             // l'algorithme a convergé
  end
  uh_old = uh;
end