$ontext We model a QVI example Scrim12 taken from the following paper: QVILIB: A library of quasi-variational inequality test problems, Francisco Facchinei, Christian Kanzow, and Simone Sagratella, 2012 Our framework finds a solution. Contributor: Youngdae Kim (06.19.2017) $offtext $set N 1200 set t / 1*%N% / i / 1*2 / j / 1*4 / ; alias(j,jj); table A(j,j) 1 2 3 4 1 4 2 0 0 2 2 10 0 4 3 0 0 2 1 4 0 2 1 5; parameter b(j) / 1 40, 2 30, 3 40, 4 30 /; table C(i,j) 1 2 3 4 1 1 1 0 0 2 0 0 1 1; table E(i,j) 1 2 3 4 1 0 0 0 0 2 0 0 0 0.5; E('1','1') = 2/3; parameter d1(i) / 1 1, 2 3 /; parameter d2(i) / 1 10, 2 4 /; positive variables y(t,j), x(t,j); equations F(t,j), M(t,i); * We have to add eps*y(t,j) + eps*x(t,j) to make sure that all variables * appear in the equation. Otherwise, JAMS will complain that the sizes of * y and x do not match. F(t,j).. sum(jj, A(j,jj)*x(t,jj)) + b(j) + eps*y(t,j) + eps*x(t,j) =N= 0; M(t,i).. sum(j, C(i,j)*y(t,j)) - sum(j, E(i,j)*x(t,j)) - d1(i) - d2(i)*(t.val - 1)/(%N% - 1) =E= 0; model scrim12 / F, M /; file empinfo / '%emp.info%' /; putclose empinfo 'qvi F y x M'; solve scrim12 using emp; display y.l, x.l;