AlfaS1= 0.0000000567* 0.7*(( y(44)+273)+(y(1)+273))*((y(44)+273)^2+(y(1)+273)^2); AlfaT1= 0.0000000567* 0.7*(( y(44)+273)+(y(25)+273))*((y(44)+273)^2+(y(25)+273)^2); AlfaS2= 0.0000000567* 0.7*(( y(45)+273)+(y(5)+273))*((y(45)+273)^2+(y(5)+273)^2); AlfaT2= 0.0000000567* 0.7*(( y(45)+273)+(y(28)+273))*((y(45)+273)^2+(y(28)+273)^2); AlfaS3= 0.0000000567* 0.7*(( y(46)+273)+(y(9)+273))*((y(46)+273)^2+(y(9)+273)^2); AlfaT3= 0.0000000567* 0.7*(( y(46)+273)+(y(31)+273))*((y(46)+273)^2+(y(31)+273)^2); AlfaS4= 0.0000000567* 0.7*(( y(47)+273)+(y(13)+273))*((y(47)+273)^2+(y(13)+273)^2); AlfaT4= 0.0000000567* 0.7*(( y(47)+273)+(y(34)+273))*((y(47)+273)^2+(y(34)+273)^2); AlfaS5= 0.0000000567* 0.7*(( y(48)+273)+(y(17)+273))*((y(48)+273)^2+(y(17)+273)^2); AlfaT5= 0.0000000567* 0.7*(( y(48)+273)+(y(37)+273))*((y(48)+273)^2+(y(37)+273)^2); AlfaS6= 0.0000000567* 0.7*(( y(49)+273)+(y(21)+273))*((y(49)+273)^2+(y(21)+273)^2); AlfaT6= 0.0000000567* 0.7*(( y(49)+273)+(y(40)+273))*((y(49)+273)^2+(y(40)+273)^2); RklS1 = 1/(AlfaS1*1); RklT1 = 1/(AlfaT1*0.81); RklS2 = 1/(AlfaS2*1); RklT2 = 1/(AlfaT2*0.81); RklS3 = 1/(AlfaS3*1); RklT3 = 1/(AlfaT3*0.81); RklS4 = 1/(AlfaS4*1); RklT4 = 1/(AlfaT4*0.81); RklS5 = 1/(AlfaS5*1); RklT5 = 1/(AlfaT5*0.81); RklS6 = 1/(AlfaS6*1); RklT6 = 1/(AlfaT6*0.81); mas_y(1,1) = y(4); mas_y(1,2) = y(8); mas_y(1,3) = y(12); mas_y(1,4) = y(16); mas_y(1,5) = y(20); mas_y(1,6) = y(24); for kk = 1:6 temp_t = mas_y(1, kk) ; if (temp_t <= 30) mas_Rkl(1, kk) = 9.8 ; end if (temp_t > 30 & temp_t <= 40) mas_Rkl(1, kk) = 10.4 ; end if (temp_t > 40 & temp_t <= 50) mas_Rkl(1, kk) = 11.2 ; end if (temp_t > 50 & temp_t <= 60) mas_Rkl(1, kk) = 12.1 ; end if (temp_t > 60 & temp_t <= 80) mas_Rkl(1, kk) = 13.3 ; end if (temp_t > 80 & temp_t <= 100) mas_Rkl(1, kk) = 14.5 ; end if (temp_t > 100 & temp_t <= 150) mas_Rkl(1, kk) = 17.4 ; end if (temp_t > 150 & temp_t <= 200) mas_Rkl(1, kk) = 20.2 ; end if (temp_t > 200 & temp_t <= 300) mas_Rkl(1, kk) = 27.0 ; end if (temp_t > 300 & temp_t <= 400) mas_Rkl(1, kk) = 35.0 ; end if (temp_t > 400) mas_Rkl(1, kk) = 36.0 ; end end Rkl1 = 1/(mas_Rkl(1,1) * 3.3856 ); Rkl2 = 1/(mas_Rkl(1,2) * 3.3856 ); Rkl3 = 1/(mas_Rkl(1,3) * 3.3856 ); Rkl4 = 1/(mas_Rkl(1,4) * 3.3856 ); Rkl5 = 1/(mas_Rkl(1,5) * 3.3856 ); Rkl6 = 1/(mas_Rkl(1,6) * 3.3856 ); Landa = 0.520+0.349*0.001*y(1); R11 = (1/Landa)*(0.01); R12 = (1/Landa)*(100); R13 = (1/Landa)*(0.00961169); R14 = (1/Landa)*(100); R15 = (1/Landa)*(100); R16 = (1/Landa)*(100); Landa = 0.520+0.349*0.001*y(2); R17 = (1/Landa)*(0.192234); R18 = (1/Landa)*(5); R19 = (1/Landa)*(0.0991867); R110 = (1/Landa)*(5); R111 = (1/Landa)*(5); R112 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(3); R113 = (1/Landa)*(0.0991867); R114 = (1/Landa)*(5); R115 = (1/Landa)*(0.0603792); R116 = (1/Landa)*(5); R117 = (1/Landa)*(5); R118 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(4); R119 = (1/Landa)*(0.00301896); R120 = (1/Landa)*(100); R121 = (1/Landa)*(0.00295369); R122 = (1/Landa)*(100); R123 = (1/Landa)*(100); R124 = (1/Landa)*(100); C11 = 10.2013 * (837+0.41*y(1)); C12 = 300.347 * (837+0.41*y(2)); C13 = 527.547 * (837+0.41*y(3)); C14 = 33.4893 * (837+0.41*y(4)); Landa = 0.520+0.349*0.001*y(5); R21 = (1/Landa)*(0.01); R22 = (1/Landa)*(100); R23 = (1/Landa)*(0.00961169); R24 = (1/Landa)*(100); R25 = (1/Landa)*(100); R26 = (1/Landa)*(100); Landa = 0.520+0.349*0.001*y(6); R27 = (1/Landa)*(0.192234); R28 = (1/Landa)*(5); R29 = (1/Landa)*(0.0991867); R210 = (1/Landa)*(5); R211 = (1/Landa)*(5); R212 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(7); R213 = (1/Landa)*(0.0991867); R214 = (1/Landa)*(5); R215 = (1/Landa)*(0.0603792); R216 = (1/Landa)*(5); R217 = (1/Landa)*(5); R218 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(8); R219 = (1/Landa)*(0.00301896); R220 = (1/Landa)*(100); R221 = (1/Landa)*(0.00295369); R222 = (1/Landa)*(100); R223 = (1/Landa)*(100); R224 = (1/Landa)*(100); C21 = 10.2013 * (837+0.41*y(5)); C22 = 300.347 * (837+0.41*y(6)); C23 = 527.547 * (837+0.41*y(7)); C24 = 33.4893 * (837+0.41*y(8)); Landa = 0.520+0.349*0.001*y(9); R31 = (1/Landa)*(0.01); R32 = (1/Landa)*(100); R33 = (1/Landa)*(0.00961169); R34 = (1/Landa)*(100); R35 = (1/Landa)*(100); R36 = (1/Landa)*(100); Landa = 0.520+0.349*0.001*y(10); R37 = (1/Landa)*(0.192234); R38 = (1/Landa)*(5); R39 = (1/Landa)*(0.0991867); R310 = (1/Landa)*(5); R311 = (1/Landa)*(5); R312 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(11); R313 = (1/Landa)*(0.0991867); R314 = (1/Landa)*(5); R315 = (1/Landa)*(0.0603792); R316 = (1/Landa)*(5); R317 = (1/Landa)*(5); R318 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(12); R319 = (1/Landa)*(0.00301896); R320 = (1/Landa)*(100); R321 = (1/Landa)*(0.00295369); R322 = (1/Landa)*(100); R323 = (1/Landa)*(100); R324 = (1/Landa)*(100); C31 = 10.2013 * (837+0.41*y(9)); C32 = 300.347 * (837+0.41*y(10)); C33 = 527.547 * (837+0.41*y(11)); C34 = 33.4893 * (837+0.41*y(12)); Landa = 0.520+0.349*0.001*y(13); R41 = (1/Landa)*(0.01); R42 = (1/Landa)*(100); R43 = (1/Landa)*(0.00961169); R44 = (1/Landa)*(100); R45 = (1/Landa)*(100); R46 = (1/Landa)*(100); Landa = 0.520+0.349*0.001*y(14); R47 = (1/Landa)*(0.192234); R48 = (1/Landa)*(5); R49 = (1/Landa)*(0.0991867); R410 = (1/Landa)*(5); R411 = (1/Landa)*(5); R412 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(15); R413 = (1/Landa)*(0.0991867); R414 = (1/Landa)*(5); R415 = (1/Landa)*(0.0603792); R416 = (1/Landa)*(5); R417 = (1/Landa)*(5); R418 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(16); R419 = (1/Landa)*(0.00301896); R420 = (1/Landa)*(100); R421 = (1/Landa)*(0.00295369); R422 = (1/Landa)*(100); R423 = (1/Landa)*(100); R424 = (1/Landa)*(100); C41 = 10.2013 * (837+0.41*y(13)); C42 = 300.347 * (837+0.41*y(14)); C43 = 527.547 * (837+0.41*y(15)); C44 = 33.4893 * (837+0.41*y(16)); Landa = 0.520+0.349*0.001*y(17); R51 = (1/Landa)*(0.01); R52 = (1/Landa)*(100); R53 = (1/Landa)*(0.00961169); R54 = (1/Landa)*(100); R55 = (1/Landa)*(100); R56 = (1/Landa)*(100); Landa = 0.520+0.349*0.001*y(18); R57 = (1/Landa)*(0.192234); R58 = (1/Landa)*(5); R59 = (1/Landa)*(0.0991867); R510 = (1/Landa)*(5); R511 = (1/Landa)*(5); R512 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(19); R513 = (1/Landa)*(0.0991867); R514 = (1/Landa)*(5); R515 = (1/Landa)*(0.0603792); R516 = (1/Landa)*(5); R517 = (1/Landa)*(5); R518 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(20); R519 = (1/Landa)*(0.00301896); R520 = (1/Landa)*(100); R521 = (1/Landa)*(0.00295369); R522 = (1/Landa)*(100); R523 = (1/Landa)*(100); R524 = (1/Landa)*(100); C51 = 10.2013 * (837+0.41*y(17)); C52 = 300.347 * (837+0.41*y(18)); C53 = 527.547 * (837+0.41*y(19)); C54 = 33.4893 * (837+0.41*y(20)); Landa = 0.520+0.349*0.001*y(21); R61 = (1/Landa)*(0.01); R62 = (1/Landa)*(100); R63 = (1/Landa)*(0.00961169); R64 = (1/Landa)*(100); R65 = (1/Landa)*(100); R66 = (1/Landa)*(100); Landa = 0.520+0.349*0.001*y(22); R67 = (1/Landa)*(0.192234); R68 = (1/Landa)*(5); R69 = (1/Landa)*(0.0991867); R610 = (1/Landa)*(5); R611 = (1/Landa)*(5); R612 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(23); R613 = (1/Landa)*(0.0991867); R614 = (1/Landa)*(5); R615 = (1/Landa)*(0.0603792); R616 = (1/Landa)*(5); R617 = (1/Landa)*(5); R618 = (1/Landa)*(5); Landa = 0.520+0.349*0.001*y(24); R619 = (1/Landa)*(0.00301896); R620 = (1/Landa)*(100); R621 = (1/Landa)*(0.00295369); R622 = (1/Landa)*(100); R623 = (1/Landa)*(100); R624 = (1/Landa)*(100); C61 = 10.2013 * (837+0.41*y(21)); C62 = 300.347 * (837+0.41*y(22)); C63 = 527.547 * (837+0.41*y(23)); C64 = 33.4893 * (837+0.41*y(24)); Landa = 0.520+0.349*0.001*y(25); Rt11 = (1/Landa)*(0.138889); Rt12 = (1/Landa)*(8.88889); Rt13 = (1/Landa)*(0.246914); Rt14 = (1/Landa)*(8.88889); Rt15 = (1/Landa)*(8.88889); Rt16 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(26); Rt17 = (1/Landa)*(0.246914); Rt18 = (1/Landa)*(8.88889); Rt19 = (1/Landa)*(0.555556); Rt110 = (1/Landa)*(8.88889); Rt111 = (1/Landa)*(8.88889); Rt112 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(27); Rt113 = (1/Landa)*(0.555556); Rt114 = (1/Landa)*(8.88889); Rt115 = (1/Landa)*(2.22222); Rt116 = (1/Landa)*(8.88889); Rt117 = (1/Landa)*(8.88889); Rt118 = (1/Landa)*(8.88889); Ct11 = 70.2422 * (837+0.41*y(25)); Ct12 = 36.0703 * (837+0.41*y(26)); Ct13 = 13.2891 * (837+0.41*y(27)); Landa = 0.520+0.349*0.001*y(28); Rt21 = (1/Landa)*(0.138889); Rt22 = (1/Landa)*(8.88889); Rt23 = (1/Landa)*(0.246914); Rt24 = (1/Landa)*(8.88889); Rt25 = (1/Landa)*(8.88889); Rt26 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(29); Rt27 = (1/Landa)*(0.246914); Rt28 = (1/Landa)*(8.88889); Rt29 = (1/Landa)*(0.555556); Rt210 = (1/Landa)*(8.88889); Rt211 = (1/Landa)*(8.88889); Rt212 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(30); Rt213 = (1/Landa)*(0.555556); Rt214 = (1/Landa)*(8.88889); Rt215 = (1/Landa)*(2.22222); Rt216 = (1/Landa)*(8.88889); Rt217 = (1/Landa)*(8.88889); Rt218 = (1/Landa)*(8.88889); Ct21 = 70.2422 * (837+0.41*y(28)); Ct22 = 36.0703 * (837+0.41*y(29)); Ct23 = 13.2891 * (837+0.41*y(30)); Landa = 0.520+0.349*0.001*y(31); Rt31 = (1/Landa)*(0.138889); Rt32 = (1/Landa)*(8.88889); Rt33 = (1/Landa)*(0.246914); Rt34 = (1/Landa)*(8.88889); Rt35 = (1/Landa)*(8.88889); Rt36 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(32); Rt37 = (1/Landa)*(0.246914); Rt38 = (1/Landa)*(8.88889); Rt39 = (1/Landa)*(0.555556); Rt310 = (1/Landa)*(8.88889); Rt311 = (1/Landa)*(8.88889); Rt312 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(33); Rt313 = (1/Landa)*(0.555556); Rt314 = (1/Landa)*(8.88889); Rt315 = (1/Landa)*(2.22222); Rt316 = (1/Landa)*(8.88889); Rt317 = (1/Landa)*(8.88889); Rt318 = (1/Landa)*(8.88889); Ct31 = 70.2422 * (837+0.41*y(31)); Ct32 = 36.0703 * (837+0.41*y(32)); Ct33 = 13.2891 * (837+0.41*y(33)); Landa = 0.520+0.349*0.001*y(34); Rt41 = (1/Landa)*(0.138889); Rt42 = (1/Landa)*(8.88889); Rt43 = (1/Landa)*(0.246914); Rt44 = (1/Landa)*(8.88889); Rt45 = (1/Landa)*(8.88889); Rt46 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(35); Rt47 = (1/Landa)*(0.246914); Rt48 = (1/Landa)*(8.88889); Rt49 = (1/Landa)*(0.555556); Rt410 = (1/Landa)*(8.88889); Rt411 = (1/Landa)*(8.88889); Rt412 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(36); Rt413 = (1/Landa)*(0.555556); Rt414 = (1/Landa)*(8.88889); Rt415 = (1/Landa)*(2.22222); Rt416 = (1/Landa)*(8.88889); Rt417 = (1/Landa)*(8.88889); Rt418 = (1/Landa)*(8.88889); Ct41 = 70.2422 * (837+0.41*y(34)); Ct42 = 36.0703 * (837+0.41*y(35)); Ct43 = 13.2891 * (837+0.41*y(36)); Landa = 0.520+0.349*0.001*y(37); Rt51 = (1/Landa)*(0.138889); Rt52 = (1/Landa)*(8.88889); Rt53 = (1/Landa)*(0.246914); Rt54 = (1/Landa)*(8.88889); Rt55 = (1/Landa)*(8.88889); Rt56 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(38); Rt57 = (1/Landa)*(0.246914); Rt58 = (1/Landa)*(8.88889); Rt59 = (1/Landa)*(0.555556); Rt510 = (1/Landa)*(8.88889); Rt511 = (1/Landa)*(8.88889); Rt512 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(39); Rt513 = (1/Landa)*(0.555556); Rt514 = (1/Landa)*(8.88889); Rt515 = (1/Landa)*(2.22222); Rt516 = (1/Landa)*(8.88889); Rt517 = (1/Landa)*(8.88889); Rt518 = (1/Landa)*(8.88889); Ct51 = 70.2422 * (837+0.41*y(37)); Ct52 = 36.0703 * (837+0.41*y(38)); Ct53 = 13.2891 * (837+0.41*y(39)); Landa = 0.520+0.349*0.001*y(40); Rt61 = (1/Landa)*(0.138889); Rt62 = (1/Landa)*(8.88889); Rt63 = (1/Landa)*(0.246914); Rt64 = (1/Landa)*(8.88889); Rt65 = (1/Landa)*(8.88889); Rt66 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(41); Rt67 = (1/Landa)*(0.246914); Rt68 = (1/Landa)*(8.88889); Rt69 = (1/Landa)*(0.555556); Rt610 = (1/Landa)*(8.88889); Rt611 = (1/Landa)*(8.88889); Rt612 = (1/Landa)*(8.88889); Landa = 0.520+0.349*0.001*y(42); Rt613 = (1/Landa)*(0.555556); Rt614 = (1/Landa)*(8.88889); Rt615 = (1/Landa)*(2.22222); Rt616 = (1/Landa)*(8.88889); Rt617 = (1/Landa)*(8.88889); Rt618 = (1/Landa)*(8.88889); Ct61 = 70.2422 * (837+0.41*y(40)); Ct62 = 36.0703 * (837+0.41*y(41)); Ct63 = 13.2891 * (837+0.41*y(42)); Ct_center = 11.3906 * (837+0.41*y(43)); Landa_center = 0.520+0.349*0.001 * y(43); Rt119 = (1/Landa_center)*(4.44444); Rt219 = (1/Landa_center)*(4.44444); Rt319 = (1/Landa_center)*(4.44444); Rt419 = (1/Landa_center)*(4.44444); Rt519 = (1/Landa_center)*(4.44444); Rt619 = (1/Landa_center)*(4.44444);