fid = fopen('matlab_opt', 'r');

m = fscanf(fid, '%d', 1);
n = fscanf(fid, '%d', 1);
c = fscanf(fid, '%d', 1);

C = fscanf(fid, '%f', [n m])';
d = fscanf(fid, '%f', m);
Aeq = fscanf(fid, '%f', [n c])';
beq = fscanf(fid, '%f', c);
lb = fscanf(fid, '%f', n);
ub = fscanf(fid, '%f', n);

disp('With bounds');
x = lsqlin(C, d, [], [], Aeq, beq, lb, ub);

x2 = lsqlin(C, d, [], [], Aeq, beq);
