數(shù)學(xué)模型課程設(shè)計(jì)

1、東北大學(xué)秦皇島分校數(shù)學(xué)建模課程設(shè)計(jì)報(bào)告企業(yè)人力計(jì)劃系 別信息與計(jì)算科學(xué)專 業(yè)學(xué) 號7080130姓 名陶韜指導(dǎo)教師張尚國 郭靜梅成 績教師評語：指導(dǎo)教師簽字： 2011年7月1日1.1 課題的背景“企業(yè)人力計(jì)劃”數(shù)學(xué)模型是為實(shí)現(xiàn)企業(yè)人力資源的合理安排，企業(yè)低成本投入的計(jì)劃模型，一個(gè)合理的人力計(jì)劃安排不僅可以保證企業(yè)的正常運(yùn)行，而且可以減少企業(yè)的額外投資，這個(gè)模型就是聯(lián)系企業(yè)經(jīng)營的實(shí)際問題和我們所學(xué)的數(shù)學(xué)知識，就企業(yè)如何合理安排人員保證每年的人力需求，怎樣才能最大限度的減少額外投資的問題，作出相應(yīng)的解答和處理。問題一：根據(jù)企業(yè)現(xiàn)有人數(shù)和未來三年的人力需求，以及企業(yè)的招工，培訓(xùn)，解聘，超員招工等情

2、況建立動態(tài)整數(shù)規(guī)劃模型，企業(yè)的人力資源在現(xiàn)有人力的基礎(chǔ)上又有人力流入和人力流出兩個(gè)方面，流入包括企業(yè)招工，超員雇傭，半日工；流出包括企業(yè)解聘，和人員的離職。同時(shí)為了滿足企業(yè)對各個(gè)職工類型（熟練，半熟練，不熟練）的需求，企業(yè)內(nèi)部通過培訓(xùn)來實(shí)現(xiàn)自身職業(yè)技能的提升和人力資源的優(yōu)化利用。我們的模型實(shí)際就是“i類職工j-1年的現(xiàn)有人數(shù)+i類職工的j年的流入-i類職工j年的流出=j年企業(yè)對i類職工的需求”（j=1,2,3；j=0是指企業(yè)開始的職工人數(shù)）在綜合考慮流入，流出所包含各個(gè)因素的約束條件，就可以求出企業(yè)為實(shí)現(xiàn)目標(biāo)所能解雇的最少人數(shù)。問題二，在已經(jīng)建立好的動態(tài)整數(shù)規(guī)劃模型下，問題一是取如何能讓企業(yè)在

3、這三年解雇的人數(shù)最少為主要矛盾，并求取可能的最小值，而問題二實(shí)際上就是將問題的主要矛盾改為如何能讓企業(yè)在這三年中支出費(fèi)用最少，在已有模型下求解。通過一二問我們可以看出這個(gè)模型可以就企業(yè)人力計(jì)劃為主導(dǎo)提出多個(gè)問題的解決方案，其關(guān)鍵點(diǎn)主要是在取那個(gè)問題為主要矛盾，在保證需求的基礎(chǔ)上實(shí)現(xiàn)最優(yōu)解。1.2問題的提出人力計(jì)劃某公司正經(jīng)歷一系列變化，這要影響到它在未來幾年中的人力需求。由于裝備了新機(jī)器，對不熟練工人的需求相對減少，對熟練和半熟練工人的需求相對增加；同時(shí)，預(yù)期下一年度的貿(mào)易量將下降，從而減少對各類人力的需求?，F(xiàn)有人數(shù)及對未來三年人力需求的估計(jì)數(shù)見下表:(表)1分類不熟練半熟練熟練現(xiàn)有人數(shù)第一年

4、需求第二年需求第三年需求20001000500015001400200025001000100015002000為此，公司希望為未來三年確定（1）招工，（2）人員再培訓(xùn)，（3）解雇和超員費(fèi)用，（4）設(shè)半日工的計(jì)劃方案。因工人自動離職和其他原因，存在自然減員問題。有不少人在受雇后不滿一年就自動離職；干滿一年后，離職的情況就少了?？紤]到這一因素，設(shè)自然減員率如下表：(表2)分類不熟練半熟練熟練工作不滿一年工作一年以上2510255105現(xiàn)在沒有招工。所有的現(xiàn)有工人都已受雇一年以上。招工 每年新招的各類工人數(shù)熟練工和不熟練工各不超過500，半熟練工不超過800名。再培訓(xùn) 每年可培訓(xùn)200名不熟練工成

