set NODI;
param sorgente symbolic in NODI;
param destinazione symbolic in NODI, <> sorgente;
set ARCHI within {NODI,NODI}; 

param cap {ARCHI} >= 0;
var x{(i,j) in ARCHI} >= 0, <= cap[i,j];

subject to equilibrio {k in NODI diff {sorgente,destinazione}}:
   sum {(i,k) in ARCHI} x[i,k] = sum {(k,j) in ARCHI} x[k,j];

maximize flusso_totale:  sum {(j,destinazione) in ARCHI} x[j,destinazione];