참조
Neculai Andrei, 슬롯 머신 기술을 사용한 비선형 최적화 애플리케이션, 스프링거 최적화 및 그 애플리케이션, 모델Cpa(4.5) 장비선형 방정식 시스템, 2013
메인파일 : cpa.gms
$onText
공기 중 프로판 연소.

Hiebert(1983)의 축소 변종.
Hiebert, K.L., 시스템을 해결하는 수학적 소프트웨어 평가
비선형 방정식. 수학 소프트웨어의 ACM 트랜잭션,
권. 8, 1983, pp.5-20.

Shacham, M., Brauner, N., Cutlip, M.B., (2002) 웹 기반 라이브러리
비선형 대수학을 해결하기 위한 수치 소프트웨어의 성능 테스트
방정식. 컴퓨터 및 화학 공학, vol. 2002년 26일, pp.547-554.

Meintjes, K., Morgan, A.P., (1990) 화학 평형 시스템
수치 테스트 문제. ACM 트랜스. 수학. 소프트웨어, 16, 1990, pp. 143-151.
$offText

스칼라 R /40/;

변수 x1,x2,x3,x4,x5,x6,x7,x8,x9,x10, obj;

방정식 e1,e2,e3,e4,e5,e6,e7,e8,e9,e10, eobj;

* 제약:

e1.. x1 + x4 - 3 =e= 0;
e2.. 2*x1 + x2 + x4 + x7 + x8 + x9 + 2*x10 - R =e= 0;
e3.. 2*x2 + 2*x5 + x6 + x7 - 8 =e= 0;
e4.. 2*x3 + x5 - 4*R =e= 0;
e5.. x1*x5 - 0.193*x2*x4 =e= 0;
e6.. x6*sqrt(x2) - 0.002597*sqrt(x2*x4*(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10)) =e= 0;
e7.. x7*sqrt(x4) - 0.003448*sqrt(x1*x4*(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10)) =e= 0;
e8.. x4*x8 - 1.799*x2*(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10)/100000 =e= 0;
e9.. x4*x9 - 0.0002155*x1*sqrt(x3*(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10)) =e= 0;
e10.. x10*sqr(x4) - 3.84 *sqr(x4)*(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10)/100000 =e= 0;

* 목표:

eobj..obj =e= 1;

*변수에 바인딩됨:
x1.lo = 0.000001;    x1.up = 100;
x2.lo = 0.000001;    x2.up = 100;
x3.lo = 0.000001;    x3.up = 100;
x4.lo = 0.000001;    x4.up = 100;
x5.lo = 0.000001;    x5.up = 100;
x6.lo = 0.000001;    x6.up = 100;
x7.lo = 0.000001;    x7.up = 100;
x8.lo = 0.000001;    x8.up = 100;
x9.lo = 0.000001;    x9.up = 100;
x10.lo = 0.000001;   x10.up = 100;

* 초기 지점:
x1.l = 2;
x2.l = 5;
x3.l = 40;
x4.l = 1;
x5.1 = 0;
x6.l = 0;
x7.l = 0;
x8.l = 0;
x9.l = 0;
x10.l= 5;

모델 cpa /all/;

nlp를 사용하여 obj를 최소화하는 cpa를 해결합니다.

$ifThenI x%mode%==xbook
파일 소품 /cpa1.dat/
소품을 넣어;
x1.l:12:7, x2.l:12:7, x3.l:12:7 x4.l:12:7, x5.l:12:7,
    x6.l:12:7, x7.l:12:7, x8.l:12:7 x9.l:12:7, x10.l:12:7/;
$endIf

* CPA 종료