A-組合投資收益和風(fēng)險問題的數(shù)學(xué)建模-數(shù)學(xué)建模課程設(shè)計

A-組合投資收益和風(fēng)險問題的數(shù)學(xué)建模-數(shù)學(xué)建模課程設(shè)計A-組合投資收益和風(fēng)險問題的數(shù)學(xué)建模-數(shù)學(xué)建模課程設(shè)計A-組合投資收益和風(fēng)險問題的數(shù)學(xué)建模-數(shù)學(xué)建模課程設(shè)計xxx公司A-組合投資收益和風(fēng)險問題的數(shù)學(xué)建模-數(shù)學(xué)建模課程設(shè)計文件編號：文件日期：修訂次數(shù)：第1.0次更改批準審核制定方案設(shè)計，管理制度安徽工程大學(xué)

數(shù)學(xué)建模

課程設(shè)計論文組合投資收益和風(fēng)險問題的數(shù)學(xué)建模姓名：何健堃學(xué)號：16班級：數(shù)學(xué)112指導(dǎo)老師：周金明成績：完成日期：2013年7月3日摘要在現(xiàn)代商業(yè)、金融投資中，投資者總是希望實現(xiàn)收益最大化，然而投資是要承擔(dān)風(fēng)險的，收益與風(fēng)險之間存在難以調(diào)和的矛盾，怎樣兼顧兩者，尋找切實可行的決策思想，是投資的收益和風(fēng)險決策的一個重要的問題。本論文主要討論并解決了在組合投資問題中的投資收益與風(fēng)險的有關(guān)問題。分別在不考慮投資項目之間的影響和考慮投資項目之間的影響以及不考慮風(fēng)險和考慮風(fēng)險的情況下，建立相應(yīng)的數(shù)學(xué)模型，來使得投資所獲得的總利潤達到最大。問題一是一典型的線性規(guī)劃問題。根據(jù)題目要求，要求第五年末的最大利潤，則建立線性規(guī)劃模型，在LINDO中編程求得第五年末的最大利潤為萬元。第一年投資項目有1，2，3，4，5，6，投資額分別為，，，，，萬元；第二年投資項目有1，2，7，投資額分別為，，萬元；第三年投資項目有1，2，8，投資額分別為，，萬元；第四年投資項目有1，2，3，4，投資額分別為，，，萬元；第五年投資項目有1，2，投資額分別為，萬元。

問題二是在問題一的基礎(chǔ)上，增加了約束條件（考慮項目間的影響）的組合投資問題。建立非線性規(guī)劃模型，在LINGO中求解得到第五年末的最大利潤為萬元。第一年投資項目有1，2，3，4，投資額分別為，，，萬元；第二年投資項目有1，2，5，投資額分別為，，

