function [A,c,b,M,B,x0,eps]=ex5(x) s=2*x; A(1,1:x-1)=2; A(1,x)=1; A(2:x+1,:)=eye(x); A(:,x+1:s+1)=eye(x+1); A=A(:,[1:x x+2:s+1]) c=zeros(2*x,1); c(x)=-1 b=[x]; b(2:x+1)=ones(x,1) b=b'; M=[1:x] B=[x:s] x0=zeros(x-1,1); x0(x)=x; x0(x+1:s)=ones(x,1); x0=x0' eps=1e-3; return;