5、為半熟練工，每培訓(xùn)一名費(fèi)用開支為4000元。培訓(xùn)半熟練工成為熟練工，培訓(xùn)一名的開支為5000元；培訓(xùn)人數(shù)不能超過所訓(xùn)崗位當(dāng)時(shí)熟練工人數(shù)的1/4.可以將工人降低熟練等級使用，這雖然不需要公司支付什么費(fèi)用，但這樣的工人有50將離職（這一減員要另外加到上述的自然減員上）。解雇 解雇一名不熟練工需付他1600元，解雇一名半熟練工或熟練工要付3000元。超員雇用 全公司可超需要多雇用150人。額外的費(fèi)用每人每年為：不熟練工15000元；半熟練工20000元；熟練工30000元。半日工 不熟練、半熟練和熟練工可以各有不超過50名作為半日工，完成半個(gè)人的生產(chǎn)任務(wù)。這樣做公司對其每人每年支付費(fèi)用為不熟練工50

6、00元，半熟練和熟練工7000元。公司提出的目標(biāo)為解雇的人數(shù)最少。為此應(yīng)如何運(yùn)轉(zhuǎn)？如果公司的目標(biāo)為費(fèi)用最少，能多節(jié)省多少費(fèi)用？導(dǎo)出每年每類崗位所節(jié)省的費(fèi)用。2.1問題的分析這是一個(gè)優(yōu)化問題，它的目標(biāo)函數(shù)有兩個(gè)：第一個(gè)目標(biāo)：公司解雇的人數(shù)最少；要做的決策是關(guān)于各個(gè)月各類工人的人數(shù)安排，即依托現(xiàn)有的工人，通過對公認(rèn)的培訓(xùn)、招聘、降級使用、提升、解雇、額外招工、招收短工的方式來保證在表1中對各個(gè)月熟練、半熟練、不熟練的工人的人數(shù)要求。即：表1分類不熟練半熟練熟練現(xiàn)有人數(shù)第一年需求第二年需求第三年需求20001000500015001400200025001000100015002000在滿足要求的基

7、礎(chǔ)上，求解出人數(shù)的最小值。決策收到以下條件的約束：1、自然減員問題；2、各類工人招收人數(shù)限制；3、培訓(xùn)現(xiàn)有各類工人的人數(shù)限制；4、額外招共工人數(shù)限制；5、解雇工人的限制第二個(gè)目標(biāo)：公司的使用的費(fèi)用最少，求得節(jié)省費(fèi)用。所需的決策為公司功過培訓(xùn)、招聘、降級使用、提升、解雇、額外招工、招收短工的方式來保證在表1（同上）中對各個(gè)月熟練、半熟練、不熟練的工人的人數(shù)要求。在滿足要求的基礎(chǔ)上再求出公司為人力安排所需的費(fèi)用最小值，在本題的費(fèi)用開銷問題上，排除了對一般工人正常工資的發(fā)放，只考慮因?yàn)榕嘤?xùn)、招聘、額外招工、招收短工而為公司帶來的費(fèi)用開銷。按照題目要求及所給條件，將決策變量、目標(biāo)函數(shù)和約束條件用數(shù)學(xué)符

8、號及式子表達(dá)來就可以得到所需最優(yōu)解的人力計(jì)劃。2.2基本假設(shè)問題一假設(shè)：1.題中給出了公司未來三年的人力需求的估計(jì)數(shù)，在現(xiàn)實(shí)人員的安排上可以適當(dāng)?shù)亩嘤诠镜囊螅嗟穆毠た梢宰鲛o退處理2.公司在運(yùn)作過程中存在著自然減員的情況，在這里我們就假定自然減員就只有職工自動離職這一情況。3.公司的職員，不論是新職工，還是具有一定資歷的職工，在被解聘，被降職使用，被提拔和被培訓(xùn)方面擁有著同等機(jī)會。4.各類職工（熟練，半熟練，不熟練）在被解聘，被降職使用，被提拔和被培訓(xùn)方面擁有著同等機(jī)會。問題二假設(shè)：1.保留問題一種的四項(xiàng)假設(shè)同時(shí)在本題目中我們所要考慮的公司的費(fèi)用節(jié)省只是針對于培訓(xùn)，解聘，超員招工，半天工的

