AlfaS1= 0.0000000567* 0.7*(( y(7)+273)+(y(1)+273))*((y(7)+273)^2+(y(1)+273)^2);
AlfaT1= 0.0000000567* 0.7*(( y(7)+273)+(y(4)+273))*((y(7)+273)^2+(y(4)+273)^2);

AlfaTr2= 0.0000000567* 0.7*(( y(7)+273)+(y(8)+273))*((y(7)+273)^2+(y(8)+273)^2);
RklS1 =  1/(AlfaS1); RklT1 =  1/(AlfaT1);
RklTr2 = 1/(AlfaTr2* 0.0157  ) ; 

AlfaTr3= 0.0000000567* 0.7*(( y(1)+273)+(y(8)+273))*((y(1)+273)^2+(y(8)+273)^2); 
AlfaTr4= 0.0000000567* 0.7*(( y(4)+273)+(y(8)+273))*((y(4)+273)^2+(y(8)+273)^2); 

RklTr3 = 1/(AlfaTr3* 0.0157  ) ; 
RklTr4 = 1/(AlfaTr4* 0.0157  ) ; 

Tr1 = -( (1/(( RTr1 + RklTr2 )*CTr1)) + (1/((RTr2+Rkl1)*CTr1))  +  (1/(( RTr1 + RklTr3 )*CTr1)) + (1/(( RTr1 + RklTr4 )*CTr1))  ) ;
TrN =    (1/(( RTr1 + RklTr2 )*CTr1)) ; 
TrS =                                                              (1/(( RTr1 + RklTr3 )*CTr1)) ; 
TrT =                                                                                             (1/(( RTr1 + RklTr4 )*CTr1)) ; 

TrS1 =  (1/(( RTr1 + RklTr3 )*C11)) ; 
TrT1 =  (1/(( RTr1 + RklTr4 )*Ct11)) ; 


S11 = -( (1/((RklS1+R11)*C11)) + (1/((R12+R21)*C11)) + (1/(( RTr1 + RklTr3 )*C11))   );
S12 =                            (1/((R12+R21)*C11)) ;
S13 = 0 ; 

S21 =                          (1/((R12+R21)*C22)) ;
S22 = -( (1/((R22+R31)*C22)) + (1/((R12+R21)*C22)) );
S23 =    (1/((R22+R31)*C22)) ; 

S31 = 0 ;
S32 =       (1/((R31+R22)*C33)) ;
S33 = -( (1/((R31+R22)*C33)) + (1/((R32+Rkl1)*C33))    )  ; 


T11 = -( (1/((RklT1+Rt11)*Ct11)) + (1/((Rt12+Rt21)*Ct11))  + (1/(( RTr1 + RklTr4 )*Ct11))        );
T12 =                              (1/((Rt12+Rt21)*Ct11)) ;
T13 = 0 ; 

T21 =                             (1/((Rt12+Rt21)*Ct22)) ;
T22 = -( (1/((Rt22+Rt31)*Ct22)) + (1/((Rt12+Rt21)*Ct22)) );
T23 =    (1/((Rt22+Rt31)*Ct22)) ; 

T31 = 0 ;
T32 =    (1/((Rt31+Rt22)*Ct33)) ;
T33 = -( (1/((Rt31+Rt22)*Ct33)) ) ; 

NG11 = (1/((RklS1 + R11)*CN1))   ; 
NG12 = (1/((RklT1 + Rt11)*CN1))   ; 
TrN1 = (1/(( RTr1 + RklTr2 )*CN1)) ; 
NG13 = -( NG11 + NG12 + TrN1) ; 

S_1  =  1/ (C11  *( RklS1 + R11  ))  ; 
S_2  =  1/ (Ct11 *( RklT1 + Rt11 ))  ; 

yp = [ 
    S11*y(1)  + S12*y(2) + S13*y(3) + 0         + 0        + 0        + S_1*y(7)  + TrS1*y(8) + 0     
    S21*y(1)  + S22*y(2) + S23*y(3) + 0         + 0        + 0        + 0         + 0         + 0           
    S31*y(1)  + S32*y(2) + S33*y(3) + 0         + 0        + 0        + 0         + 0         + 0
    0         + 0        + 0        + T11*y(4)  + T12*y(5) + T13*y(6) + S_2*y(7)  + TrT1*y(8) + 0
    0         + 0        + 0        + T21*y(4)  + T22*y(5) + T23*y(6) + 0         + 0         + 0
    0         + 0        + 0        + T31*y(4)  + T32*y(5) + T33*y(6) + 0         + 0         + 0
    NG11*y(1) + 0        + 0        + NG12*y(4) + 0        + 0        + NG13*y(7) + TrN1*y(8) + (q1/CN1)
    TrS*y(1)  + 0        + 0        + TrT*y(4)  + 0        + 0        + TrN*y(7)  + Tr1*y(8)  + 0
] ;