% Generate feasible but nonsolvable LCP
n= input('no of variables: ');
q=rand(n-1,1)*10;
tmp=min(q);
q(n)=tmp/(1+rand(1,1));
q=-q;
M=zeros(n,n);
M(1:n-2,1:n-1)=rand(n-2,n-1);
M(1:n-2,n)=-rand(n-2,1);
M(n-1,n)=1.0;
for i=1:n-1
	M(n,i)=1.0;
end;
fprintf('data','%d\n',n);
fprintf('data','%f %f %f %f %f \n',M');
fprintf('data','\n');
fprintf('data','%f %f %f %f %f \n',q);
fprintf('data','\n');
save dat n;
save dat M;
save dat q;


