* Test problem 10.1.4 in the web page
* Test problem 10.2.5 in test problem collection
* Test problem from Cottle et al 1992
* The third solution is obtained with the first solution
* values used as the starting point

SETS J /1*4/;

alias(I,J);

PARAMETER M(I,J) matrix/

1.1  0.
1.2  0.
1.3  10.
1.4  20.
2.1  0.
2.2  0.
2.3  30.
2.4  15.
3.1  10.
3.2  20.
3.3  0.
3.4  0.
4.1  30.
4.2  15.
4.3  0.
4.4  0.
/;

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 COTTLE /eq.x/;
x.l('1') = 0.1;
x.l('2') = 0;
x.l('3') = 0.1;
x.l('4') = 0;

SOLVE COTTLE using mcp;