9、招聘的費(fèi)用節(jié)省，而正常職員的基本薪金，獎(jiǎng)金等并不包含在內(nèi)。2.在建模中為了保證模型的科學(xué)性和準(zhǔn)確性，所有的決策變量都假設(shè)為非負(fù)實(shí)數(shù)，因而造成的費(fèi)用可以忽略。2.3、定義符號說明Xij：第i年的擁有的j型工人（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）XCij：第i年辭退的j型工人（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）XDij: 第i年招聘短工的j型工人（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）XEij：第i年額外招聘的j型工人（i表示從0到3的

10、年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）XZij：第i年招聘的j型工人（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）XPij：第i年培訓(xùn)的j型工人（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）XJij：第i年降級使用的j型工人（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表半熟練被降使用的，2代表熟練被降為不熟練使用的，3代表熟練降為半熟練使用的）XHij：第i年降級后再提升使用的j型工人(i表示從0到3的年數(shù)：j表示工人的代表類型，1代表原半熟練被降后恢復(fù)使用的，

11、2代表熟練被降為不熟練后提升半熟練的使用的，3代表熟練降為半熟練后恢復(fù)使用的，4、代表熟練被降為不熟練后提升為熟練的使用的）2.4模型的建立 在建模中的相關(guān)問題基本模型1、決策變量：、設(shè)將表1中的數(shù)據(jù)分別設(shè)為Xij（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）X01=2000, X02=1500， X03=1000，X11=1000， X12=1400， X13=1000，X21=500， X22=2000， X23=1500，X31=0， X32=2500， X33=2000設(shè)招工人數(shù)為第一年招不熟練的XZ11,半熟練的XZ12,熟練地XZ13；

12、第二年招不熟練的XZ21,半熟練的XZ22,熟練地XZ23；第三年招不熟練的XZ31,半熟練的XZ32,熟練地XZ33；由已知得: XZi1=500；XZi2=800；XZi3=500 (i=1、2、3)設(shè)培訓(xùn)的人數(shù)：第一年培訓(xùn)不熟練的XP11，半熟練的XP12；第二年培訓(xùn)不熟練的XP21，半熟練的XP22；第三年培訓(xùn)不熟練的XP31，半熟練的XP32；由已知得： XPi1=200,XPi2=1/4Xi1.培訓(xùn)費(fèi)用為 P=（4000XPi1+5000XPi2） (i=1、2、3)設(shè)被降級使用的工人XJij（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表半熟練被降使用的，2代表熟練被降為不熟

13、練使用的，3代表熟練降為半熟練使用的）由已知得： XJij=Xi295%+XZi275%XJi2+ XJi3=Xi195%+XZi390%. (i=1、2、3)設(shè)被降級的工人中又被提升的人數(shù)為Xhij(i表示從0到3的年數(shù)：j表示工人的代表類型，1代表原半熟練被降后恢復(fù)使用的，2代表熟練被降為不熟練后提升半熟練的使用的，3代表熟練降為半熟練后恢復(fù)使用的，4、代表熟練被降為不熟練后提升為熟練的使用的）由已知得：XH1j=0XH21= XJ1150%90%XH22+XH24-XJ1250%90%=0XH23-XJ1350%95%=0XH31-(XJ2150%+XJ1150%90%-XH21) 90

14、%=0XH32+XH34=(XJ2250%+XJ1250%90%-XH22-XH24)90%XH33=(XJ2350%+XJ1350%95%-XH23+XH22)95%注：提升只是對于上一年降級的工人而言設(shè)辭退的工人人數(shù)為XCij（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）由已知得：XC11=200090%-XP11XC12=150095%-XP12XC13=100095%-XP13XC21=100090%-XP21XC22=140095%-XP22XC23=100095%-XP23XC31=50090%-XP31XC32=200095%-XP3

15、2XC33=150095%-XP33辭退費(fèi)用為C=（1600Ci1+3000XCi2+3000XCi3）設(shè)每年招收額外一年各類工人人數(shù)為XEij（i表示從0到3的年數(shù)：j表示工人的代表類型，1代表不熟練的，2代表半熟練的，3代表熟練的）由已知得：XEij=150額外招工的費(fèi)用為E=（1500XEi1+2000XEi2+3000XEi3）設(shè)招收短工分別為XDij由已知得：XDi1，XDi2，XDi3=50招收短工的費(fèi)用為D=（5000XDi1+7000XDi2+7000XDi3）將整個(gè)的決策變量整合，建立綜合的初等線性規(guī)劃模型：200090%+XZ1175%+(XJ11+XJ12)50%-XP1

16、1-XC11+XE11+XD11=1000；150095%+XZ1275%+XJ1350%-XJ11+XP11-XP12-XC12+XE12+XD12=1400；100095%+XZ1390%-XJ12-XJ13+XP12-XC13+XE13+XD13=1000；100090%+XZ1175%+(XJ21+XJ22) 50%-XP21-XC21+XE21+XD21-XH21-XH22-XH24=500；140095%+XZ2275%+XJ2350%-XJ21+XP21-XP22-XC22+XE22+XD22+XH21+XH22-XH23=2000；100095%+XZ2390%+XP22-XJ

17、22-XJ23-XC23+XE23+XD23+XH23+XH24=1500；50090%+XZ3175%+(XJ31+XJ32) 50%-XP31-XC31+XE31+XD31-XH31-XH32-XH34=0；200095%+XZ3275%+XJ3350%-XJ31+XP31-XP32-XC32+XE32+XD32+XH31+XH32-XH33=2500；150095%+XZ3390%+XP32-XJ32-XJ33-XC33+XE33+XD33+XH33+XH34=2000注：同時(shí)以上所有的變量都為非負(fù)的.2、目標(biāo)函數(shù)：對于（1）的要求，盡量少的辭退工人 MinXCij對于（2）的要求，盡量

18、少的費(fèi)用 Min（P+C+E+D）3.1 模型的求解以下是邏輯代碼，并不能再軟件上運(yùn)行，運(yùn)行代碼在附件；Lingo運(yùn)行：MODELMin=XC11+XC12+XC13+XC21+XC22+XC23+XC31+XC32+XC33Min=4000*XP11+5000*XP12+4000*XP21+5000*XP22+4000*XP31+5000*XP32+XC11*1600+3000*XC12+3000*XC13+1600*XC21+3000*XC22+3000*XC23+1600*XC31+3000*XC32+3000*XC33+1500*XE11+2000*XE12+3000*XE13+150

19、0*XE21+2000*XE22+3000*XE23+1500*XE31+2000*XE32+3000*XE33+5000*XD11+7000*XD12+7000*XD13+5000*XD21+7000*XD22+7000*XD23+5000*XD31+7000*XD32+7000*XD33;X01=2000X02=1500X03=1000X11=1000X12=1400X13=1000X21=500X22=2000X23=1500X31=0X32=2500X33=2000XZ11=500XZ21=500XZ31=500XZ12=800XZ22=800XZ32=800XZ13=500XZ23=

20、500XZ33=500XP11=200XP21=200XP31=200XP12=200XP22=200XP32=200XP13=200XP23=200XP33=200XP12=X13/4XP22=X23/4XP32=X33/4XJ11=X12*95%+XZ12*75%XJ21=X22*95%+XZ22*75%XJ31=X32*95%+XZ32*75%XJ12+XJ13=X11*95%+XZ13*90%XJ22+XJ23=X21*95%+XZ23*90%XJ32+XJ33=X31*95%+XZ33*90%XH11=0XH12=0XH13=0XH14=0XH21=XJ11*50%*90%XH22+

21、XH24=XJ12*50%*90%XH23=XJ13*50%*95%XH31=(XJ21*50%+XJ11*50%*90%-XH21)*90%XH32+XH34=(XJ22*50%+XJ12*50%*90%-XH22-XH24)*90%XH33=(XJ23*50%+XJ13*50%*95%-XH23+XH22)*95%XC11=2000*90%-XP11XC12=1500*95%-XP12XC13=1000*95%XC21=1000*90%-XP21XC22=1400*95%-XP22XC23=1000*95%XC31=500*90%-XP31XC32=2000*95%-XP32XC33=15

22、00*95%XE11+XE12+XE13+XE21+XE22+XE23+XE31+XE32+XE33=150XD11+XD21+XD31=50XD12+XD22+XD32=50XD13+XD23+XD33=502000*90%+XZ11*75%+(XJ11+XJ12)*50%-XP11-XC11+XE11+XD11=10001500*95%+XZ12*75%+XJ13*50%-XJ11+XP11-XP12-XC12+XE12+XD12=14001000*95%+XZ13*90%-XJ12-XJ13+XP12-XC13+XE13+XD13=10001000*90%+XZ11*75%+(XJ21+

23、XJ22)*50%-XP21-XC21+XE21+XD21-XH21-XH22-XH24=5001400*95%+XZ22*75%+XJ23*50%-XJ21+XP21-XP22-XC22+XE22+XD22+XH21+XH22-XH23=20001000*95%+XZ23*90%+XP22-XJ22-XJ23-XC23+XE23+XD23+XH23+XH24=1500500*90%+XZ31*75%+(XJ31+XJ32)*50%-XP31-XC31+XE31+XD31-XH31-XH32-XH34=02000*95%+XZ32*75%+XJ33*50%-XJ31+XP31-XP32-XC3

24、2+XE32+XD32+XH31+XH32-XH33=25001500*95%+XZ33*90%+XP32-XJ32-XJ33-XC33+XE33+XD33+XH33+XH34=2000END4.1、模型的評價(jià)與改進(jìn)該模型屬于數(shù)學(xué)規(guī)劃模型，為按公司要求最合理安排人力而設(shè)計(jì)的模型，根據(jù)不同年份，公司對熟練程度不同的工人的安排需求不一樣，因?yàn)闆Q策變量較多，所以分析較難，得考慮的方面也較多。公司每年都需根據(jù)不同的市場經(jīng)濟(jì)的需求，和自身公司人員數(shù)量，采取招聘、培訓(xùn)、降級使用、提升、解雇、額外招工、招收短工等方式進(jìn)行內(nèi)部的調(diào)整。因?yàn)楣久磕陮λ璨煌愋偷墓と擞兴?guī)劃，但要滿足最優(yōu)解的所有條件，所求得的

