StA = [ S1(1)*y(1)+ S1(2)*y(2)+ 0+ 0+ S2(1)*y(5)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S4(1)*y(13)+ 0+ 0+ 0+ S5(1)*y(17)+ 0+ 0+ 0+ S6(1)*y(21)+ 0+ 0+ 0 S7(1)*y(1)+ S7(2)*y(2)+ S7(3)*y(3)+ 0+ 0+ S8(2)*y(6)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S10(2)*y(14)+ 0+ 0+ 0+ S11(2)*y(18)+ 0+ 0+ 0+ S12(2)*y(22)+ 0+ 0 0+ S13(2)*y(2)+ S13(3)*y(3)+ S13(4)*y(4)+ 0+ 0+ S14(3)*y(7)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S16(3)*y(15)+ 0+ 0+ 0+ S17(3)*y(19)+ 0+ 0+ 0+ S18(3)*y(23)+ 0 0+ 0+ S19(3)*y(3)+ S19(4)*y(4)+ 0+ 0+ 0+ S20(4)*y(8)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S22(4)*y(16)+ 0+ 0+ 0+ S23(4)*y(20)+ 0+ 0+ 0+ S24(4)*y(24) S25(1)*y(1)+ 0+ 0+ 0+ S26(1)*y(5)+ S26(2)*y(6)+ 0+ 0+ S27(1)*y(9)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S29(1)*y(17)+ 0+ 0+ 0+ S30(1)*y(21)+ 0+ 0+ 0 0+ S31(2)*y(2)+ 0+ 0+ S32(1)*y(5)+ S32(2)*y(6)+ S32(3)*y(7)+ 0+ 0+ S33(2)*y(10)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S35(2)*y(18)+ 0+ 0+ 0+ S36(2)*y(22)+ 0+ 0 0+ 0+ S37(3)*y(3)+ 0+ 0+ S38(2)*y(6)+ S38(3)*y(7)+ S38(4)*y(8)+ 0+ 0+ S39(3)*y(11)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S41(3)*y(19)+ 0+ 0+ 0+ S42(3)*y(23)+ 0 0+ 0+ 0+ S43(4)*y(4)+ 0+ 0+ S44(3)*y(7)+ S44(4)*y(8)+ 0+ 0+ 0+ S45(4)*y(12)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S47(4)*y(20)+ 0+ 0+ 0+ S48(4)*y(24) 0+ 0+ 0+ 0+ S50(1)*y(5)+ 0+ 0+ 0+ S51(1)*y(9)+ S51(2)*y(10)+ 0+ 0+ S52(1)*y(13)+ 0+ 0+ 0+ S53(1)*y(17)+ 0+ 0+ 0+ S54(1)*y(21)+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ S56(2)*y(6)+ 0+ 0+ S57(1)*y(9)+ S57(2)*y(10)+ S57(3)*y(11)+ 0+ 0+ S58(2)*y(14)+ 0+ 0+ 0+ S59(2)*y(18)+ 0+ 0+ 0+ S60(2)*y(22)+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ S62(3)*y(7)+ 0+ 0+ S63(2)*y(10)+ S63(3)*y(11)+ S63(4)*y(12)+ 0+ 0+ S64(3)*y(15)+ 0+ 0+ 0+ S65(3)*y(19)+ 0+ 0+ 0+ S66(3)*y(23)+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ S68(4)*y(8)+ 0+ 0+ S69(3)*y(11)+ S69(4)*y(12)+ 0+ 0+ 0+ S70(4)*y(16)+ 0+ 0+ 0+ S71(4)*y(20)+ 0+ 0+ 0+ S72(4)*y(24) S73(1)*y(1)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S75(1)*y(9)+ 0+ 0+ 0+ S76(1)*y(13)+ S76(2)*y(14)+ 0+ 0+ S77(1)*y(17)+ 0+ 0+ 0+ S78(1)*y(21)+ 0+ 0+ 0 0+ S79(2)*y(2)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S81(2)*y(10)+ 0+ 0+ S82(1)*y(13)+ S82(2)*y(14)+ S82(3)*y(15)+ 0+ 0+ S83(2)*y(18)+ 0+ 0+ 0+ S84(2)*y(22)+ 0+ 0 0+ 0+ S85(3)*y(3)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S87(3)*y(11)+ 0+ 0+ S88(2)*y(14)+ S88(3)*y(15)+ S88(4)*y(16)+ 0+ 0+ S89(3)*y(19)+ 0+ 0+ 0+ S90(3)*y(23)+ 0 0+ 0+ 0+ S91(4)*y(4)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S93(4)*y(12)+ 0+ 0+ S94(3)*y(15)+ S94(4)*y(16)+ 0+ 0+ 0+ S95(4)*y(20)+ 0+ 0+ 0+ S96(4)*y(24) S97(1)*y(1)+ 0+ 0+ 0+ S98(1)*y(5)+ 0+ 0+ 0+ S99(1)*y(9)+ 0+ 0+ 0+ S100(1)*y(13)+ 0+ 0+ 0+ S101(1)*y(17)+ S101(2)*y(18)+ 0+ 0+ 0+ 0+ 0+ 0 0+ S103(2)*y(2)+ 0+ 0+ 0+ S104(2)*y(6)+ 0+ 0+ 0+ S105(2)*y(10)+ 0+ 0+ 0+ S106(2)*y(14)+ 0+ 0+ S107(1)*y(17)+ S107(2)*y(18)+ S107(3)*y(19)+ 0+ 0+ 0+ 0+ 0 0+ 0+ S109(3)*y(3)+ 0+ 0+ 0+ S110(3)*y(7)+ 0+ 0+ 0+ S111(3)*y(11)+ 0+ 0+ 0+ S112(3)*y(15)+ 0+ 0+ S113(2)*y(18)+ S113(3)*y(19)+ S113(4)*y(20)+ 0+ 0+ 0+ 0 0+ 0+ 0+ S115(4)*y(4)+ 0+ 0+ 0+ S116(4)*y(8)+ 0+ 0+ 0+ S117(4)*y(12)+ 0+ 0+ 0+ S118(4)*y(16)+ 0+ 0+ S119(3)*y(19)+ S119(4)*y(20)+ 0+ 0+ 0+ 0 S121(1)*y(1)+ 0+ 0+ 0+ S122(1)*y(5)+ 0+ 0+ 0+ S123(1)*y(9)+ 0+ 0+ 0+ S124(1)*y(13)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ S126(1)*y(21)+ S126(2)*y(22)+ 0+ 0 0+ S127(2)*y(2)+ 0+ 0+ 0+ S128(2)*y(6)+ 0+ 0+ 0+ S129(2)*y(10)+ 0+ 0+ 0+ S130(2)*y(14)+ 0+ 0+ 0+ 0+ 0+ 0+ S132(1)*y(21)+ S132(2)*y(22)+ S132(3)*y(23)+ 0 0+ 0+ S133(3)*y(3)+ 0+ 0+ 0+ S134(3)*y(7)+ 0+ 0+ 0+ S135(3)*y(11)+ 0+ 0+ 0+ S136(3)*y(15)+ 0+ 0+ 0+ 0+ 0+ 0+ S138(2)*y(22)+ S138(3)*y(23)+ S138(4)*y(24) 0+ 0+ 0+ S139(4)*y(4)+ 0+ 0+ 0+ S140(4)*y(8)+ 0+ 0+ 0+ S141(4)*y(12)+ 0+ 0+ 0+ S142(4)*y(16)+ 0+ 0+ 0+ 0+ 0+ 0+ S144(3)*y(23)+ S144(4)*y(24) ]; StT = [ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 ]; StN = [ S_1*y(44)+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ S_2*y(45)+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ S_3*y(46)+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ S_4*y(47)+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ S_5*y(48)+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ S_6*y(49)+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 ]; TqS = [ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 ]; TqA = [ T1(1)*y(25)+ T1(2)*y(26)+ 0+ T2(1)*y(28)+ 0+ 0+ 0+ 0+ 0+ T4(1)*y(34)+ 0+ 0+ T5(1)*y(37)+ 0+ 0+ T6(1)*y(40)+ 0+ 0+ 0 T7(1)*y(25)+ T7(2)*y(26)+ T7(3)*y(27)+ 0+ T8(2)*y(29)+ 0+ 0+ 0+ 0+ 0+ T10(2)*y(35)+ 0+ 0+ T11(2)*y(38)+ 0+ 0+ T12(2)*y(41)+ 0+ 0 0+ T13(2)*y(26)+ T13(3)*y(27)+ 0+ 0+ T14(3)*y(30)+ 0+ 0+ 0+ 0+ 0+ T16(3)*y(36)+ 0+ 0+ T17(3)*y(39)+ 0+ 0+ T18(3)*y(42)+ Cen_tqlo_S1*y(43) T19(1)*y(25)+ 0+ 0+ T20(1)*y(28)+ T20(2)*y(29)+ 0+ T21(1)*y(31)+ 0+ 0+ 0+ 0+ 0+ T23(1)*y(37)+ 0+ 0+ T24(1)*y(40)+ 0+ 0+ 0 0+ T25(2)*y(26)+ 0+ T26(1)*y(28)+ T26(2)*y(29)+ T26(3)*y(30)+ 0+ T27(2)*y(32)+ 0+ 0+ 0+ 0+ 0+ T29(2)*y(38)+ 0+ 0+ T30(2)*y(41)+ 0+ 0 0+ 0+ T31(3)*y(27)+ 0+ T32(2)*y(29)+ T32(3)*y(30)+ 0+ 0+ T33(3)*y(33)+ 0+ 0+ 0+ 0+ 0+ T35(3)*y(39)+ 0+ 0+ T36(3)*y(42)+ Cen_tqlo_S2*y(43) 0+ 0+ 0+ T38(1)*y(28)+ 0+ 0+ T39(1)*y(31)+ T39(2)*y(32)+ 0+ T40(1)*y(34)+ 0+ 0+ T41(1)*y(37)+ 0+ 0+ T42(1)*y(40)+ 0+ 0+ 0 0+ 0+ 0+ 0+ T44(2)*y(29)+ 0+ T45(1)*y(31)+ T45(2)*y(32)+ T45(3)*y(33)+ 0+ T46(2)*y(35)+ 0+ 0+ T47(2)*y(38)+ 0+ 0+ T48(2)*y(41)+ 0+ 0 0+ 0+ 0+ 0+ 0+ T50(3)*y(30)+ 0+ T51(2)*y(32)+ T51(3)*y(33)+ 0+ 0+ T52(3)*y(36)+ 0+ 0+ T53(3)*y(39)+ 0+ 0+ T54(3)*y(42)+ Cen_tqlo_S3*y(43) T55(1)*y(25)+ 0+ 0+ 0+ 0+ 0+ T57(1)*y(31)+ 0+ 0+ T58(1)*y(34)+ T58(2)*y(35)+ 0+ T59(1)*y(37)+ 0+ 0+ T60(1)*y(40)+ 0+ 0+ 0 0+ T61(2)*y(26)+ 0+ 0+ 0+ 0+ 0+ T63(2)*y(32)+ 0+ T64(1)*y(34)+ T64(2)*y(35)+ T64(3)*y(36)+ 0+ T65(2)*y(38)+ 0+ 0+ T66(2)*y(41)+ 0+ 0 0+ 0+ T67(3)*y(27)+ 0+ 0+ 0+ 0+ 0+ T69(3)*y(33)+ 0+ T70(2)*y(35)+ T70(3)*y(36)+ 0+ 0+ T71(3)*y(39)+ 0+ 0+ T72(3)*y(42)+ Cen_tqlo_S4*y(43) T73(1)*y(25)+ 0+ 0+ T74(1)*y(28)+ 0+ 0+ T75(1)*y(31)+ 0+ 0+ T76(1)*y(34)+ 0+ 0+ T77(1)*y(37)+ T77(2)*y(38)+ 0+ 0+ 0+ 0+ 0 0+ T79(2)*y(26)+ 0+ 0+ T80(2)*y(29)+ 0+ 0+ T81(2)*y(32)+ 0+ 0+ T82(2)*y(35)+ 0+ T83(1)*y(37)+ T83(2)*y(38)+ T83(3)*y(39)+ 0+ 0+ 0+ 0 0+ 0+ T85(3)*y(27)+ 0+ 0+ T86(3)*y(30)+ 0+ 0+ T87(3)*y(33)+ 0+ 0+ T88(3)*y(36)+ 0+ T89(2)*y(38)+ T89(3)*y(39)+ 0+ 0+ 0 + Cen_tqlo_S5*y(43) T91(1)*y(25)+ 0+ 0+ T92(1)*y(28)+ 0+ 0+ T93(1)*y(31)+ 0+ 0+ T94(1)*y(34)+ 0+ 0+ 0+ 0+ 0+ T96(1)*y(40)+ T96(2)*y(41)+ 0+ 0 0+ T97(2)*y(26)+ 0+ 0+ T98(2)*y(29)+ 0+ 0+ T99(2)*y(32)+ 0+ 0+ T100(2)*y(35)+ 0+ 0+ 0+ 0+ T102(1)*y(40)+ T102(2)*y(41)+ T102(3)*y(42)+ 0 0+ 0+ T103(3)*y(27)+ 0+ 0+ T104(3)*y(30)+ 0+ 0+ T105(3)*y(33)+ 0+ 0+ T106(3)*y(36)+ 0+ 0+ 0+ 0+ T108(2)*y(41)+ T108(3)*y(42)+ Cen_tqlo_S6*y(43) 0+ 0+ Ctq_1*y(27)+ 0+ 0+ Ctq_2*y(30)+ 0+ 0+ Ctq_3*y(33)+ 0+ 0+ Ctq_4*y(36)+ 0+ 0+ Ctq_5*y(39)+ 0+ 0+ Ctq_6*y(42)+ Ctq_7*y(43) ]; TqN = [ T_1*y(44)+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ T_2*y(45)+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ T_3*y(46)+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ T_4*y(47)+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ T_5*y(48)+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ T_6*y(49)+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0 ]; NaS = [ NG11*y(1)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ NG21*y(5)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG31*y(9)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG41*y(13)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG51*y(17)+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG61*y(21)+ 0+ 0+ 0 ]; NaT = [ NG12*y(25)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ NG22*y(28)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ NG32*y(31)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG42*y(34)+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG52*y(37)+ 0+ 0+ 0+ 0+ 0+ 0 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ 0+ NG62*y(40)+ 0+ 0+ 0 ]; NaA = [ NG13*y(44)+ 0+ 0+ 0+ 0+ 0+ 60 0+ NG23*y(45)+ 0+ 0+ 0+ 0+ 60 0+ 0+ NG33*y(46)+ 0+ 0+ 0+ 60 0+ 0+ 0+ NG43*y(47)+ 0+ 0+ 60 0+ 0+ 0+ 0+ NG53*y(48)+ 0+ 60 0+ 0+ 0+ 0+ 0+ NG63*y(49)+ 60 ]; yp = [ StA + StT + StN TqS + TqA + TqN NaS + NaT + NaA ];