Warning in get_ep_min: 2 solutions
Error using  ^ 
Inputs must be a scalar and a square matrix.
To compute elementwise POWER, use POWER (.^) instead.

Error in get_lb (line 71)
  n = 1000 + round(1000 * max(0, k - 1)); xb = g/ (g + eb); xb3 = xb ^ (1/3);

Error in get_ue0 (line 49)
    [lb, info] = get_lb(p, eb);