25、值并不一定是最完美的，下一年的結(jié)果都是根據(jù)上一年的人員再幅度增加或減少的，肯定會存在一定的誤差，所以根據(jù)所求的結(jié)果再進(jìn)行一定的修改，再和公司的現(xiàn)實(shí)情況結(jié)合進(jìn)一步確定，從而求得最優(yōu)解.4.2、個(gè)人體會 個(gè)人心得體會思考并快樂著。通過這次和我的隊(duì)友一起完成學(xué)期的課程設(shè)計(jì)，我從中學(xué)到了很多，主要有兩個(gè)方面：一方面就是通過對問題的深入思考和將自己的所學(xué)賦于實(shí)踐讓我在專業(yè)知識的理解上有了很大的提升；另一方面，通過和隊(duì)友的相互協(xié)作讓我體會到團(tuán)隊(duì)合作的重要性。首先，在關(guān)于規(guī)劃模型方面的知識，自己通過查閱資料和自己的一些體會有以下一些總結(jié)，雖然自己在模型的設(shè)計(jì)和求解中并未將這些知識盡數(shù)用到，但因?yàn)檠芯康男枰?/p>

26、己也認(rèn)真的做了一些總結(jié)，包括以下幾個(gè)方面：（1）線性規(guī)劃問題，其實(shí)我們所做的就是一道特殊的線性規(guī)劃問題，做這類型題應(yīng)該知道的基本知識點(diǎn)是： 將給定的線性規(guī)劃問題化為標(biāo)準(zhǔn)型 并能根據(jù)簡單的實(shí)際問題，建立線性規(guī)劃問題的數(shù)學(xué)模型，并用單純形法求解；同時(shí)要知道1.若線性規(guī)劃問題存在最優(yōu)解，它一定在可行域的某個(gè)頂點(diǎn)得到。 2 .若在兩個(gè)頂點(diǎn)同時(shí)得到最優(yōu)解，則它們連線上的任意一點(diǎn)都是最優(yōu)解，即有無窮多最優(yōu)解。 3 .線性規(guī)劃問題的每個(gè)基可行解對應(yīng)可行域的一個(gè)頂點(diǎn)。 4 .線性規(guī)劃問題的可行解如為最優(yōu)解，則該可行解一定是基可行解。（2）線性規(guī)劃問題和非線性規(guī)劃問題： 主要是要掌握如何將實(shí)際問題轉(zhuǎn)化為線性或非

