gistfile1.py

M = Model("minimal sphere enclosing a set of points dual")

print('Declaring the variables...')
y= M.variable('y',NDSet(k,n+1),Domain.inQCone(k,n+1))

c=[0. for i in range(n+1)]
c[0]=1.

print('Defining the constraints...')
M.constraint('equalities', Expr.mul(DenseMatrix(1,k,1.0), y), Domain.equalsTo(c) )

print('Defining the objective function...')
M.objective('dual',ObjectiveSense.Maximize, Expr.sum(Expr.mulDiag(DenseMatrix(p.tolist()), y.slice([0,1],[k,n+1]).transpose())))

M.solve()