var eps >= 1.0e-8;

param n := 2;		# dimension
param m := 3;		# number of points

param a {1..m, 1..n};

var x {1..n};
var t {1..m};

minimize sumLone: sum {i in 1..m} t[i];

subject to t_def {i in 1..m}: 
    eps + 
    sqrt( sum {j in 1..n} (eps^2 + (x[j] - a[i,j])^2) ) <= t[i];

let eps := 1;

let a[1,1] := 0;  let a[1,2] := 0;
let a[2,1] := 4;  let a[2,2] := 0;
let a[3,1] := 2;  let a[3,2] := 4;

let x[1] := 1; let x[2] := 0;
display {i in 1..m}: sqrt( sum {j in 1..n} (x[j] - a[i,j])^2 );

solve;

display x;

let eps := 0;
display sumLone;