27、線性模型，并會用LINDO或者LINGO對問題進(jìn)行求解。（3）動態(tài)規(guī)劃問題，要學(xué)會動態(tài)規(guī)劃的遞推方法并要學(xué)會將實(shí)際的問題求解轉(zhuǎn)化為相應(yīng)的LINDO代碼進(jìn)行求解。其次，通過這次課設(shè)讓我了解到自己所學(xué)專業(yè)在實(shí)際生活中的重要性。同時(shí)，也認(rèn)識到了自己的知識面的狹窄，不足以滿足實(shí)際生活的需要。附件（包括代碼和結(jié)果分析）（1） 運(yùn)行代碼：MODEL:min=XC11+XC12+XC13+XC21+XC22+XC23+XC31+XC32+XC33(問題一的目標(biāo)函數(shù))min=4000*XP11+5000*XP12+4000*XP21+5000*XP22+4000*XP31+5000*XP32+XC11*160

28、0+3000*XC12+3000*XC13+1600*XC21+3000*XC22+3000*XC23+1600*XC31+3000*XC32+3000*XC33+1500*XE11+2000*XE12+3000*XE13+1500*XE21+2000*XE22+3000*XE23+1500*XE31+2000*XE32+3000*XE33+5000*XD11+7000*XD12+7000*XD13+5000*XD21+7000*XD22+7000*XD23+5000*XD31+7000*XD32+7000*XD33;（問題二的目標(biāo)函數(shù)）注：兩個(gè)問題的決策變量和約束條件一致，以下代碼一致。X0

