참조

Neculai Andrei, 슬롯 사이트 기술을 사용한 비선형 최적화 애플리케이션, 스프링거 최적화 및 그 애플리케이션, 모델CpaR(4.4) 장비선형 방정식 시스템, 2013

메인파일 : cpar.gms

$onText
공기 중 프로판 연소, 제형 감소. 미지수는 다음을 나타냅니다.
프로판 1몰당 생성된 생성물의 몰수.

환원된 제제는 Meintjes, K. 및 Morgan, A.P.,
수치 테스트 문제로서의 화학 평형 시스템. ACM 트랜스. 수학.
소프트웨어, 16, 1990, pp. 143-151.

다음을 참조하세요:
Neculai Andrei, "모델, 테스트 문제 및 응용 프로그램
수학적 프로그래밍". 기술 언론, 부쿠레슈티, 2003.
애플리케이션 U86, 페이지 67. 애플리케이션 A23, 페이지 376.
$offText

스칼라 R5 /0.193/
        R6 /0.4106217541e-3/
        R7 /0.5451766686e-3/
        R8 /0.44975e-6/
        R9 /0.3407354178e-4/
        R10 /0.9615e-6/
        R /10/;

변수 x1,x2,x3,x4,x5, obj;
방정식 e1,e2,e3,e4,e5, e;

e1.. x1*x2+x1-3*x5 =e= 0;
e2.. 2*x1*x2+x1+2*R10*파워(x2,2)+x2*파워(x3,2)+R7*x2*x3+
     R9*x2*x4+R8*x2-R*x5 =e= 0;
e3.. 2*x2*파워(x3,2)+R7*x2*x3+2*R5*파워(x3,2)+R6*x3-8*x5 =e= 0;
e4.. R9*x2*x4+2*POWER(x4,2)-4*R*x5 =e= 0;
e5.. x1*x2+x1+R10*파워(x2,2)+x2*파워(x3,2)+R7*x2*x3+R9*x2*x4+R8*x2+
     R5*파워(x3,2)+R6*x3+파워(x4,2)-1 =e= 0;

예..obj =e= 1;

* 변수의 한계
x1.lo = 0.0001;
x1.up = 100;
x2.lo = 0.0001;
x2.up = 100;
x3.lo = 0.0001;
x3.up = 100;
x4.lo = 0.0001;
x4.up = 100;
x5.lo = 0.0001;
x5.up = 100;

* 초기점
x1.l = 10;
x2.l = 10;
x3.l = 0.05;
x4.1 = 50.5;
x5.l = 0.05;

모델 cpaR /all/;

obj를 최소화하는 nlp를 사용하여 cpaR을 해결합니다.
* cpaR 종료