AlfaL1 = 3.969e-008*(( y(1)+273)+(293))*((y(1)+273)^2+(85849) );
AlfaL2 = 3.969e-008*(( y(4)+273)+(293))*((y(4)+273)^2+(85849) );
AlfaL3 = 3.969e-008*(( y(7)+273)+(293))*((y(7)+273)^2+(85849) );
AlfaL4 = 3.969e-008*(( y(10)+273)+(293))*((y(10)+273)^2+(85849) );
AlfaL5 = 3.969e-008*(( y(13)+273)+(293))*((y(13)+273)^2+(85849) );
AlfaL6 = 3.969e-008*(( y(16)+273)+(293))*((y(16)+273)^2+(85849) );
AlfaK(1,6) = 0 ; 
Temp_powyrhnost (1,6) = 0 ; 
Temp_powyrhnost (1,1) = y(1); 
Temp_powyrhnost (1,2) = y(4); 
Temp_powyrhnost (1,3) = y(7); 
Temp_powyrhnost (1,4) = y(10); 
Temp_powyrhnost (1,5) = y(13); 
Temp_powyrhnost (1,6) = y(16); 
for k = 1:6
 	 Tpow = Temp_powyrhnost (1,k) ; 

 	 if (Tpow < -137.15 ) 
 		 Pran = 0.77;
 		 koef_Landa = 0.00925 ;
 		 koef_my = 0.000006924 ;
 		 koef_ro = 3.601 ;
 		 koef_cp = 1027 ; 
 		 koef_beta = 0.01 ;
 	elseif (Tpow >= -137.15 & Tpow < -87.15 ) 
 		 Pran = 0.77;
 		 koef_Landa = 0.00925 ;
 		 koef_my = 0.000006924 ;
 		 koef_ro = 3.601 ;
 		 koef_cp = 1027 ; 
 		 koef_beta = 0.01 ;
 	 elseif (Tpow >= -87.15 & Tpow < -37.15 ) 
 		 Pran = 0.753;
 		 koef_Landa = 0.01374 ;
 		 koef_my = 0.000010283 ;
 		 koef_ro = 2.3675 ;
 		 koef_cp = 1010 ; 
 		 koef_beta = 0.00667 ;
 	 elseif (Tpow >= -37.15 & Tpow < 12.85 ) 
 		 Pran = 0.739 ;
 		 koef_Landa = 0.01809 ;
 		 koef_my = 0.000013289 ;
 		 koef_ro = 1.7684 ;
 		 koef_cp = 1006 ; 
 		 koef_beta = 0.005 ;
 	 elseif (Tpow >= 12.85 & Tpow < 62.85 ) 
 		 Pran = 0.722;
 		 koef_Landa = 0.02227 ;
 		 koef_my = 0.00001488 ;
 		 koef_ro = 1.4128 ;
 		 koef_cp = 1005 ; 
 		 koef_beta = 0.004 ;
 	 elseif (Tpow >= 62.85 & Tpow < 112.85 ) 
 		 Pran = 0.708;
 		 koef_Landa = 0.02624 ;
 		 koef_my = 0.00001983 ;
 		 koef_ro = 1.1774 ;
 		 koef_cp = 1006 ; 
 		 koef_beta = 0.00333 ;
 	 elseif (Tpow >= 112.85 & Tpow < 162.85 ) 
 		 Pran = 0.697;
 		 koef_Landa = 0.03003 ;
 		 koef_my = 0.00002075 ;
 		 koef_ro = 0.998 ;
 		 koef_cp = 1009 ; 
 		 koef_beta = 0.00286 ;
 	 elseif (Tpow >= 162.85 & Tpow < 212.85 ) 
 		 Pran = 0.689;
 		 koef_Landa = 0.03365 ;
 		 koef_my = 0.00002286 ;
 		 koef_ro = 0.8826 ;
 		 koef_cp = 1014 ; 
 		 koef_beta = 0.0025 ;
 	 elseif (Tpow >= 212.85 & Tpow < 262.85 ) 
 		 Pran = 0.683;
 		 koef_Landa = 0.03707 ;
 		 koef_my = 0.00002484 ;
 		 koef_ro = 0.7833 ;
 		 koef_cp = 1021 ; 
 		 koef_beta = 0.00222 ;
 	 elseif (Tpow >= 262.85 & Tpow < 312.85 ) 
 		 Pran = 0.680;
 		 koef_Landa = 0.04038 ;
 		 koef_my = 0.00002671 ;
 		 koef_ro = 0.7048 ;
 		 koef_cp = 1030 ; 
 		 koef_beta = 0.002 ;
 	 elseif (Tpow >= 312.85 & Tpow < 362.85 ) 
 		 Pran = 0.680;
 		 koef_Landa = 0.04360 ;
 		 koef_my = 0.00002848 ;
 		 koef_ro = 0.6423 ;
 		 koef_cp = 1039 ; 
 		 koef_beta = 0.00182 ;
 	 elseif (Tpow >= 362.85 & Tpow < 412.85 ) 
 		 Pran = 0.680;
 		 koef_Landa = 0.04659 ;
 		 koef_my = 0.00003018 ;
 		 koef_ro = 0.5879 ;
 		 koef_cp = 1055 ; 
 		 koef_beta = 0.00167 ;
 	 elseif (Tpow >= 412.85 & Tpow < 462.85 ) 
 		 Pran = 0.682;
 		 koef_Landa = 0.04953 ;
 		 koef_my = 0.00003177 ;
 		 koef_ro = 0.543 ;
 		 koef_cp = 1064 ; 
 		 koef_beta = 0.00154 ;
 	 elseif (Tpow >= 462.85 & Tpow < 512.85 ) 
 		 Pran = 0.684;
 		 koef_Landa = 0.0523 ;
 		 koef_my = 0.00003332 ;
 		 koef_ro = 0.503 ;
 		 koef_cp = 1075 ; 
 		 koef_beta = 0.00143 ;
 	 elseif (Tpow >= 512.85 & Tpow < 562.85 ) 
 		 Pran = 0.686;
 		 koef_Landa = 0.05509 ;
 		 koef_my = 0.00003481 ;
 		 koef_ro = 0.4709 ;
 		 koef_cp = 1086 ; 
 		 koef_beta = 0.00133 ;
 	 elseif (Tpow >= 562.85 & Tpow < 612.85 ) 
 		 Pran = 0.689;
 		 koef_Landa = 0.05779 ;
 		 koef_my = 0.00003625 ;
 		 koef_ro = 0.4405 ;
 		 koef_cp = 1098 ; 
 		 koef_beta = 0.00125 ;
 	 elseif (Tpow >= 612.85 & Tpow < 662.85 ) 
 		 Pran = 0.692;
 		 koef_Landa = 0.06028 ;
 		 koef_my = 0.00003765 ;
 		 koef_ro = 0.4149 ;
 		 koef_cp = 1110 ; 
 		 koef_beta = 0.00118 ;
 	 elseif (Tpow >= 662.85 & Tpow < 712.85 ) 
 		 Pran = 0.696;
 		 koef_Landa = 0.06279 ;
 		 koef_my = 0.00003899 ;
 		 koef_ro = 0.3925 ;
 		 koef_cp = 1121 ; 
 		 koef_beta = 0.00111 ;
 	 elseif (Tpow >= 712.85 & Tpow < 762.85 ) 
 		 Pran = 0.699;
 		 koef_Landa = 0.06525 ;
 		 koef_my = 0.00004023 ;
 		 koef_ro = 0.3716 ;
 		 koef_cp = 1132 ; 
 		 koef_beta = 0.00105 ;
 	 elseif (Tpow >= 762.85 & Tpow < 862.85 ) 
 		 Pran = 0.702;
 		 koef_Landa = 0.06752 ;
 		 koef_my = 0.00004152 ;
 		 koef_ro = 0.3524 ;
 		 koef_cp = 1142 ; 
 		 koef_beta = 0.001 ;
 	 elseif (Tpow >= 862.85 & Tpow < 962.85 ) 
 		 Pran = 0.704;
 		 koef_Landa = 0.0732 ;
 		 koef_my = 0.0000444 ;
 		 koef_ro = 0.3204 ;
 		 koef_cp = 1160 ; 
 		 koef_beta = 0.000909 ;
 	 elseif (Tpow >= 962.85 & Tpow < 1062.85 ) 
 		 Pran = 0.707;
 		 koef_Landa = 0.0782 ;
 		 koef_my = 0.0000469 ;
 		 koef_ro = 0.2947 ;
 		 koef_cp = 1179 ; 
 		 koef_beta = 0.000833 ;
 	 elseif (Tpow >= 1062.85 & Tpow < 1162.85 ) 
 		 Pran = 0.705;
 		 koef_Landa = 0.0837 ;
 		 koef_my = 0.0000493 ;
 		 koef_ro = 0.2707 ;
 		 koef_cp = 1197 ; 
 		 koef_beta = 0.000769 ;
 	 elseif (Tpow >= 1162.85 & Tpow < 1262.85 ) 
 		 Pran = 0.705;
 		 koef_Landa = 0.0891 ;
 		 koef_my = 0.0000517 ;
 		 koef_ro = 0.2515 ;
 		 koef_cp = 1214 ; 
 		 koef_beta = 0.000714 ;
 	 elseif (Tpow >= 1262.85 & Tpow < 1362.85 ) 
 		 Pran = 0.705;
 		 koef_Landa = 0.0946 ;
 		 koef_my = 0.000054 ;
 		 koef_ro = 0.2355 ;
 		 koef_cp = 1230 ; 
 		 koef_beta = 0.000667 ;
 	 elseif (Tpow >= 1362.85 & Tpow < 1462.85 ) 
 		 Pran = 0.705;
 		 koef_Landa = 0.1 ;
 		 koef_my = 0.0000563 ;
 		 koef_ro = 0.2211 ;
 		 koef_cp = 1248 ; 
 		 koef_beta = 0.000625 ;
 	 elseif (Tpow >= 1462.85 & Tpow < 1562.85 ) 
 		 Pran = 0.705;
 		 koef_Landa = 0.105 ;
 		 koef_my = 0.0000585 ;
 		 koef_ro = 0.2082 ;
 		 koef_cp = 1267 ; 
 		 koef_beta = 0.000588 ;
 	 elseif (Tpow >= 1562.85 & Tpow < 1662.85 ) 
 		 Pran = 0.704;
 		 koef_Landa = 0.111 ;
 		 koef_my = 0.0000607 ;
 		 koef_ro = 0.197 ;
 		 koef_cp = 1287 ; 
 		 koef_beta = 0.000556 ;
 	 elseif (Tpow >= 1662.85 & Tpow < 1762.85 ) 
 		 Pran = 0.704;
 		 koef_Landa = 0.177 ;
 		 koef_my = 0.0000629 ;
 		 koef_ro = 0.1858 ;
 		 koef_cp = 1309 ; 
 		 koef_beta = 0.000526 ;
 	 elseif (Tpow >= 1762.85 & Tpow < 1862.85 ) 
 		 Pran = 0.702;
 		 koef_Landa = 0.124 ;
 		 koef_my = 0.000065 ;
 		 koef_ro = 0.1762 ;
 		 koef_cp = 1338 ; 
 		 koef_beta = 0.0005 ;
 	 elseif (Tpow >= 1862.85 & Tpow < 1962.85 ) 
 		 Pran = 0.7;
 		 koef_Landa = 0.131 ;
 		 koef_my = 0.0000672 ;
 		 koef_ro = 0.1682 ;
 		 koef_cp = 1372 ; 
 		 koef_beta = 0.000476 ;
 	 elseif (Tpow >= 1962.85 & Tpow < 2062.85 ) 
 		 Pran = 0.707;
 		 koef_Landa = 0.139 ;
 		 koef_my = 0.0000693 ;
 		 koef_ro = 0.1602 ;
 		 koef_cp = 1419 ; 
 		 koef_beta = 0.000454 ;
 	 elseif (Tpow >= 2062.85 & Tpow < 2162.85 ) 
 		 Pran = 0.71;
 		 koef_Landa = 0.149 ;
 		 koef_my = 0.0000714 ;
 		 koef_ro = 0.1538 ;
 		 koef_cp = 1482 ; 
 		 koef_beta = 0.000435 ;
 	 elseif (Tpow >= 2162.85 & Tpow < 2262.85 ) 
 		 Pran = 0.718;
 		 koef_Landa = 0.161 ;
 		 koef_my = 0.0000735 ;
 		 koef_ro = 0.1458 ;
 		 koef_cp = 1574 ; 
 		 koef_beta = 0.000417 ;
 	 else (Tpow >= 2262.85 ) 
 		 Pran = 0.73;
 		 koef_Landa = 0.175 ;
 		 koef_my = 0.0000757 ;
 		 koef_ro = 0.1394 ;
 		 koef_cp = 1688 ; 
 		 koef_beta = 0.0004 ;
 	 end
 	 delta_t = ( Tpow - TayOS) ; 

 	 if ( k <= 4 ) 
 		 koef_X = 0.8;
 		 Pran_sqrt = sqrt (Pran) ;
 		 koef_K = (1 + (1+(1/Pran_sqrt))^2 )^ -0.25 ; 
 		 Grastgof = (9.81*koef_beta*(koef_ro^2)*delta_t*(koef_X^3))/(koef_my^2) ;
 		 Nuselt = 0.8 * ((Grastgof*Pran)^0.25) * koef_K ;
 		 alfa_ = (Nuselt*koef_Landa)/koef_X ;
 	 end 

 	 if ( k == 5 )
 		 koef_X = 0.8;
 		 Grastgof = (9.81*koef_beta*(koef_ro^2)*delta_t*(koef_X^3))/(koef_my^2) ;
 		 Nuselt = 0.14 * ((Grastgof*Pran)^0.333) ;
 		 alfa_ = (Nuselt*koef_Landa)/koef_X ;
 	 end 

 	 if ( k == 6 )
 		 koef_X = 0.8;
 		 Grastgof = (9.81*koef_beta*(koef_ro^2)*delta_t*(koef_X^3))/(koef_my^2) ;
 		 Nuselt = 0.27 * ((Grastgof*Pran)^0.25) ;
 		 alfa_ = (Nuselt*koef_Landa)/koef_X ;
 	 end 
 	 AlfaK(1,k) = alfa_  ; 
end  

RklT1 = 1/( (AlfaK(1,1)+AlfaL1) * 0.64) ;
RklT2 = 1/( (AlfaK(1,2)+AlfaL2) * 0.64) ;
RklT3 = 1/( (AlfaK(1,3)+AlfaL3) * 0.64) ;
RklT4 = 1/( (AlfaK(1,4)+AlfaL4) * 0.64) ;
RklT5 = 1/( (AlfaK(1,5)+AlfaL5) * 0.64) ;
RklT6 = 1/( (AlfaK(1,6)+AlfaL6) * 0.64) ;