* Test problem 10.1.1 in the web page
* Test problem 10.2.1 in test problem collection
* Test problem from MCPLIB Dirkse and Ferris 1995

SETS J /1*16/;

alias(I,J);

PARAMETER M(I,J) matrix;

M(I,J) = 2$(ord(J) gt ord(I))+ 1$(ord(J) eq ord(I));

PARAMETER q(J) vector;

q(J) = -1;

VARIABLES x(J);
EQUATIONS eq(I);

eq(I).. sum(J, M(I,J)*x(J)) + q(i) =e= 0.0;
MODEL DIRKSE /eq.x/;
x.lo(J) = 0.0;

SOLVE DIRKSE using mcp;