$ontext We model a QVI example BiLin1B 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 y*=(0.001,0.008,-0.008,-0.055,-0.130). Contributor: Youngdae Kim (06.18.2017) $offtext set i / 1*5 / k / 1*3 /; alias(i,j); table A(i,j) 1 2 3 4 5 1 19.8699 0.5369 2.9482 0.3358 7.1239 2 4.1819 16.3484 -5.2030 5.4332 2.7143 3 -5.6554 0.9422 19.0981 7.1556 -7.3810 4 -1.8770 0.1918 -5.3596 18.3565 -7.8847 5 -6.0303 -3.6171 -1.4658 4.6238 15.4085; parameter b(i) / 1*5 10 /; parameter lo(i) / 1 -0.1202, 2 -1.7418, 3 -2.7064, 4 -2.0502, 5 -4.4616 /; parameter c(k) / 1 0.3070, 2 1.1186, 3 2.6149 /; table leftQ1(i,j) 1 2 3 4 5 1 1.9073 0.2403 0.2352 -0.4903 -0.2651 2 0.2403 1.1319 1.2087 -0.3268 0.2540 3 0.2352 1.2087 1.6862 0.2941 0.6732 4 -0.4903 -0.3268 0.2941 1.8258 0.1363 5 -0.2651 0.2540 0.6732 0.1363 1.5527; table leftQ2(i,j) 1 2 3 4 5 1 2.7307 0.5988 1.5728 1.4072 -0.3082 2 0.5988 2.2435 0.7546 1.3632 1.5852 3 1.5728 0.7546 2.3809 1.2625 1.0403 4 1.4072 1.3632 1.2625 1.7612 0.3071 5 -0.3082 1.5852 1.0403 0.3071 2.6305; table leftQ3(i,j) 1 2 3 4 5 1 2.5189 2.1947 1.7697 2.2753 1.9885 2 2.1947 3.8143 1.3839 1.5636 1.8451 3 1.7697 1.3839 3.3655 1.6441 1.9946 4 2.2753 1.5636 1.6441 3.6885 2.3272 5 1.9885 1.8451 1.9946 2.3272 2.2883; table M(i,j) 1 2 3 4 5 1 1 0 0 0 0 2 1 1 0 0 0 3 1 1 1 0 0 4 1 1 1 1 0 5 1 1 1 1 1; parameter alpha / 10 /; parameter Q(k,i,j); Q('1',i,j) = leftQ1(i,j) + alpha*M(i,j); Q('2',i,j) = leftQ2(i,j) + alpha*M(i,j); Q('3',i,j) = leftQ3(i,j) + alpha*M(i,j); variable y(j), x(j); equations F(i), g(k); F(i).. sum(j, A(i,j)*y(j)) + b(i) =N= 0; g(k).. sum(i, x(i)*sum(j, Q(k,i,j)*y(j))) - c(k) =L= 0; model bilin1a / F, g /; file empinfo / '%emp.info%' /; putclose empinfo 'qvi F y x g'; y.lo(j) = lo(j); y.up(j)= -lo(j); solve bilin1a using emp; display y.l, x.l;