function [r,center] = meb_not3(P,acc)

% Input: P is a matrix whose columns are the coordinates of each of
% the points. P is a (d) by n matrix.

[d,n] = size(P);
c = [];
S = -sparse(eye(d+1));
At = [];
b = sparse([-1;zeros(d,1)]);
K.q = [];
blk{1,1} = 'q';
blk{1,2} = [];

for i = 1:n
    At = sparse([At;S]);
    c = sparse([c;0; -P(:,i)]);
    blk{1,2} = [blk{1,2} d+1];
end;
   OPTIONS.vers          = 1;
   OPTIONS.gam           = 0;
   OPTIONS.predcorr      = 1;
   OPTIONS.expon         = [1 1];
   OPTIONS.gaptol        = acc;
   OPTIONS.inftol        = 1e-8;
   OPTIONS.steptol       = 1e-6;
   OPTIONS.maxit         = 50;
   OPTIONS.printyes      = 1;
   OPTIONS.scale_data    = 0;
   OPTIONS.randnstate    = 0;
   OPTIONS.spdensity     = 0.5;
   OPTIONS.rmdepconstr   = 0;
A{1} = At;
C{1} = c;
[obj,X,y,Z] = sqlp(blk,A,C,b,[],[],[],OPTIONS);

r = -obj(1);
center = y(2:d+1);