萬元；第三年投資項目有1，2，6，投資額分別為，，萬元；第四年投資項目有1，2，3，4，投資額分別為，，，萬元；第五年投資項目有1，2，投資額分別為，萬元。問題三在問題二的基礎(chǔ)上考慮投資風(fēng)險，要求風(fēng)險最小，收益最大，是一雙目標函數(shù)問題，求解過程中假設(shè)投資者能承受的最大風(fēng)險度，進而將其轉(zhuǎn)化為單目標問題。在LINGO中求解得到第五年末的最大利潤為萬元。第一年投資項目有1，2，3，4，6投資額分別為，，，，萬元；第二年投資項目有1，2，3，5，6投資額分別為，，，，萬元；第三年投資項目有1，2，3，4，6，投資額分別為，，，，萬元；第四年投資項目有1，2，3，投資額分別為，，萬元；第五年投資項目有1，2，投資額分別為，萬元。關(guān)鍵詞：線性規(guī)劃非線性規(guī)劃LINGO風(fēng)險1、問題重述由題意知，現(xiàn)有一筆20億的資金作為未來5年的投資資金，不是將其投資于一個項目上，而是有8個項目可供選擇投資，并要求：項目1、項目2每年初投資，當年年末回收本利（本金和利潤）；項目3、項目4每年初投資，要到第二年末才可回收本利；項目5、項目6每年初投資，要到第三年末才可回收本利；項目7只能在第二年年初投資，到第五年末回收本利；項目8只能在第三年年初投資，到第五年末回收本利。1.公司財務(wù)分析人員的實驗數(shù)據(jù)，求怎樣安排投資，使得第五年末的利潤最大，也即是五年累計的利潤最大？2.考慮各項目投資的相互影響的同時又考慮了其他的因素：項目5的投資額固定，為500萬，可重復(fù)投資；對投資項目1，公司管理層爭取到一筆資金捐贈，若在項目1中投資超過20000萬，則同時可獲得該筆投資金額的1%的捐贈，用于當年對各項目的投資。每年對各個項目投資的到期利潤率各不相同且受到某些項目同時投資的影響，如同時對項目3和項目4投資時的到期利潤率與項目3和項目4獨立投資不同。所以應(yīng)在模型一中加以更多的限制來求解第五年末所得最大利潤。3.本題考慮投資風(fēng)險，并用這筆資金投資若干種項目時，總體風(fēng)險可用所投資的項目中最大的一個風(fēng)險來度量，來求如何投資，使得總體風(fēng)險最小，第五年末的利潤最大。其它約束條件與二問相同。2、條件假設(shè)1.假設(shè)市場經(jīng)濟在五年內(nèi)基本穩(wěn)定。2.無交易費和投資費用等的費用開支。3.假設(shè)在五年內(nèi)產(chǎn)生的利潤可作為下一年的投資資金。4.對問題一：（1）假設(shè)不考慮各投資項目間的影響以及存在的風(fēng)險。（2）假設(shè)在五年內(nèi)各投資項目的到期利潤率不變。5.對問題二：（1）假設(shè)不考慮投資風(fēng)險。（2）假設(shè)仍然考慮一問條件，如項目3、項目4每年初投資，要到第二年末才可回收本利。6.對問題三：（1）某些項目單獨投資或共同投資的風(fēng)險損失率（標準差）每年基本相同。（2）假設(shè)仍然考慮一問條件，如項目3、項目4每年初投資，要到第二年末才可回收本利。（3）投資者能承受的最大風(fēng)險度。3、符號說明下面將會用到以下符號：----第五年末所得的利潤。----第年末回收的本利。----每年年初的可用來投資的資金。----第年初對第個項目的投資金額。----問題一中第個項目的預(yù)計到期利潤率（）。----第個項目的投資上限。-----問題二中各項目單獨投資時，第年第個項目的到期利潤率。-----問題二中某些兩個項目同時投資時，第年第個項目的到期利潤率。-----問題二中五六八項目同時投資時的到期利潤率。-----問題二中實際采用的第年第個項目的到期利潤率。-----問題二中第一項目獲得的捐贈資金。-----問題二中第年項目五可以重復(fù)投資的次數(shù)。-----問題三中的總體投資風(fēng)險。-----問題三中實際采用的風(fēng)險損失率。-----問題三中第個項目單獨投資時的風(fēng)險損失率。-----問題三中某些項目同時投資時的風(fēng)險損失率。-----問題三中投資者能夠承受的最大風(fēng)險度。4、問題分析題目背景：隨著市場經(jīng)濟的發(fā)展，越來越多的人把目光轉(zhuǎn)移到投資上。為了獲得最大的利益，組合投資的風(fēng)險與收益成為投資者必須考慮的問題。投資者在考慮投資項目的預(yù)期收益率的同時，也要考慮投資項目存在的風(fēng)險。如何使利潤最大化，風(fēng)險最小化是投資者首要考慮的問題。下面是對該問題的分析與解答。在分析解答之前給出以下定義：預(yù)期收益率：預(yù)期收益率也稱為期望收益率，是指如果事件不發(fā)生的話可以預(yù)計到的收益率。本題中到期收益率與之類似，即到期利潤率是指對某項目的一次投資中，到期回收利潤與本金的比值。風(fēng)險：風(fēng)險可以用收益的方差（或標準差）來進行衡量：方差越大，則認為風(fēng)險越大；方差越小，則認為風(fēng)險越小。在一定的假設(shè)下，有收益的方差（或標準差）來衡量風(fēng)險確實是合適的。風(fēng)險損失率和風(fēng)險一樣也可以用收益的方差（或標準差）來進行衡量。1.題一，根據(jù)題意，在不考慮各項投資之間的相互影響以及投資存在的風(fēng)險情況下，求第五年末的最大利潤，即五年內(nèi)利潤的和的最大值。則第五年末的最大利潤就等于五年中各項投資項目的金額與該項目的到期利潤率的乘積之和。還應(yīng)考慮的限制條件為：限制條件1：每年可用來投資的最大金額。第一年為20億元。限制條件2：對每項項目投資的最大金額。2.題二,根據(jù)題意,在考慮各項投資之間的相互影響且不考慮在投資中的風(fēng)險情況下，求第五年末的最大利潤，即五年內(nèi)利潤的和的最大值。首先，應(yīng)根據(jù)附表數(shù)據(jù)求得各年對各項目的單獨投資的到期的利潤率，以及各年對某些項目同時投資的到期利潤率，如同時對項目5、項目6和項目8投資的到期利潤率與它們單獨投資時的到期利潤率不同，應(yīng)分情況考慮，即它們是否同時投資。然后，對投資項目1應(yīng)單獨考慮，因為當項目一投資大于20000萬時，可以獲得1%的捐贈，這是我們必須考慮的。其次，項目5的投資額固定，為500萬，可重復(fù)投資。所以應(yīng)考慮項目5在五年內(nèi)的投資次數(shù)。限制條件1：每年可用來投資的最大金額。限制條件2：對每項項目投資的最大金額。限制條件3：同時投資時的到期利潤率不同。3.題三,根據(jù)題意考慮各投資項目的風(fēng)險情況下，求第五年末的最大利潤，即五年內(nèi)利潤的和的最大值，并且使得風(fēng)險最小。此問題為一雙目標規(guī)劃問題。首先，應(yīng)該把總的風(fēng)險表示出來，即總體風(fēng)險可用所投資的項目中最大的一個風(fēng)險來度量。風(fēng)險可以用收益的方差或標準差來衡量。其次，在考慮風(fēng)險的基礎(chǔ)上來求解如何組合投資使得風(fēng)險最小而利潤最大。限制條件1：每年可用來投資的最大金額。限制條件2：對每項項目投資的最大金額。限制條件3：同時投資時的到期利潤率不同。限制條件4：同時投資與非同時投資時的風(fēng)險損失率（標準差）不同。5、模型建立問題一：不考慮投資風(fēng)險和各投資項目間的相互影響，建立線形規(guī)劃模型求解。（1）目標函數(shù)：第五年末所得利潤，即五年內(nèi)各項目產(chǎn)生的總利潤（2）約束條件1：每年年初可用來投資的金額應(yīng)大于等于每年年初投資各項目的總額：第年年初可用來投資的資金等于第年初未投資的資金加上第年末回收的本利。第一年初可用來投資的資金第一年末回收第一年投資中一二項目的本利第二年初可用來投資的資金第二年末回收第一年三四項目和第二年一二項目的本利第三年初可用來投資的資金第三年末回收第一年五六項目、第二年三四項目和第三年一二項目的本利第四年初可用來投資的資金第四年末回收第二年五六項目、第三年三四項目和第四年一二項目的本利第五年初可用來投資的資金第五年末回收第二年七項目、第三年五六八項目、第四年三四項目、第五年一二項目的本利約束條件2：每項項目在投資期間不能超過項目最大金額的限制第一二項目每年均可以投資，可以投五次：第三四項目每年初投資，第二年末才可回收本利，可投資四次：第五六項目每年初投資，第三年末才可回收本利，可投資三次：第七項目第二年初投資，到第五年末回收本利，可投資一次：第八項目第三年初投資，到第五年末回收本利，可投資一次：所以可以建立線形規(guī)劃模型，利用LINGO軟件或LINGDO軟件求解。問題二：不考慮投資風(fēng)險時，建立非線性規(guī)劃模型求解。根據(jù)題目中表二、表三分別求得各投資項目獨立投資時的投資額及到期利潤率（見表一）以及一些投資項目同時投資時的投資額及到期利潤（見表二）。項目12345678到期利潤率20082009201020112012表一.各投資項目獨立投資時的到期利潤率項目同時投資項目3、4同時投資項目5、6同時投資項目5、6、83456568到期利潤率20082009201020112012表二.一些投資項目同時投資時到期利潤率（2）目標函數(shù)：第五年末所得利潤，即五年內(nèi)各項目產(chǎn)生的總利潤（項目單獨投資時，某些項目同時投資時）（3）約束條件1：某些項目同時投資時的到期利潤率。1）當萬元時，第年可得捐贈資金，作為第年初的投資資金，當萬元時，第年可得捐贈資金。2）當三四項目同時投資時，，此時，否則。3）同理：當五六項目同時投資時，，此時，否則。4）當?shù)谌瓿跷辶隧椖客瑫r投資時，，此時，否則。5）項目五為固定投資500萬且可以重復(fù)投資，，。約束條件2：每年年初可用來投資的金額應(yīng)大于等于每年年初投資各項目的總額：第一年初可用來投資的資金第一年末回收第一年投資中一二項目的本利第二年初可用來投資的資金第二年末回收第一年三四項目和第二年一二項目的本利第三年初可用來投資的資金第三年末回收第一年五六項目、第二年三四項目和第三年一二項目的本利第四年初可用來投資的資金第四年末回收第二年五六項目、第三年三四項目和第四年一二項目的本利第五年初可用來投資的資金第五年末回收第二年七項目、第三年五六八項目、第四年三四項目、第五年一二項目的本利約束條件2：每項項目在投資期間不能超過項目最大金額的限制第一二項目每年均可以投資，可以投五次：第三四項目每年初投資，第二年末才可回收本利，可投資四次：第五六項目每年初投資，第三年末才可回收本利，可投資三次：第七項目第二年初投資，到第五年末回收本利，可投資一次：第八項目第三年初投資，到第五年末回收本利，可投資一次：所以可以建立非線性規(guī)劃模型，利用LINGO軟件或LINDO軟件求解。問題三：考慮投資風(fēng)險的情況下，求第五年末的最大利潤，并使風(fēng)險最小。為此建立雙目標規(guī)劃，但為計算簡便，通過一定的假設(shè)條件將其轉(zhuǎn)化為單目標規(guī)劃問題。根據(jù)題目，我們用預(yù)計到期收益率的標準差來衡量風(fēng)險損失率。根據(jù)各投資項目獨立投資時的到期利潤率（表一）和一些投資項目同時投資時到期利潤率（表二），利用EXCEL中的STDEV函數(shù)，求得各項目單獨投資時的風(fēng)險損失率（表三）以及某些項目同時投資的風(fēng)險損失率（表四）：項目1項目2項目3項目4項目5項目6項目7項目8標準差表三.各項目單獨投資時的風(fēng)險損失率項目項目3、4影響項目5、6影響項目5、6、8影響項目3項目4項目5項目6項目5項目6項目8標準差表四.某些項目同時投資時的風(fēng)險損失率在這里我們定義項目風(fēng)險度=項目風(fēng)險損失率*項目投資額/投資資金總額。則可以得到目標函數(shù)：風(fēng)險最小利潤最大為了簡化模型，我們給定投資者能承受的最大風(fēng)險，則將雙目標函數(shù)轉(zhuǎn)化為單目標函數(shù)。也就是說：在滿足投資者能承受的最大風(fēng)險的情況下，如何使收益最大。目標函數(shù)：約束條件：約束條件1：與問題二中約束條件完全一樣，這里不再重復(fù)。約束條件2：總體投資風(fēng)險不能超過投資者能承受的最大風(fēng)險則項目的風(fēng)險度應(yīng)小于投資者能承受的最大風(fēng)險。即6模型求解問題一：根據(jù)上面的分析和建立的模型知：模型化解:實際在本題中可以使計算簡便。由以上分析并在LINDO中編程（源程序見附件一）得到第五年末所得利潤的最大值為萬元。各年的投資項目見下表（單位：萬元）：第一年第二年第三年第四年第五年項目一項目二項目三00--項目四00--項目五000--項目六000--項目七--------項目八--------由上表可知：第一年投資項目有1，2，3，4，5，6，投資額分別為，，，，，萬元；第二年投資項目有1，2，7，投資額分別為，，萬元；第三年投資項目有1，2，8，投資額分別為，，萬元；第四年投資項目有1，2，3，4，投資額分別為，，，萬元；第五年投資項目有1，2，投資額分別為，萬元。問題二：根據(jù)上面的分析和建立的模型知：由以上分析并在LINGO中編程（源程序見附件二）得到第五年末所得利潤的最大值為萬元。各年的投資項目見下表（單位：萬元）：第一年第二年第三年第四年第五年項目一項目二項目三00--項目四00--項目五00----項目六00----項目七--0------項目八----0----由上表可知：第一年投資項目有1，2，3，4，投資額分別為，，，萬元；第二年投資項目有1，2，5，投資額分別為，，萬元；第三年投資項目有1，2，6，投資額分別為，，萬元；第四年投資項目有1，2，3，4，投資額分別為，，，萬元；第五年投資項目有1，2，投資額分別為，萬元。問題三：根據(jù)上面的分析和建立的模型知：由以上分析并在LINGO中編程（源程序見附件三）得到第五年末所得利潤的最大值為萬元。各年的投資項目見下表（單位：萬元）：第一年第二年第三年第四年第五年項目一項目二項目三--項目四00--項目五00----項目六----項目七--0------項目八----0----由上表可知：第一年投資項目有1，2，3，4，6投資額分別為，，，，萬元；第二年投資項目有1，2，3，5，6投資額分別為，，，，萬元；第三年投資項目有1，2，3，4，6，投資額分別為，，，，萬元；第四年投資項目有1，2，3，投資額分別為，，萬元；第五年投資項目有1，2，投資額分別為，萬元。7、模型檢驗與評價本文均采用規(guī)劃模型實現(xiàn)了組合投資的收益和風(fēng)險問題的優(yōu)化投資方案。問題一：問題一是一典型的線性規(guī)劃問題，在不考慮投資項目間的影響和投資風(fēng)險時，不難得到目標函數(shù)和約束條件，在LINDO中很快得到全局最優(yōu)解。問題二：問題二在問題一的基礎(chǔ)上，把投資項目間的影響考慮了進來，增加了約束條件。由于各項目各年的收益率不同，增加了求解模型的難度。且將模型變得非線性化。在本模型中用LINGO的if選擇語句實現(xiàn)其在組合投資中對個項目是否同時投資的選擇，在一定程度上實現(xiàn)了非線性模型的建立,并借助第一問的約束條件，在LINGO中編程，求得局部最優(yōu)解。問題三：問題三在以上兩問的基礎(chǔ)上，又把投資風(fēng)險考慮進來，要求在風(fēng)險最小的情況下建立模型求得投資方案以獲得最大利潤，為雙目標規(guī)劃模型，但在求解上有一定的難度，于是在此轉(zhuǎn)化為單目標規(guī)劃模型。即假設(shè)投資者能承受的最大風(fēng)險度，這樣將雙目標問題轉(zhuǎn)化為與二問類似的問題，并且投資越分散說明投資風(fēng)險越小。進而得到各年的投資計劃，從投資計劃中看出，該投資計劃相對比一二問的投資項目分散，這說明該模型的假設(shè)還是比較合理的。當改變假設(shè)的最大風(fēng)險度時，第五年末的最大利潤也發(fā)生相應(yīng)的變化。時，最大利潤；時，最大利潤；時，最大利潤；由上面的數(shù)據(jù)分析可知，投資者能承受的最大風(fēng)險度越大，所得的利潤就越大。但這還是有一定的局限性，對此還可以考慮以下方案：（1）在投資者要求的最低利潤條件下，來求風(fēng)險的最小值。即給定一個最小利潤，在以此為約束的條件下來求風(fēng)險的最小值，這樣也能在一定程度上將雙目標函數(shù)優(yōu)化為單目標函數(shù)。（2）采用加權(quán)法，對風(fēng)險和投資利潤加于權(quán)重因子構(gòu)造出目標函數(shù)，再以一定的約束條件以及權(quán)重因子的約束為條件，用LINGO軟件編程求解相應(yīng)的投資方案。8、參考文獻[1]韓中庚。數(shù)學(xué)建模方法及其應(yīng)用。北京：高等教育出版社，2006.[2]謝金星，薛毅。優(yōu)化建模與LINDO/LINGO軟件。北京：清華大學(xué)出版社，2006.9、附錄附源程序一：（LINDO中求解）max+++++++++++++++++++++++++STx11<60000x21<60000x31<60000x41<60000x51<60000x12<30000x22<30000x32<30000x42<30000x52<30000x23+x13<40000x33+x23<40000x43+x33<40000x24+x14<30000x34+x24<30000x44+x34<30000x25+x15+x35<30000x26+x16+x36<20000x27<40000x38<30000x11+x12+x13+x14+x15+x16<200000x21+x22+x23+x24+x25+x26+附源程序二：（LINGO中求解）max=*x11+*x12+b13*x13+b14*x14+b15*x15+b16*x16+*x21+*x22+b23*x23+b24*x24+b25*x25+b26*x26+*x27+*x31+*x32+b33*x33+b34*x34+b35*x35+b36*x36+b38*x38+*x41+*x42+b43*x43+b44*x44+*x51+*x52;x15=500*n1;x25=500*n2;x35=500*n3;b13=@if(x13*x14#EQ#0,,;b14=@if(x13*x14#EQ#0,,;b15=@if(x15*x16#EQ#0,,;b16=@if(x15*x16#EQ#0,,;b23=@if(x23*x24#EQ#0,,;b24=@if(x23*x24#EQ#0,,;b25=@if(x25*x26#EQ#0,,;b26=@if(x25*x26#EQ#0,,;b33=@if(x33*x34#EQ#0,,;b34=@if(x33*x34#EQ#0,,;b35=@if(x35*x36#EQ#0,,@if(x38#EQ#0,,);b36=@if(x35*x36#EQ#0,,@if(x38#EQ#0,,);b38=@if(x35*x36*x38#NE#0,,;b43=@if(x43*x44#EQ#0,,;b44=@if(x43*x44#EQ#0,,;x11<60000;x21<60000;x31<60000;x41<60000;x51<60000;x12<60000;x22<60000;x32<60000;x42<60000;x52<60000;x23+x13<35000;x33+x23<35000;x43+x33<35000;x24+x14<30000;x34+x24<30000;x44+x34<30000;x25+x15+x35<30000;x26+x16+x36<40000;x27<30000;x38<30000;g11=@if(x11#LT#20000,0,x11*;g21=@if(x21#LT#20000,0,x21*;g31=@if(x31#LT#20000,0,x31*;g41=@if(x41#LT#20000,0,x41*;g51=@if(x51#LT#20000,0,x51*;x11+x12+x13+x14+x15+x16-g11<200000;x21+x22+x23+x24+x25+x26+**x12+x13+x14+x15+x16-g21<200000;x31+x32+x33+x34+x35+x36+****x12-b13*x13-b14*x14+x23+x24+x15+x16+x25+x26+x27-g31<200000;x41+x42+x43+******x12-b23*x23-b24*x24-b13*x13-b14*x14-b15*x15-b16*x16+x33+x34+x25+x26+x35+x36+x27+x38-g41<200000;x51+********x12-x33*b33-x34*b34-b23*x23-b24*x24-b13*x13-b14*x14-x25*b25-x26*b26-x15*b15-x16*b16+x43+x44+x35+x36+x27+x38-g51<200000;@gin(n1);@gin(n2);@gin(n3);@free(b13);@free(b23);@free(b33);@free(b43);@free(b14);@free(b24);@free(b34);@free(b44);@free(b15);@free(b25);@free(b35);@free(b16);@free(b26);@free(b36);@free(b38);@free(b27);end附源程序三：（LINGO中求解）max=*x11+*x12+b13*x13+b14*x14+b15*x15+b16*x16+*x21+*x22+b23*x23+b24*x24+b25*x25+b26*x26+*x27+*x31+*x32+b33*x33+b34*x34+b35*x35+b36*x36+b38*x38+*x41+*x42+b43*x43+b44*x44+*x51+*x52;y0=200000;x15=500*n1;x25=500*n2;x35=500*n3;b13=@if(x13*x14#EQ#0,,;b14=@if(x13*x14#EQ#0,,;b15=@if(x15*x16#EQ#0,,;b16=@if(x15*x16#EQ#0,,;b23=@if(x23*x24#EQ#0,,;b24=@if(x23*x24#EQ#0,,;b25=@if(x25*x26#EQ#0,,;b26=@if(x25*x26#EQ#0,,;b33=@if(x33*x34#EQ#0,,;b34=@if(x33*x34#EQ#0,,;b35=@if(x35*x36#EQ#0,,@if(x38#EQ#0,,);b36=@if(x35*x36#EQ#0,,@if(x38#EQ#0,,);b38=@if(x35*x36*x38#NE#0,,;b43=@if(x43*x44#EQ#0,,;b44=@if(x43*x44#EQ#0,,;o13=@if(x13*x14#EQ#0,,;o23=@if(x23*x24#EQ#0,,;o33=@if(x33*x34#EQ#0,,;o43=@if(x43*x44#EQ#0,,;o14=@if(x13*x14#EQ#0,,;o24=@if(x23*x24#EQ#0,,;o34=@if(x33*x34#EQ#0,,;o44=@if(x43*x44#EQ#0,,;o15=@if(x15*x16#EQ#0,,;o25=@if(x25*x26#EQ#0,,;o35=@if(x35*x36#EQ#0,,@if(x38#EQ#0,,);o16=@if(x15*x16#EQ#0,,;o26=@if(x25*x26#EQ#0,,;o36=@if(x35*x36#EQ#0,,@if(x38#EQ#0,,);o38=@if(x35*x36*x38#NE#0,,;o11=;o21=;o31=;o41=;o51=;o12=;o22=;o32=;o42=;o52=;o27=;a=;o11*x11/(y0+g11)<a;o21*x21/(y0+*x11+*x12-x13-x14-x15-x16+g21)<a;o31*x31/(y0+*x21+*x22+*x11+*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o41*x41/(y0+*x31+*x32+*x21+*x22+*x11+*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36-x27-x38+g41)<a;o51*x51/(y0+*x41+*x42+*x31+*x32+*x21+*x22+*x11+*x12+x33*b33+x34*b34+b23*x23+b24*x24+b13*x13+b14*x14+x25*b25+x26*b26+x15*b15+x16*b16-x43-x44-x35-x36-x27-x38+g51)<a;o12*x12/(y0+g11)<a;o22*x22/(y0+*x11+*x12-x13-x14-x15-x16+g21)<a;o32*x32/(y0+*x21+*x22+*x11+*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o42*x42/(y0+*x31+*x32+*x21+*x22+*x11+*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36-x27-x38+g41)<a;o52*x52/(y0+*x41+*x42+*x31+*x32+*x21+*x22+*x11+*x12+x33*b33+x34*b34+b23*x23+b24*x24+b13*x13+b14*x14+x25*b25+x26*b26+x15*b15+x16*b16-x43-x44-x35-x36-x27-x38+g51)<a;o13*x13/(y0+g11)<a;o23*x23/(y0+*x11+*x12-x13-x14-x15-x16+g21)<a;o33*x33/(y0+*x21+*x22+*x11+*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o43*x43/(y0+*x31+*x32+*x21+*x22+*x11+*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36-x27-x38+g41)<a;o14*x14/(y0+g11)<a;o24*x24/(y0+*x11+*x12-x13-x14-x15-x16+g21)<a;o34*x34/(y0+*x21+*x22+*x11+*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o44*x44/(y0+*x31+*x32+*x21+*x22+*x11+*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36-x27-x38+g41)<a;o15*x15/(y0+g11)<a;o25*x25/(y0+*x11+*