using Convex, SCS n = 25; m = 15; A = randn(m, n); b = rand(m, 1); x = Variable(n); problem = maximize(entropy(x), sum(x) == 1, A * x <= b) solve!(problem, SCSSolver(verbose=0)) println(problem.optval) println(x.value)