clear;
n = 200; m = 100;

A = zeros(1,n); b = 0;

Aeq = rand(m,n);
xfeas = 1.0+ 2*rand(n,1);
beq = Aeq * xfeas;

lb = zeros(n,1);
% ub=inf*ones(n,1);
ub=[];

f = -rand(n,1);

options = optimset('Display','iter');
x0=rand(n,1);
tic;
[x,fval,exitflag, output, lambda] = linprog(f,A,b,Aeq,beq,lb,ub,x0,options);
toc

% x
% lambda.eqlin
% lambda.lower