29、1=2000;X02=1500;X03=1000;X11=1000;X12=1400;X13=1000;X21=500;X22=2000;X23=1500;X31=0;X32=2500;X33=2000;XZ11=500;XZ21=500;XZ31=500;XZ12=800;XZ22=800;XZ32=800;XZ13=500;XZ23=500;XZ33=500;XP11=200;XP21=200;XP31=200;XP12=250;XP22=375;XP32=500;XJ11-XZ12*0.75=1330;XJ21-X22*0.95-XZ22*0.75=0;XJ31-X32*0.95-XZ3

30、2*0.75=0;XJ12+XJ13-X11*0.95-XZ13*0.9=0;XJ22+XJ23-X21*0.95-XZ23*0.9=0;XJ32+XJ33-X31*0.95-XZ33*0.9=0;XH11=0;XH12=0;XH13=0;XH14=0;XH21-XJ11*0.45=0;XH22+XH24-XJ12*0.45=0;XH23-XJ13*0.48=0;XH31-XJ21*0.45-XJ11*0.41+XH21*0.9=0;XH32+XH34-XJ22*0.45-XJ12*0.41+XH22*0.9+XH24*0.9=0;XH33-XJ23*0.48-XJ13*0.46+XH23*0.95-XH22*0.95=0;XC11+XP11=1800;XC12+XP12=1425;XC13=950;XC21+XP21=900;XC22+XP22=1330;XC23=950;XC31+XP31=450;XC32+XP32=1900;XC33=1425;XE11+XE12+XE13+