  a_i^Tx &lt;= b_i   i=1,...,m,</span>
<span class="comment">%           Fx = g,</span>
<span class="comment">% is the solution of the unconstrained minimization problem</span>
<span class="comment">%           minimize    -sum_{i=1}^m log(b_i-a_i^Tx).</span>

<span class="comment">% Input data</span>
randn(<span class="string">'state'</span>, 0);
rand(<span class="string">'state'</span>, 0);
n = 10;
m = 50;
p = 5;
tmp = randn(n,1);
A = randn(m,n);
b = A*tmp + 10*rand(m,1);
F = randn(p,n);
g = F*tmp;

<span