Error using call_myequations Too many input arguments.

Hi I divided my script into 4 script:

parameters.m

function k= parameters
     %constantes :::
    %temperature
    T=17000;
    %constante k: recombinaison radiative
    k(1)=2.36 *(1e-12)*(T/300)^(-0.29)*exp(17.60/T); 
    k(2) =3.24 *(1e-12) *(T/300)^(-0.66);
    k(3)=3.50 *(1e-12)*(T/300)^(-0.53)* exp(3.20/T);
    %constante k::Reconmbinaison 
    k(4)=1*10^(-10); 
    k(5)=4.54*10^(-10);
    k(6)=9.10*10^(-10);
    k(7)=8.98*10^(-9)*(T/300)^(-0.5);
    k(8)=5.56*10^(-11)*(T/300)^(0.41)*exp(26.90/T);
    k(9)=4.98*10^(-10)*exp(-18116/T);
    k(10)=6*10^(-11)*(T/300)^(-0.16);
    k(11)=8.69*10^(-11)*exp(-22600/T);
    k(12)=1*10^(-10);
    k(13)=2.94*10^(-11)*(T/300)^(0.5)*exp(-58025/T);
    k(14)=1*10^(-9);
    k(15)=4.80*10^(-10);
    k(16)=1.00*10^(-9)*(T/300)^(-0.5);
    k(17)=8.30*10^(-10);
    k(18)=2.42*10^(-12)*(T/300)^(-0.21)*exp(44/T);
    k(19)=1.18*10^(-11)*exp(-20413/T);
    k(20)=9.82*10^(-12)*(T/300)^(-0.21)*exp(-5.20/T);
    k(21)=1.66*10^(-10)*exp(-14100/T);
    k(22)=1.15*10^(-10)*exp(-13400/T);
    k(23)=2.51*10^(-10)*exp(-38602/T);
    k(24)=1.00*10^(-10);
    k(25)=5.37*10^(-11)*exp(-13800/T);
    k(26)=5.00*10^(-11)*exp(-200/T);
    k(27)=5.00*10^(-12)*exp(-900/T);
    k(28)=6.00*10^(-12);
    k(29)=2.00*10^(-10)*(T/300)^(-0.12);
    k(30)=1.30*10^(-10);
    k(31)=5.00*10^(-10);
    k(32)=2.90*10^(-10);
    k(33)=3.10*10^(-10);
    k(34)=3.66*10^(-11);
    k(35)=2.26*10^(-12)*(T/300)^(0.86)*exp(-3134.0/T);
    k(36)=1.80*10^(-10);
    k(37)=3.38*10^(-11)*(T/300)^(-0.17)*exp(2.80/T);
    k(38)=3*10^(-12);
    k(39)=1.00*10^(-10)*(T/300)^(0.4);
    k(40)=6.10*10^(-10);
    k(41)=1*10^(-13);
    k(42)=5*10^(-11);
    k(43)=5*10^(-11);
    k(44)=1.15*10^(-10);
    k(45)=4*10^(-11);
    k(46)=1*10^(-10);
    k(47)=1.50*10^(-11)*exp(-4300/T);
    k(48)=5.3137*10^(-10);
    k(49)=4.10*10^(-10);
    k(50)=8*10^(-10);
    k(51)=1.32*10^(-12);
    k(52)=2.8*10^(-12)*exp(23400/T);
    k(53)=1*10^(-10)*exp(-55200/T);
    k(54)=8.60*10^(-11);
    k(55)=5.99*10^(-12)*exp(-24075/T);
    k(56)=2.02*10^(-11)*(T/300)^(-0.19)*exp(31.90/T);
    k(57)=2.51*10^(-11)*exp(-30653/T);
    k(58)=3.47*10^(-11)*(T/300)^(-1.33)*exp(-242/T);
    
    %constante k::Recombinaison dissociative
    k(59)=3*10^(-7)*(T/300)^(-0.5);
    k(60)=1.70*10^(-7)*(T/300)^(-0.3);
     
    %constante k::Association radiative
    k(61)=4.01*10^(-18)*(T/300)^(0.17)*exp(-101.50/T);
    k(62)=1.08*10^(-18)*(T/300)^(0.07)*exp(-57.50/T);
    k(63)=3.14*10^(-18)*(T/300)^(-0.15)*exp(-68/T);
    k(64)=4.00*10^(-14)*(T/300)^(-1);
    k(65)=4.69*10^(-19)*(T/300)^(1.52)*exp(50.50/T);
    k(66)=3*10^(-16)*(T/300)^(-1);
    k(67)=5*10^(-10)*(T/300)^(-3.7)*exp(-800/T);
    k(68)=4.36*10^(-18)*(T/300)^(0.35)*exp(-161.30/T);
    k(69)=5.72*10^(-19)*(T/300)^(0.37)*exp(-51/T);
    k(70)=4.90*10^(-20)*(T/300)^(1.58);
    k(71)=3.71*10^(-18)*(T/300)^(0.24)*exp(-26.10/T);
    
    %constante k::Transfert de charge
    k(72)=7.05*10^(-10)*(T/300)^(-0.03)*exp(16.70/T);
    k(73)=1*10^(-10);
    k(74)=1.10*10^(-10);
    k(75)=1.10*10^(-10);
    k(76)=1.10*10^(-10);
    k(77)=5.20*10^(-11);
    k(78)=4.80*10^(-10);
    k(79)=4.90*10^(-12)*(T/300)^(0.5)*exp(-4580/T);
    k(80)=6.30*10^(-10);
    k(81)=1.90*10^(-11);
    k(82)=1.00*10^(-11);
    k(83)=1.40*10^(-10);
    k(84)=7.30*10^(-10)*exp(-890/T);
    k(85)=1.00*10^(-9);
    k(86)=3.11*10^(-10);
    k(87)=1.00*10^(-9);
    k(88)=8.25*10^(-10);
    k(89)=1.10*10^(-9)*(T/300)^(-0.50);
    k(90)=4.51*10^(-10);
    k(91)=1* 10^(-11);
    k(92)=8.40*10^(-10);
    k(93)=8.50^(-10);
    k(94)=3.40*10^(-10);
    k(95)=5*10^(-11);
    k(96)=1.20*10^(-10);
    k(97)=6.60*10^(-10);
    k(98)=4.60*10^(-10);
    k(99)=1.00*10^(-10)*(T/300)^(-0.5);
    k(100)=7.40*10^(-11);
    k(101)=6.30*10^(-10);
    k(102)=3.30*10^(-10);
    k(103)=5.70*10^(-10);
    
    %constante k::Neutralisation
    k(104)=7.51*10^(-8)*(T/300)^(-0.5);
    k(105)=7.51*10^(-8)*(T/300)^(-0.5);
    k(106)=7.51*10^(-8)*(T/300)^(-0.5);
    k(107)=7.51*10^(-8)*(T/300)^(-0.5);
    k(108)=7.51*10^(-8)*(T/300)^(-0.5);
    k(109)=7.51*10^(-8)*(T/300)^(-0.5);
    k(110)=7.51*10^(-8)*(T/300)^(-0.50);
    k(111)=7.51*10^(-8)*(T/300)^(-0.50);
    k(112)=7.51*10^(-8)*(T/300)^(-0.50);
    k(113)=7.51*10^(-8)*(T/300)^(-0.50);
    k(114)=7.51*10^(-8)*(T/300)^(-0.5);
    k(115)=7.51*10^(-8)*(T/300)^(-0.5);
    k(116)=7.51*10^(-8)*(T/300)^(-0.5);
    k(117)=7.51*10^(-8)*(T/300)^(-0.5);
    %constante k::Photoionisation
    k(118)=3.10*10^(-10)*exp(-3.3/T);
    k(119)=4.90*10^(-8)*exp(-0.50/T);
    k(120)=1.09*10^(-8)*exp(-0.5/T);
    k(121)=4.10*10^(-10)*exp(-3.80/T);
    k(122)=1*10^(-11)*exp(-1.70/T);
    k(123)=2.96*10^(-9)*exp(-2/T);
    k(124)=5*10^(-9)*exp(-2.1/T);
    k(125)=6.88*10^(-9)*exp(-1.5/T);
    k(126)=7.90*10^(-10)*exp(-2.1/T);
    k(127)=6.10*10^(-9)*exp(-0.5/T);
    k(128)=3.50*10^(-11)*exp(-2/T);
    k(129)=2.30*10^(-10)*exp(-3.90/T);
    k(130)=2.10*10^(-10)*exp(-3.50/T);
    k(131)=4.70*10^(-10)*exp(-2.10/T);
    %constante k::Attachement radiatif
    
    k(132)=2.25*10^(-15);
    k(133)=1.50*10^(-15);
    k(134)=2.00*10^(-15)*(T/300)^(-0.5);
    k(135)=1.70*10^(-14)*(T/300)^(-0.5);
    %constante k::Détachement associatif
    
    k(136)=5*10^(-10);
    k(137)=1*10^(-9);
    k(138)=1*10^(-9);
    k(139)=5*10^(-10);
    k(140)=5*10^(-10);
    k(141)=5*10^(-11);
    k(142)=5*10^(-10);
    k(143)=1.90*10^(-10);
    k(144)=2.90*10^(-10);
    k(145)=3.10*10^(-10)*(T/300)^(-0.83);
    k(146)=6.50*10^(-10);
    k(147)=2.20*10^(-10);
    k(148)=1.00*10^(-9);
    k(149)=1*10^(-9);
    
    %constante k::Excitation/ Désexcitation
    k(150)=2.0425*10^(-7)*T^(-0.3)*exp(-29489.55/T);
    k(151)=2.6*10^(-10)*T^(0.89)*exp(-11372.9/T);
    k(152)=1.2*10^(-10)*T^(0.59)*exp(-18916.15/T);
    k(153)=6.0278*10^(-7)*T^(-0.4)*exp(-19880.36/T);
    k(154)=1.6306*10^(-17);
    k(155)=3.6532*10^(-11);
    
    %constante k::Détachement électronique
    k(156)=7*10^(-10);
    k(157)=2*10^(-10);
    k(158)=3.6*10^(-10);
    
    %constante k::Dissociation par impact électronique
    k(159)= 6.2*10^(-10) * T^(1.75) * exp(-113032.7/T );
    
    %constante k::Ionisation par impact électronique
    k(160)= 7.1*10^(-9) * T^(0.76) * exp(-181038/T );
    
    %constante k::Attachement électronique cn(1)6/s
    k(161)=1*10^(-31);
    k(162)=1*10^(-31);
    
    %constante k::Processus de recombinaison cn(1)6/s
    k(163)=5.70*10^(-34)*T^(-2.6);
    k(164)=2.5919*10^(-80);
    k(165)=1.31*10^(-31)*T^(-1.5);
end 
    
    

myequations.m

 function dndt = myequations(t,n,k)
    dndt = zeros(43,1);
    %Chemin C
    dndt(1)= k(1)*n(2)*n(42)+k(15)*n(5)*n(9)+k(16)*n(5)*n(26)+k(25)*n(4)*n(26)...
    +k(29)*n(4)*n(9)+k(33)*n(4)*n(11)+k(39)*n(7)*n(26)+k(40)*n(7)*n(27)...
    +k(43)*n(7)*n(9)+k(44)*n(7)*n(12)+k(48)*n(9)*n(9)+k(59)*n(11)*n(42)...
    +k(72)*n(2)*n(24)+k(104)*n(3)*n(5)+k(105)*n(3)*n(8)+k(106)*n(3)*n(25)...
    +k(117)*n(18)*n(2)+k(119)*n(3)*n(43)+k(122)*n(11)*n(43)+k(124)*n(13)*n(43)...
    +k(130)*n(22)*n(43)-k(4)*n(1)*n(37)-k(8)*n(1)*n(15)-k(9)*n(1)*n(26)...
    -k(10)*n(1)*n(24)-k(11)*n(1)*n(20)-k(12)*n(1)*n(35)-k(13)*n(1)*n(22)...
    -k(61)*n(2)*n(1)-k(64)*n(1)*n(13)-k(65)*n(1)*n(4)-k(66)*n(1)*n(9)...
    -k(67)*n(1)*n(5)-k(68)*n(1)*n(1)-k(69)*n(1)*n(7)-k(73)*n(1)*n(27)...
    -k(74)*n(1)*n(11)-k(75)*n(1)*n(23)-k(76)*n(1)*n(21)-k(77)*n(1)*n(19)...
    -k(118)*n(1)*n(43)-k(132)*n(1)*n(42)-k(136)*n(1)*n(3)-k(137)*n(1)*n(12)...
    -k(138)*n(1)*n(14)-k(139)*n(1)*n(6)-2*k(164)*n(7)*n(1)*n(1)...
    +k(164)*n(7)*n(1)^(2)-k(163)*n(1)*n(15)*n(4)+k(163)*n(1)*n(15)*n(4)...
    -k(165)*n(1)*n(32)*n(4)+k(165)*n(1)*n(32)*n(4);
    %Chemin C+:
    dndt(2)= -k(1)*n(2)*n(42)-k(5)*n(2)*n(15)-k(6)*n(2)*n(30)-k(7)*n(2)*n(36)...
    -k(61)*n(2)*n(1)-k(62)*n(2)*n(7)-k(63)*n(2)*n(4)-k(72)*n(2)*n(24)...
    -k(117)*n(18)*n(2)+k(45)*n(7)*n(11)+k(73)*n(1)*n(27)...
    +k(74)*n(1)*n(11)+k(75)*n(1)*n(23)+k(76)*n(1)*n(21)...
    +k(77)*n(1)*n(19)+k(118)*n(1)*n(43)+k(122)*n(11)*n(43);
    %Chemin C-:
    dndt(3)= k(32)*n(4)*n(12)+k(132)*n(1)*n(42)-k(14)*n(3)*n(24)-k(104)*n(3)*n(5)...
    -k(105)*n(3)*n(8)-k(106)*n(3)*n(25)-k(119)*n(3)*n(43)...
    -k(136)*n(1)*n(3)-k(140)*n(3)*n(4)-k(141)*n(3)*n(15)-k(142)*n(3)*n(7);
    %Chemin O:
    dndt(4)= k(2)*n(5)*n(42)+k(8)*n(1)*n(15)+k(10)*n(1)*n(24)+k(13)*n(1)*n(22)...
    +k(14)*n(3)*n(24)+k(35)*n(7)*n(15)+k(36)*n(7)*n(19)...
    +k(37)*n(7)*n(24)+k(38)*n(7)*n(32)+k(52)*n(15)*n(24)...
    +k(53)*n(15)*n(20)+k(55)*n(15)*n(22)+k(56)*n(15)*n(26)...
    +k(78)*n(5)*n(9)+k(79)*n(5)*n(22)+k(80)*n(5)*n(30)...
    +k(81)*n(5)*n(15)+k(84)*n(6)*n(15)+k(85)*n(6)*n(26)...
    +k(104)*n(3)*n(5)+2*k(107)*n(5)*n(6)+k(108)*n(5)*n(18)...
    +k(109)*n(6)*n(8)+k(120)*n(6)*n(43)+2*k(126)*n(15)*n(43)...
    +k(128)*n(19)*n(43)+k(130)*n(22)*n(43)+k(131)*n(24)*n(43)...
    +k(162)*n(42)*n(4)*n(15)-k(19)*n(4)*n(24)-k(20)*n(4)*n(32)...
    -k(21)*n(4)*n(30)-k(22)*n(4)*n(30)-k(23)*n(4)*n(20)...
    -k(24)*n(4)*n(36)-k(25)*n(4)*n(26)-k(26)*n(4)*n(26)...
    -k(27)*n(4)*n(13)-k(28)*n(4)*n(35)-k(29)*n(4)*n(9)...
    -k(30)*n(4)*n(21)-k(31)*n(4)*n(14)-k(32)*n(4)*n(12)...
    -k(33)*n(4)*n(11)-k(63)*n(2)*n(4)-k(65)*n(1)*n(4)...
    -2*k(70)*n(4)*n(4)-k(82)*n(4)*n(21)-k(83)*n(4)*n(23)...
    -k(133)*n(4)*n(42)-k(140)*n(3)*n(4)-k(143)*n(4)*n(6)...
    -k(144)*n(4)*n(12)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15)-k(163)*n(1)*n(15)*n(4)-k(165)*n(1)*n(32)*n(4);
    %Chemin O+:
    dndt(5)= +k(5)*n(2)*n(15)+k(34)*n(8)*n(15)+k(82)*n(4)*n(21)...
    +k(83)*n(4)*n(23)+k(128)*n(19)*n(43)-k(2)*n(5)*n(42)-k(15)*n(5)*n(9)...
    -k(16)*n(5)*n(26)-k(17)*n(5)*n(32)-k(18)*n(5)*n(20)-k(67)*n(1)*n(5)...
    -k(78)*n(5)*n(9)-k(79)*n(5)*n(22)-k(80)*n(5)*n(30)-k(81)*n(5)*n(15)...
    -k(104)*n(3)*n(5)-k(107)*n(5)*n(6)-k(108)*n(5)*n(18);
    %Chemin O-:
    dndt(6)= k(133)*n(4)*n(42)+k(161)*n(42)*n(4)*n(15)-k(84)*n(6)*n(15)...
    -k(85)*n(6)*n(26)-k(107)*n(5)*n(6)-k(109)*n(6)*n(8)-k(120)*n(6)*n(43)...
    -k(139)*n(1)*n(6)-k(143)*n(4)*n(6)-k(145)*n(6)*n(24)...
    -k(146)*n(6)*n(22)-k(147)*n(6)*n(7);
    %Chemin N:
    dndt(7)= k(3)*n(8)*n(42)+k(9)*n(1)*n(26)+k(11)*n(1)*n(20)...
    +k(18)*n(5)*n(20)+k(19)*n(4)*n(24)+k(23)*n(4)*n(20)...
    +k(26)*n(4)*n(26)+k(30)*n(4)*n(21)+2*k(60)*n(21)*n(42)...
    +k(86)*n(8)*n(15)+k(87)*n(8)*n(9)+k(88)*n(8)*n(22)...
    +k(89)*n(8)*n(26)+k(90)*n(8)*n(24)+k(105)*n(3)*n(8)...
    +k(109)*n(6)*n(8)+k(110)*n(8)*n(18)+k(111)*n(8)*n(12)...
    +k(112)*n(8)*n(28)+k(113)*n(8)*n(14)+2*k(129)*n(20)*n(43)...
    +k(131)*n(24)*n(43)+2*k(159)*n(20)*n(42)-k(35)*n(7)*n(15)...
    -k(36)*n(7)*n(19)-k(37)*n(7)*n(24)-k(38)*n(7)*n(32)...
    -k(39)*n(7)*n(26)-k(40)*n(7)*n(27)-k(41)*n(7)*n(13)...
    -k(42)*n(7)*n(14)-k(43)*n(7)*n(9)-k(44)*n(7)*n(12)...
    -k(45)*n(7)*n(11)-k(46)*n(7)*n(35)-k(62)*n(2)*n(7)...
    -k(69)*n(1)*n(7)-k(71)*n(8)*n(7)-k(91)*n(7)*n(21)...
    -k(142)*n(3)*n(7)-k(147)*n(6)*n(7)-k(164)*n(7)*n(1)*n(1);
    %Chemin N+:
    dndt(8)=k(91)*n(7)*n(21)-k(3)*n(8)*n(42)-k(34)*n(8)*n(15)...
    -k(71)*n(8)*n(7)-k(86)*n(8)*n(15)-k(87)*n(8)*n(9)...
    -k(88)*n(8)*n(22)-k(89)*n(8)*n(26)-k(90)*n(8)*n(24)...
    -k(105)*n(3)*n(8)-k(109)*n(6)*n(8)-k(110)*n(8)*n(18)...
    -k(111)*n(8)*n(12)-k(112)*n(8)*n(28)-k(113)*n(8)*n(14);
    %Chemin C2:
    dndt(9)= k(9)*n(1)*n(26)+k(12)*n(1)*n(35)+k(13)*n(1)*n(22)...
    +k(27)*n(4)*n(13)+k(41)*n(7)*n(13)+k(68)*n(1)*n(1)...
    +k(74)*n(1)*n(11)+k(94)*n(11)*n(24)+k(111)*n(8)*n(12)...
    +k(115)*n(12)*n(25)+k(121)*n(9)*n(43)+k(123)*n(12)*n(43)...
    +k(124)*n(13)*n(43)+k(136)*n(1)*n(3)+k(153)*n(10)*n(42)...
    +k(154)*n(10)^(2)+k(155)*n(10)*n(1)-k(15)*n(5)*n(9)...
    -k(29)*n(4)*n(9)-k(43)*n(7)*n(9)-k(47)*n(9)*n(15)...
    -2*k(48)*n(9)*n(9)-k(49)*n(9)*n(19)-k(66)*n(1)*n(9)...
    -k(78)*n(5)*n(9)-k(87)*n(8)*n(9)-k(92)*n(9)*n(23)...
    -k(93)*n(9)*n(27)-k(121)*n(9)*n(43)-k(134)*n(9)*n(42)...
    -k(148)*n(9)*n(12)-k(150)*n(9)*n(42);
    %Chemin C2(d)
    dndt(10)= k(150)*n(9)*n(42)-k(153)*n(10)*n(42)-k(154)*n(10)^(2)-k(155)*n(10)*n(1);
    %Chemin C2+:
    dndt(11)= k(61)*n(2)*n(1)+k(78)*n(5)*n(9)+k(87)*n(8)*n(9)...
    +k(92)*n(9)*n(23)+k(93)*n(9)*n(27)-k(33)*n(4)*n(11)...
    -k(45)*n(7)*n(11)-k(50)*n(11)*n(15)-k(59)*n(11)*n(42)...
    -k(74)*n(1)*n(11)-k(94)*n(11)*n(24)-k(122)*n(11)*n(43);
    %Chemin C2-:
    dndt(12)= -k(32)*n(4)*n(12)-k(44)*n(7)*n(12)-k(111)*n(8)*n(12)...
    -k(115)*n(12)*n(25)-k(123)*n(12)*n(43)-k(137)*n(1)*n(12)...
    -k(144)*n(4)*n(12)-k(148)*n(9)*n(12)-k(149)*n(12)*n(13)...
    +k(31)*n(4)*n(14)+k(42)*n(7)*n(14)+k(134)*n(9)*n(42);
    %Chemin C3:
    dndt(13)= k(48)*n(9)*n(9)+k(66)*n(1)*n(9)+k(113)*n(8)*n(14)...
    +k(114)*n(14)*n(25)+k(125)*n(14)*n(43)...
    +k(137)*n(1)*n(12)-k(27)*n(4)*n(13)-k(41)*n(7)*n(13)...
    -k(64)*n(1)*n(13)-k(124)*n(13)*n(43)...
    -k(135)*n(13)*n(42)-k(149)*n(12)*n(13);
    %Chemin C3-:
    dndt(14)= k(135)*n(13)*n(42)-k(31)*n(4)*n(14)...
    -k(42)*n(7)*n(14)-k(113)*n(8)*n(14)...
    -k(114)*n(14)*n(25)-k(125)*n(14)*n(43)-k(138)*n(1)*n(14);
    %Chemin O2:
    dndt(15)= k(17)*n(5)*n(32)+k(19)*n(4)*n(24)+k(20)*n(4)*n(32)...
    +k(21)*n(4)*n(30)+k(57)*n(24)*n(24)+k(70)*n(4)*n(4)...
    +k(77)*n(1)*n(19)+k(97)*n(19)*n(32)+k(98)*n(19)*n(24)...
    +k(108)*n(5)*n(18)+k(110)*n(8)*n(18)+k(116)*n(18)*n(25)...
    +k(117)*n(18)*n(2)+k(127)*n(18)*n(43)+k(143)*n(4)*n(6)...
    +2*k(156)*n(18)*n(15)+2*k(157)*n(18)*n(16)...
    +2*k(158)*n(18)*n(17)+k(161)*n(42)*n(4)*n(15)...
    -k(5)*n(2)*n(15)-k(8)*n(1)*n(15)-k(34)*n(8)*n(15)...
    -k(35)*n(7)*n(15)-k(47)*n(9)*n(15)-k(50)*n(11)*n(15)...
    -k(51)*n(15)*n(37)-k(52)*n(15)*n(24)-k(53)*n(15)*n(20)...
    -k(54)*n(15)*n(27)-k(55)*n(15)*n(22)-k(56)*n(15)*n(26)...
    -k(81)*n(5)*n(15)-k(84)*n(6)*n(15)-k(86)*n(8)*n(15)...
    -k(95)*n(15)*n(21)-k(96)*n(15)*n(23)-k(126)*n(15)*n(43)...
    -k(141)*n(3)*n(15)-k(151)*n(15)*n(42)-k(152)*n(15)*n(42)...
    -k(156)*n(18)*n(15)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15)-k(163)*n(1)*n(15)*n(4);
    %Chemin O2(a):
    dndt(16)= k(151)*n(42)*n(15)-k(157)*n(16)*n(18);
    %Chemin O2(b):
    dndt(17)= k(152)*n(42)*n(15)-k(158)*n(17)*n(18);
    %Chemin O2-:
    dndt(18)= k(84)*n(6)*n(15)+k(162)*n(42)*n(4)*n(15)...
    -k(108)*n(5)*n(18)-k(110)*n(8)*n(18)...
    -k(116)*n(18)*n(25)-k(117)*n(18)*n(2)-k(127)*n(18)*n(43)...
    -k(156)*n(18)*n(15)-k(157)*n(18)*n(16)-k(158)*n(18)*n(17);
    %Chemin O2+:
    dndt(19)= k(81)*n(5)*n(15)+k(86)*n(8)*n(15)+k(95)*n(15)*n(21)...
    +k(96)*n(15)*n(23)-k(36)*n(7)*n(19)-k(49)*n(9)*n(19)-k(77)*n(1)*n(19)...
    -k(97)*n(19)*n(32)-k(98)*n(19)*n(24)-k(128)*n(19)*n(43);
    %Chemin N2:
    dndt(20)= k(21)*n(4)*n(30)+k(37)*n(7)*n(24)+k(39)*n(7)*n(26)...
    +k(57)*n(24)*n(24)+k(76)*n(1)*n(21)+k(82)*n(4)*n(21)...
    +k(91)*n(7)*n(21)+k(95)*n(15)*n(21)+k(99)*n(21)*n(26)...
    +k(100)*n(21)*n(22)-k(11)*n(1)*n(20)-k(18)*n(5)*n(20)...
    -k(23)*n(4)*n(20)-k(53)*n(15)*n(20)...
    -k(129)*n(20)*n(43)-k(159)*n(20)*n(42)-k(160)*n(20)*n(42);
    %Chemin N2+:
    dndt(21)= k(40)*n(7)*n(27)+k(71)*n(8)*n(7)+k(160)*n(20)*n(42)...
    -k(30)*n(4)*n(21)-k(60)*n(21)*n(42)-k(76)*n(1)*n(21)...
    -k(82)*n(4)*n(21)-k(91)*n(7)*n(21)...
    -k(95)*n(15)*n(21)-k(99)*n(21)*n(26)-k(100)*n(21)*n(22);
    %Chemin CO:
    dndt(22)= k(4)*n(1)*n(37)+k(5)*n(2)*n(15)+k(8)*n(1)*n(15)...
    +k(24)*n(4)*n(36)+k(26)*n(4)*n(26)+k(27)*n(4)*n(13)...
    +k(28)*n(4)*n(35)+k(29)*n(4)*n(9)+k(31)*n(4)*n(14)...
    +k(32)*n(4)*n(12)+2*k(47)*n(9)*n(15)+k(49)*n(9)*n(19)...
    +k(50)*n(11)*n(15)+k(54)*n(15)*n(27)+k(58)*n(24)*n(37)...
    +k(65)*n(1)*n(4)+k(75)*n(1)*n(23)+k(83)*n(4)*n(23)...
    +k(92)*n(9)*n(23)+k(96)*n(15)*n(23)+k(102)*n(24)*n(23)...
    +k(139)*n(1)*n(6)+k(140)*n(3)*n(4)-k(13)*n(1)*n(22)...
    -k(55)*n(15)*n(22)-k(79)*n(5)*n(22)-k(88)*n(8)*n(22)...
    -k(100)*n(21)*n(22)-k(101)*n(22)*n(27)...
    -k(130)*n(22)*n(43)-k(146)*n(6)*n(22);
    %Chemin CO+:
    dndt(23)= k(7)*n(2)*n(36)+k(15)*n(5)*n(9)+k(33)*n(4)*n(11)...
    +k(49)*n(9)*n(19)+k(50)*n(11)*n(15)+k(63)*n(2)*n(4)...
    +k(67)*n(1)*n(5)+k(79)*n(5)*n(22)+k(88)*n(8)*n(22)...
    +k(100)*n(21)*n(22)+k(101)*n(22)*n(27)-k(75)*n(1)*n(23)...
    -k(83)*n(4)*n(23)-k(92)*n(9)*n(23)...
    -k(96)*n(15)*n(23)-k(102)*n(24)*n(23);
    %Chemin NO:
    dndt(24)= k(20)*n(4)*n(32)+k(22)*2*n(4)*n(30)+k(23)*n(4)*n(20)...
    +k(24)*n(4)*n(36)+k(25)*n(4)*n(26)+k(34)*n(8)*n(15)...
    +k(35)*n(7)*n(15)+k(51)*n(15)*n(37)+k(106)*n(3)*n(25)...
    +k(114)*n(14)*n(25)+k(115)*n(12)*n(25)+k(116)*n(18)*n(25)...
    +k(147)*n(6)*n(7)-k(10)*n(1)*n(24)-k(14)*n(3)*n(24)-k(19)*n(4)*n(24)...
    -k(37)*n(7)*n(24)-k(52)*n(15)*n(24)-k(57)*2*n(24)*n(24)...
    -k(58)*n(24)*n(37)-k(72)*n(2)*n(24)-k(90)*n(8)*n(24)...
    -k(94)*n(11)*n(24)-k(98)*n(19)*n(24)-k(102)*n(24)*n(23)...
    -k(103)*n(24)*n(27)-k(131)*n(24)*n(43)-k(145)*n(6)*n(24);
    %Chemin NO+:
    dndt(25)= k(6)*n(2)*n(30)+k(16)*n(5)*n(26)+k(17)*n(5)*n(32)...
    +k(18)*n(5)*n(20)+k(30)*n(4)*n(21)+k(36)*n(7)*n(19)...
    +k(54)*n(15)*n(27)+k(72)*n(2)*n(24)+k(90)*n(8)*n(24)...
    +k(94)*n(11)*n(24)+k(98)*n(19)*n(24)+k(102)*n(24)*n(23)...
    +k(103)*n(24)*n(27)-k(106)*n(3)*n(25)-k(114)*n(14)*n(25)...
    -k(115)*n(12)*n(25)-k(116)*n(18)*n(25);
    %Chemin CN:
    dndt(26)= k(4)*n(1)*n(37)+k(6)*n(2)*n(30)+k(7)*n(2)*n(36)...
    +k(10)*n(1)*n(24)+k(11)*n(1)*n(20)+k(12)*n(1)*n(35)...
    +k(28)*n(4)*n(35)+k(41)*n(7)*n(13)+k(42)*n(7)*n(14)...
    +k(43)*n(7)*n(9)+k(45)*n(7)*n(11)+k(46)*2*n(7)*n(35)...
    +k(69)*n(1)*n(7)+k(73)*n(1)*n(27)+k(93)*n(9)*n(27)...
    +k(101)*n(22)*n(27)+k(103)*n(24)*n(27)...
    +k(112)*n(8)*n(28)+k(142)*n(3)*n(7)...
    +k(164)*n(7)*n(1)*n(1)-k(9)*n(1)*n(26)-k(16)*n(5)*n(26)...
    -k(25)*n(4)*n(26)-k(26)*n(4)*n(26)-k(39)*n(7)*n(26)...
    -k(56)*n(15)*n(26)-k(85)*n(6)*n(26)...
    -k(89)*n(8)*n(26)-k(99)*n(21)*n(26);
    %Chemin CN+:
    dndt(27)= k(62)*n(2)*n(7)+k(89)*n(8)*n(26)+k(99)*n(21)*n(26)...
    -k(40)*n(7)*n(27)-k(54)*n(15)*n(27)...
    -k(73)*n(1)*n(27)-k(93)*n(9)*n(27)...
    -k(101)*n(22)*n(27)-k(103)*n(24)*n(27);
    %Chemin CN-:
    dndt(28)= k(14)*n(3)*n(24)+k(44)*n(7)*n(12)+k(85)*n(6)*n(26)-k(112)*n(8)*n(28);
    %Chemin O3:
    dndt(29)= k(163)*n(1)*n(15)*n(4);
    %Chemin N2O:
    dndt(30)= k(38)*n(7)*n(32)+k(53)*n(15)*n(20)+k(58)*n(24)*n(37)...
    -k(6)*n(2)*n(30)-k(21)*n(4)*n(30)...
    -k(22)*n(4)*n(30)-k(80)*n(5)*n(30);
    %Chemin N2O+:
    dndt(31)= k(80)*n(5)*n(30);
    %Chemin NO2:
    dndt(32)= k(52)*n(15)*n(24)+k(145)*n(6)*n(24)-k(17)*n(5)*n(32)...
    -k(20)*n(4)*n(32)-k(38)*n(7)*n(32)...
    -k(97)*n(19)*n(32)-k(165)*n(1)*n(32)*n(4);
    %Chemin NO2+:
    dndt(33)= k(97)*n(19)*n(32);
    %Chemin CO2:
    dndt(34)= k(51)*n(15)*n(37)+k(55)*n(15)*n(22)...
    +k(141)*n(3)*n(15)+k(146)*n(6)*n(22);
    %Chemin C2N:
    dndt(35)= -k(12)*n(1)*n(35)-k(28)*n(4)*n(35)-k(46)*n(7)*n(35);
    %Chemin CNO:
    dndt(36)= -k(7)*n(2)*n(36)-k(24)*n(4)*n(36);
    %Chemin OCN
    dndt(37)= -k(4)*n(1)*n(37)-k(51)*n(15)*n(37)...
    -k(58)*n(24)*n(37)+k(56)*n(15)*n(26);
    %Chemin NO3:
    dndt(38)= k(165)*n(1)*n(32)*n(4);
    %Chemin C2O:
    dndt(39)= k(144)*n(4)*n(12);
    %Chemin C4:
    dndt(40)= k(64)*n(1)*n(13)+k(138)*n(1)*n(14)+k(148)*n(9)*n(12);
    %Chemin C5:
    dndt(41)= k(149)*n(12)*n(13);
    %Chemin e:
    dndt(42)= k(118)*n(1)*n(43)+k(119)*n(3)*n(43)+k(120)*n(6)*n(43)...
    +k(121)*n(12)*n(43)+k(123)*n(12)*n(43)+k(125)*n(14)*n(43)...
    +k(127)*n(18)*n(43)+k(136)*n(1)*n(3)+k(137)*n(1)*n(12)...
    +k(138)*n(1)*n(14)+k(139)*n(1)*n(6)+k(140)*n(3)*n(4)...
    +k(141)*n(3)*n(15)+k(142)*n(3)*n(7)+k(143)*n(4)*n(6)...
    +k(144)*n(4)*n(12)+k(145)*n(6)*n(24)+k(146)*n(6)*n(22)...
    +k(147)*n(6)*n(7)+k(148)*n(9)*n(12)+k(149)*n(12)*n(13)...
    +k(150)*n(9)*n(42)+k(151)*n(15)*n(42)+k(152)*n(15)*n(42)...
    +k(153)*n(10)*n(42)+k(156)*n(18)*n(15)+k(157)*n(18)*n(16)...
    +k(158)*n(18)*n(17)+k(159)*n(20)*n(42)...
    +2*k(160)*n(20)*n(42)-k(1)*n(2)*n(42)-k(2)*n(5)*n(42)...
    -k(3)*n(8)*n(42)-k(59)*n(11)*n(42)-k(60)*n(21)*n(42)...
    -k(132)*n(1)*n(42)-k(133)*n(4)*n(42)-k(134)*n(9)*n(42)...
    -k(135)*n(13)*n(42)-k(150)*n(9)*n(42)-k(151)*n(15)*n(42)...
    -k(152)*n(15)*n(42)-k(153)*n(10)*n(42)-k(159)*n(20)*n(42)...
    -k(160)*n(20)*n(42)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15);
    %Chemin hv:
    dndt(43)= -k(118)*n(1)*n(43)-k(119)*n(3)*n(43)-k(120)*n(6)*n(43)...
    -k(121)*n(9)*n(43)-k(122)*n(11)*n(43)-k(123)*n(12)*n(43)...
    -k(124)*n(13)*n(43)-k(125)*n(14)*n(43)-k(126)*n(15)*n(43)...
    -k(127)*n(18)*n(43)-k(128)*n(19)*n(43)-k(129)*n(20)*n(43)...
    -k(130)*n(22)*n(43)-k(131)*n(24)*n(43)+k(1)*n(2)*n(42)...
    +k(2)*n(5)*n(42)+k(3)*n(8)*n(42)+k(61)*n(2)*n(1)...
    +k(62)*n(2)*n(7)+k(63)*n(2)*n(4)+k(64)*n(1)*n(13)...
    +k(65)*n(1)*n(4)+k(66)*n(1)*n(9)+k(67)*n(1)*n(5)...
    +k(68)*n(1)*n(1)+k(69)*n(1)*n(7)+k(70)*n(4)*n(4)...
    +k(71)*n(8)*n(7)+k(132)*n(1)*n(42)+k(133)*n(4)*n(42)...
    +k(134)*n(9)*n(42)+k(135)*n(13)*n(42);
 end 
 

call_myequations.m

 function [t,n]= call_myequations()
    %condition intial 
    %IC=IC1=IC2=....=IC43=10^6cm^3
    IC =10^19*ones(1,43);
    tic,[t,n] = ode45( t,n,k,[0 10e-3],IC);
    toc
    plot(t,n)
 end
 

main.m

 function main
k = parameters();
[t,n] = call_myequations(k);

I got this error: when ecuting main.m

Error using call_myequations
Too many input arguments.

Error in main (line 3)
[t,n] = call_myequations(k);

Answer 1

[Star Strider]

The ‘call_my_equations’ function needs to be changed to call ‘myequations’ correctly, and to reference ‘parameters’ in order to supply ‘k’.

The online Run feature times out before completing this, however this runs without error.

[t,n]= call_myequations;

Elapsed time is 0.880277 seconds.

function [t,n]= call_myequations()
    %condition intial 
    %IC=IC1=IC2=....=IC43=10^6cm^3
    k= parameters;
    IC =10^19*ones(1,43);
    tic,[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 10e-3],IC);
    toc
    plot(t,n)
 end

EDIT — (1 Jun 2021 at 13:42)

Changed to ode15s.

function k= parameters
     %constantes :::
    %temperature
    T=17000;
    %constante k: recombinaison radiative
    k(1)=2.36 *(1e-12)*(T/300)^(-0.29)*exp(17.60/T); 
    k(2) =3.24 *(1e-12) *(T/300)^(-0.66);
    k(3)=3.50 *(1e-12)*(T/300)^(-0.53)* exp(3.20/T);
    %constante k::Reconmbinaison 
    k(4)=1*10^(-10); 
    k(5)=4.54*10^(-10);
    k(6)=9.10*10^(-10);
    k(7)=8.98*10^(-9)*(T/300)^(-0.5);
    k(8)=5.56*10^(-11)*(T/300)^(0.41)*exp(26.90/T);
    k(9)=4.98*10^(-10)*exp(-18116/T);
    k(10)=6*10^(-11)*(T/300)^(-0.16);
    k(11)=8.69*10^(-11)*exp(-22600/T);
    k(12)=1*10^(-10);
    k(13)=2.94*10^(-11)*(T/300)^(0.5)*exp(-58025/T);
    k(14)=1*10^(-9);
    k(15)=4.80*10^(-10);
    k(16)=1.00*10^(-9)*(T/300)^(-0.5);
    k(17)=8.30*10^(-10);
    k(18)=2.42*10^(-12)*(T/300)^(-0.21)*exp(44/T);
    k(19)=1.18*10^(-11)*exp(-20413/T);
    k(20)=9.82*10^(-12)*(T/300)^(-0.21)*exp(-5.20/T);
    k(21)=1.66*10^(-10)*exp(-14100/T);
    k(22)=1.15*10^(-10)*exp(-13400/T);
    k(23)=2.51*10^(-10)*exp(-38602/T);
    k(24)=1.00*10^(-10);
    k(25)=5.37*10^(-11)*exp(-13800/T);
    k(26)=5.00*10^(-11)*exp(-200/T);
    k(27)=5.00*10^(-12)*exp(-900/T);
    k(28)=6.00*10^(-12);
    k(29)=2.00*10^(-10)*(T/300)^(-0.12);
    k(30)=1.30*10^(-10);
    k(31)=5.00*10^(-10);
    k(32)=2.90*10^(-10);
    k(33)=3.10*10^(-10);
    k(34)=3.66*10^(-11);
    k(35)=2.26*10^(-12)*(T/300)^(0.86)*exp(-3134.0/T);
    k(36)=1.80*10^(-10);
    k(37)=3.38*10^(-11)*(T/300)^(-0.17)*exp(2.80/T);
    k(38)=3*10^(-12);
    k(39)=1.00*10^(-10)*(T/300)^(0.4);
    k(40)=6.10*10^(-10);
    k(41)=1*10^(-13);
    k(42)=5*10^(-11);
    k(43)=5*10^(-11);
    k(44)=1.15*10^(-10);
    k(45)=4*10^(-11);
    k(46)=1*10^(-10);
    k(47)=1.50*10^(-11)*exp(-4300/T);
    k(48)=5.3137*10^(-10);
    k(49)=4.10*10^(-10);
    k(50)=8*10^(-10);
    k(51)=1.32*10^(-12);
    k(52)=2.8*10^(-12)*exp(23400/T);
    k(53)=1*10^(-10)*exp(-55200/T);
    k(54)=8.60*10^(-11);
    k(55)=5.99*10^(-12)*exp(-24075/T);
    k(56)=2.02*10^(-11)*(T/300)^(-0.19)*exp(31.90/T);
    k(57)=2.51*10^(-11)*exp(-30653/T);
    k(58)=3.47*10^(-11)*(T/300)^(-1.33)*exp(-242/T);
    
    %constante k::Recombinaison dissociative
    k(59)=3*10^(-7)*(T/300)^(-0.5);
    k(60)=1.70*10^(-7)*(T/300)^(-0.3);
     
    %constante k::Association radiative
    k(61)=4.01*10^(-18)*(T/300)^(0.17)*exp(-101.50/T);
    k(62)=1.08*10^(-18)*(T/300)^(0.07)*exp(-57.50/T);
    k(63)=3.14*10^(-18)*(T/300)^(-0.15)*exp(-68/T);
    k(64)=4.00*10^(-14)*(T/300)^(-1);
    k(65)=4.69*10^(-19)*(T/300)^(1.52)*exp(50.50/T);
    k(66)=3*10^(-16)*(T/300)^(-1);
    k(67)=5*10^(-10)*(T/300)^(-3.7)*exp(-800/T);
    k(68)=4.36*10^(-18)*(T/300)^(0.35)*exp(-161.30/T);
    k(69)=5.72*10^(-19)*(T/300)^(0.37)*exp(-51/T);
    k(70)=4.90*10^(-20)*(T/300)^(1.58);
    k(71)=3.71*10^(-18)*(T/300)^(0.24)*exp(-26.10/T);
    
    %constante k::Transfert de charge
    k(72)=7.05*10^(-10)*(T/300)^(-0.03)*exp(16.70/T);
    k(73)=1*10^(-10);
    k(74)=1.10*10^(-10);
    k(75)=1.10*10^(-10);
    k(76)=1.10*10^(-10);
    k(77)=5.20*10^(-11);
    k(78)=4.80*10^(-10);
    k(79)=4.90*10^(-12)*(T/300)^(0.5)*exp(-4580/T);
    k(80)=6.30*10^(-10);
    k(81)=1.90*10^(-11);
    k(82)=1.00*10^(-11);
    k(83)=1.40*10^(-10);
    k(84)=7.30*10^(-10)*exp(-890/T);
    k(85)=1.00*10^(-9);
    k(86)=3.11*10^(-10);
    k(87)=1.00*10^(-9);
    k(88)=8.25*10^(-10);
    k(89)=1.10*10^(-9)*(T/300)^(-0.50);
    k(90)=4.51*10^(-10);
    k(91)=1* 10^(-11);
    k(92)=8.40*10^(-10);
    k(93)=8.50^(-10);
    k(94)=3.40*10^(-10);
    k(95)=5*10^(-11);
    k(96)=1.20*10^(-10);
    k(97)=6.60*10^(-10);
    k(98)=4.60*10^(-10);
    k(99)=1.00*10^(-10)*(T/300)^(-0.5);
    k(100)=7.40*10^(-11);
    k(101)=6.30*10^(-10);
    k(102)=3.30*10^(-10);
    k(103)=5.70*10^(-10);
    
    %constante k::Neutralisation
    k(104)=7.51*10^(-8)*(T/300)^(-0.5);
    k(105)=7.51*10^(-8)*(T/300)^(-0.5);
    k(106)=7.51*10^(-8)*(T/300)^(-0.5);
    k(107)=7.51*10^(-8)*(T/300)^(-0.5);
    k(108)=7.51*10^(-8)*(T/300)^(-0.5);
    k(109)=7.51*10^(-8)*(T/300)^(-0.5);
    k(110)=7.51*10^(-8)*(T/300)^(-0.50);
    k(111)=7.51*10^(-8)*(T/300)^(-0.50);
    k(112)=7.51*10^(-8)*(T/300)^(-0.50);
    k(113)=7.51*10^(-8)*(T/300)^(-0.50);
    k(114)=7.51*10^(-8)*(T/300)^(-0.5);
    k(115)=7.51*10^(-8)*(T/300)^(-0.5);
    k(116)=7.51*10^(-8)*(T/300)^(-0.5);
    k(117)=7.51*10^(-8)*(T/300)^(-0.5);
    %constante k::Photoionisation
    k(118)=3.10*10^(-10)*exp(-3.3/T);
    k(119)=4.90*10^(-8)*exp(-0.50/T);
    k(120)=1.09*10^(-8)*exp(-0.5/T);
    k(121)=4.10*10^(-10)*exp(-3.80/T);
    k(122)=1*10^(-11)*exp(-1.70/T);
    k(123)=2.96*10^(-9)*exp(-2/T);
    k(124)=5*10^(-9)*exp(-2.1/T);
    k(125)=6.88*10^(-9)*exp(-1.5/T);
    k(126)=7.90*10^(-10)*exp(-2.1/T);
    k(127)=6.10*10^(-9)*exp(-0.5/T);
    k(128)=3.50*10^(-11)*exp(-2/T);
    k(129)=2.30*10^(-10)*exp(-3.90/T);
    k(130)=2.10*10^(-10)*exp(-3.50/T);
    k(131)=4.70*10^(-10)*exp(-2.10/T);
    %constante k::Attachement radiatif
    
    k(132)=2.25*10^(-15);
    k(133)=1.50*10^(-15);
    k(134)=2.00*10^(-15)*(T/300)^(-0.5);
    k(135)=1.70*10^(-14)*(T/300)^(-0.5);
    %constante k::Détachement associatif
    
    k(136)=5*10^(-10);
    k(137)=1*10^(-9);
    k(138)=1*10^(-9);
    k(139)=5*10^(-10);
    k(140)=5*10^(-10);
    k(141)=5*10^(-11);
    k(142)=5*10^(-10);
    k(143)=1.90*10^(-10);
    k(144)=2.90*10^(-10);
    k(145)=3.10*10^(-10)*(T/300)^(-0.83);
    k(146)=6.50*10^(-10);
    k(147)=2.20*10^(-10);
    k(148)=1.00*10^(-9);
    k(149)=1*10^(-9);
    
    %constante k::Excitation/ Désexcitation
    k(150)=2.0425*10^(-7)*T^(-0.3)*exp(-29489.55/T);
    k(151)=2.6*10^(-10)*T^(0.89)*exp(-11372.9/T);
    k(152)=1.2*10^(-10)*T^(0.59)*exp(-18916.15/T);
    k(153)=6.0278*10^(-7)*T^(-0.4)*exp(-19880.36/T);
    k(154)=1.6306*10^(-17);
    k(155)=3.6532*10^(-11);
    
    %constante k::Détachement électronique
    k(156)=7*10^(-10);
    k(157)=2*10^(-10);
    k(158)=3.6*10^(-10);
    
    %constante k::Dissociation par impact électronique
    k(159)= 6.2*10^(-10) * T^(1.75) * exp(-113032.7/T );
    
    %constante k::Ionisation par impact électronique
    k(160)= 7.1*10^(-9) * T^(0.76) * exp(-181038/T );
    
    %constante k::Attachement électronique cn(1)6/s
    k(161)=1*10^(-31);
    k(162)=1*10^(-31);
    
    %constante k::Processus de recombinaison cn(1)6/s
    k(163)=5.70*10^(-34)*T^(-2.6);
    k(164)=2.5919*10^(-80);
    k(165)=1.31*10^(-31)*T^(-1.5);
end 

function dndt = myequations(t,n,k)
    dndt = zeros(43,1);
    %Chemin C
    dndt(1)= k(1)*n(2)*n(42)+k(15)*n(5)*n(9)+k(16)*n(5)*n(26)+k(25)*n(4)*n(26)...
    +k(29)*n(4)*n(9)+k(33)*n(4)*n(11)+k(39)*n(7)*n(26)+k(40)*n(7)*n(27)...
    +k(43)*n(7)*n(9)+k(44)*n(7)*n(12)+k(48)*n(9)*n(9)+k(59)*n(11)*n(42)...
    +k(72)*n(2)*n(24)+k(104)*n(3)*n(5)+k(105)*n(3)*n(8)+k(106)*n(3)*n(25)...
    +k(117)*n(18)*n(2)+k(119)*n(3)*n(43)+k(122)*n(11)*n(43)+k(124)*n(13)*n(43)...
    +k(130)*n(22)*n(43)-k(4)*n(1)*n(37)-k(8)*n(1)*n(15)-k(9)*n(1)*n(26)...
    -k(10)*n(1)*n(24)-k(11)*n(1)*n(20)-k(12)*n(1)*n(35)-k(13)*n(1)*n(22)...
    -k(61)*n(2)*n(1)-k(64)*n(1)*n(13)-k(65)*n(1)*n(4)-k(66)*n(1)*n(9)...
    -k(67)*n(1)*n(5)-k(68)*n(1)*n(1)-k(69)*n(1)*n(7)-k(73)*n(1)*n(27)...
    -k(74)*n(1)*n(11)-k(75)*n(1)*n(23)-k(76)*n(1)*n(21)-k(77)*n(1)*n(19)...
    -k(118)*n(1)*n(43)-k(132)*n(1)*n(42)-k(136)*n(1)*n(3)-k(137)*n(1)*n(12)...
    -k(138)*n(1)*n(14)-k(139)*n(1)*n(6)-2*k(164)*n(7)*n(1)*n(1)...
    +k(164)*n(7)*n(1)^(2)-k(163)*n(1)*n(15)*n(4)+k(163)*n(1)*n(15)*n(4)...
    -k(165)*n(1)*n(32)*n(4)+k(165)*n(1)*n(32)*n(4);
    %Chemin C+:
    dndt(2)= -k(1)*n(2)*n(42)-k(5)*n(2)*n(15)-k(6)*n(2)*n(30)-k(7)*n(2)*n(36)...
    -k(61)*n(2)*n(1)-k(62)*n(2)*n(7)-k(63)*n(2)*n(4)-k(72)*n(2)*n(24)...
    -k(117)*n(18)*n(2)+k(45)*n(7)*n(11)+k(73)*n(1)*n(27)...
    +k(74)*n(1)*n(11)+k(75)*n(1)*n(23)+k(76)*n(1)*n(21)...
    +k(77)*n(1)*n(19)+k(118)*n(1)*n(43)+k(122)*n(11)*n(43);
    %Chemin C-:
    dndt(3)= k(32)*n(4)*n(12)+k(132)*n(1)*n(42)-k(14)*n(3)*n(24)-k(104)*n(3)*n(5)...
    -k(105)*n(3)*n(8)-k(106)*n(3)*n(25)-k(119)*n(3)*n(43)...
    -k(136)*n(1)*n(3)-k(140)*n(3)*n(4)-k(141)*n(3)*n(15)-k(142)*n(3)*n(7);
    %Chemin O:
    dndt(4)= k(2)*n(5)*n(42)+k(8)*n(1)*n(15)+k(10)*n(1)*n(24)+k(13)*n(1)*n(22)...
    +k(14)*n(3)*n(24)+k(35)*n(7)*n(15)+k(36)*n(7)*n(19)...
    +k(37)*n(7)*n(24)+k(38)*n(7)*n(32)+k(52)*n(15)*n(24)...
    +k(53)*n(15)*n(20)+k(55)*n(15)*n(22)+k(56)*n(15)*n(26)...
    +k(78)*n(5)*n(9)+k(79)*n(5)*n(22)+k(80)*n(5)*n(30)...
    +k(81)*n(5)*n(15)+k(84)*n(6)*n(15)+k(85)*n(6)*n(26)...
    +k(104)*n(3)*n(5)+2*k(107)*n(5)*n(6)+k(108)*n(5)*n(18)...
    +k(109)*n(6)*n(8)+k(120)*n(6)*n(43)+2*k(126)*n(15)*n(43)...
    +k(128)*n(19)*n(43)+k(130)*n(22)*n(43)+k(131)*n(24)*n(43)...
    +k(162)*n(42)*n(4)*n(15)-k(19)*n(4)*n(24)-k(20)*n(4)*n(32)...
    -k(21)*n(4)*n(30)-k(22)*n(4)*n(30)-k(23)*n(4)*n(20)...
    -k(24)*n(4)*n(36)-k(25)*n(4)*n(26)-k(26)*n(4)*n(26)...
    -k(27)*n(4)*n(13)-k(28)*n(4)*n(35)-k(29)*n(4)*n(9)...
    -k(30)*n(4)*n(21)-k(31)*n(4)*n(14)-k(32)*n(4)*n(12)...
    -k(33)*n(4)*n(11)-k(63)*n(2)*n(4)-k(65)*n(1)*n(4)...
    -2*k(70)*n(4)*n(4)-k(82)*n(4)*n(21)-k(83)*n(4)*n(23)...
    -k(133)*n(4)*n(42)-k(140)*n(3)*n(4)-k(143)*n(4)*n(6)...
    -k(144)*n(4)*n(12)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15)-k(163)*n(1)*n(15)*n(4)-k(165)*n(1)*n(32)*n(4);
    %Chemin O+:
    dndt(5)= +k(5)*n(2)*n(15)+k(34)*n(8)*n(15)+k(82)*n(4)*n(21)...
    +k(83)*n(4)*n(23)+k(128)*n(19)*n(43)-k(2)*n(5)*n(42)-k(15)*n(5)*n(9)...
    -k(16)*n(5)*n(26)-k(17)*n(5)*n(32)-k(18)*n(5)*n(20)-k(67)*n(1)*n(5)...
    -k(78)*n(5)*n(9)-k(79)*n(5)*n(22)-k(80)*n(5)*n(30)-k(81)*n(5)*n(15)...
    -k(104)*n(3)*n(5)-k(107)*n(5)*n(6)-k(108)*n(5)*n(18);
    %Chemin O-:
    dndt(6)= k(133)*n(4)*n(42)+k(161)*n(42)*n(4)*n(15)-k(84)*n(6)*n(15)...
    -k(85)*n(6)*n(26)-k(107)*n(5)*n(6)-k(109)*n(6)*n(8)-k(120)*n(6)*n(43)...
    -k(139)*n(1)*n(6)-k(143)*n(4)*n(6)-k(145)*n(6)*n(24)...
    -k(146)*n(6)*n(22)-k(147)*n(6)*n(7);
    %Chemin N:
    dndt(7)= k(3)*n(8)*n(42)+k(9)*n(1)*n(26)+k(11)*n(1)*n(20)...
    +k(18)*n(5)*n(20)+k(19)*n(4)*n(24)+k(23)*n(4)*n(20)...
    +k(26)*n(4)*n(26)+k(30)*n(4)*n(21)+2*k(60)*n(21)*n(42)...
    +k(86)*n(8)*n(15)+k(87)*n(8)*n(9)+k(88)*n(8)*n(22)...
    +k(89)*n(8)*n(26)+k(90)*n(8)*n(24)+k(105)*n(3)*n(8)...
    +k(109)*n(6)*n(8)+k(110)*n(8)*n(18)+k(111)*n(8)*n(12)...
    +k(112)*n(8)*n(28)+k(113)*n(8)*n(14)+2*k(129)*n(20)*n(43)...
    +k(131)*n(24)*n(43)+2*k(159)*n(20)*n(42)-k(35)*n(7)*n(15)...
    -k(36)*n(7)*n(19)-k(37)*n(7)*n(24)-k(38)*n(7)*n(32)...
    -k(39)*n(7)*n(26)-k(40)*n(7)*n(27)-k(41)*n(7)*n(13)...
    -k(42)*n(7)*n(14)-k(43)*n(7)*n(9)-k(44)*n(7)*n(12)...
    -k(45)*n(7)*n(11)-k(46)*n(7)*n(35)-k(62)*n(2)*n(7)...
    -k(69)*n(1)*n(7)-k(71)*n(8)*n(7)-k(91)*n(7)*n(21)...
    -k(142)*n(3)*n(7)-k(147)*n(6)*n(7)-k(164)*n(7)*n(1)*n(1);
    %Chemin N+:
    dndt(8)=k(91)*n(7)*n(21)-k(3)*n(8)*n(42)-k(34)*n(8)*n(15)...
    -k(71)*n(8)*n(7)-k(86)*n(8)*n(15)-k(87)*n(8)*n(9)...
    -k(88)*n(8)*n(22)-k(89)*n(8)*n(26)-k(90)*n(8)*n(24)...
    -k(105)*n(3)*n(8)-k(109)*n(6)*n(8)-k(110)*n(8)*n(18)...
    -k(111)*n(8)*n(12)-k(112)*n(8)*n(28)-k(113)*n(8)*n(14);
    %Chemin C2:
    dndt(9)= k(9)*n(1)*n(26)+k(12)*n(1)*n(35)+k(13)*n(1)*n(22)...
    +k(27)*n(4)*n(13)+k(41)*n(7)*n(13)+k(68)*n(1)*n(1)...
    +k(74)*n(1)*n(11)+k(94)*n(11)*n(24)+k(111)*n(8)*n(12)...
    +k(115)*n(12)*n(25)+k(121)*n(9)*n(43)+k(123)*n(12)*n(43)...
    +k(124)*n(13)*n(43)+k(136)*n(1)*n(3)+k(153)*n(10)*n(42)...
    +k(154)*n(10)^(2)+k(155)*n(10)*n(1)-k(15)*n(5)*n(9)...
    -k(29)*n(4)*n(9)-k(43)*n(7)*n(9)-k(47)*n(9)*n(15)...
    -2*k(48)*n(9)*n(9)-k(49)*n(9)*n(19)-k(66)*n(1)*n(9)...
    -k(78)*n(5)*n(9)-k(87)*n(8)*n(9)-k(92)*n(9)*n(23)...
    -k(93)*n(9)*n(27)-k(121)*n(9)*n(43)-k(134)*n(9)*n(42)...
    -k(148)*n(9)*n(12)-k(150)*n(9)*n(42);
    %Chemin C2(d)
    dndt(10)= k(150)*n(9)*n(42)-k(153)*n(10)*n(42)-k(154)*n(10)^(2)-k(155)*n(10)*n(1);
    %Chemin C2+:
    dndt(11)= k(61)*n(2)*n(1)+k(78)*n(5)*n(9)+k(87)*n(8)*n(9)...
    +k(92)*n(9)*n(23)+k(93)*n(9)*n(27)-k(33)*n(4)*n(11)...
    -k(45)*n(7)*n(11)-k(50)*n(11)*n(15)-k(59)*n(11)*n(42)...
    -k(74)*n(1)*n(11)-k(94)*n(11)*n(24)-k(122)*n(11)*n(43);
    %Chemin C2-:
    dndt(12)= -k(32)*n(4)*n(12)-k(44)*n(7)*n(12)-k(111)*n(8)*n(12)...
    -k(115)*n(12)*n(25)-k(123)*n(12)*n(43)-k(137)*n(1)*n(12)...
    -k(144)*n(4)*n(12)-k(148)*n(9)*n(12)-k(149)*n(12)*n(13)...
    +k(31)*n(4)*n(14)+k(42)*n(7)*n(14)+k(134)*n(9)*n(42);
    %Chemin C3:
    dndt(13)= k(48)*n(9)*n(9)+k(66)*n(1)*n(9)+k(113)*n(8)*n(14)...
    +k(114)*n(14)*n(25)+k(125)*n(14)*n(43)...
    +k(137)*n(1)*n(12)-k(27)*n(4)*n(13)-k(41)*n(7)*n(13)...
    -k(64)*n(1)*n(13)-k(124)*n(13)*n(43)...
    -k(135)*n(13)*n(42)-k(149)*n(12)*n(13);
    %Chemin C3-:
    dndt(14)= k(135)*n(13)*n(42)-k(31)*n(4)*n(14)...
    -k(42)*n(7)*n(14)-k(113)*n(8)*n(14)...
    -k(114)*n(14)*n(25)-k(125)*n(14)*n(43)-k(138)*n(1)*n(14);
    %Chemin O2:
    dndt(15)= k(17)*n(5)*n(32)+k(19)*n(4)*n(24)+k(20)*n(4)*n(32)...
    +k(21)*n(4)*n(30)+k(57)*n(24)*n(24)+k(70)*n(4)*n(4)...
    +k(77)*n(1)*n(19)+k(97)*n(19)*n(32)+k(98)*n(19)*n(24)...
    +k(108)*n(5)*n(18)+k(110)*n(8)*n(18)+k(116)*n(18)*n(25)...
    +k(117)*n(18)*n(2)+k(127)*n(18)*n(43)+k(143)*n(4)*n(6)...
    +2*k(156)*n(18)*n(15)+2*k(157)*n(18)*n(16)...
    +2*k(158)*n(18)*n(17)+k(161)*n(42)*n(4)*n(15)...
    -k(5)*n(2)*n(15)-k(8)*n(1)*n(15)-k(34)*n(8)*n(15)...
    -k(35)*n(7)*n(15)-k(47)*n(9)*n(15)-k(50)*n(11)*n(15)...
    -k(51)*n(15)*n(37)-k(52)*n(15)*n(24)-k(53)*n(15)*n(20)...
    -k(54)*n(15)*n(27)-k(55)*n(15)*n(22)-k(56)*n(15)*n(26)...
    -k(81)*n(5)*n(15)-k(84)*n(6)*n(15)-k(86)*n(8)*n(15)...
    -k(95)*n(15)*n(21)-k(96)*n(15)*n(23)-k(126)*n(15)*n(43)...
    -k(141)*n(3)*n(15)-k(151)*n(15)*n(42)-k(152)*n(15)*n(42)...
    -k(156)*n(18)*n(15)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15)-k(163)*n(1)*n(15)*n(4);
    %Chemin O2(a):
    dndt(16)= k(151)*n(42)*n(15)-k(157)*n(16)*n(18);
    %Chemin O2(b):
    dndt(17)= k(152)*n(42)*n(15)-k(158)*n(17)*n(18);
    %Chemin O2-:
    dndt(18)= k(84)*n(6)*n(15)+k(162)*n(42)*n(4)*n(15)...
    -k(108)*n(5)*n(18)-k(110)*n(8)*n(18)...
    -k(116)*n(18)*n(25)-k(117)*n(18)*n(2)-k(127)*n(18)*n(43)...
    -k(156)*n(18)*n(15)-k(157)*n(18)*n(16)-k(158)*n(18)*n(17);
    %Chemin O2+:
    dndt(19)= k(81)*n(5)*n(15)+k(86)*n(8)*n(15)+k(95)*n(15)*n(21)...
    +k(96)*n(15)*n(23)-k(36)*n(7)*n(19)-k(49)*n(9)*n(19)-k(77)*n(1)*n(19)...
    -k(97)*n(19)*n(32)-k(98)*n(19)*n(24)-k(128)*n(19)*n(43);
    %Chemin N2:
    dndt(20)= k(21)*n(4)*n(30)+k(37)*n(7)*n(24)+k(39)*n(7)*n(26)...
    +k(57)*n(24)*n(24)+k(76)*n(1)*n(21)+k(82)*n(4)*n(21)...
    +k(91)*n(7)*n(21)+k(95)*n(15)*n(21)+k(99)*n(21)*n(26)...
    +k(100)*n(21)*n(22)-k(11)*n(1)*n(20)-k(18)*n(5)*n(20)...
    -k(23)*n(4)*n(20)-k(53)*n(15)*n(20)...
    -k(129)*n(20)*n(43)-k(159)*n(20)*n(42)-k(160)*n(20)*n(42);
    %Chemin N2+:
    dndt(21)= k(40)*n(7)*n(27)+k(71)*n(8)*n(7)+k(160)*n(20)*n(42)...
    -k(30)*n(4)*n(21)-k(60)*n(21)*n(42)-k(76)*n(1)*n(21)...
    -k(82)*n(4)*n(21)-k(91)*n(7)*n(21)...
    -k(95)*n(15)*n(21)-k(99)*n(21)*n(26)-k(100)*n(21)*n(22);
    %Chemin CO:
    dndt(22)= k(4)*n(1)*n(37)+k(5)*n(2)*n(15)+k(8)*n(1)*n(15)...
    +k(24)*n(4)*n(36)+k(26)*n(4)*n(26)+k(27)*n(4)*n(13)...
    +k(28)*n(4)*n(35)+k(29)*n(4)*n(9)+k(31)*n(4)*n(14)...
    +k(32)*n(4)*n(12)+2*k(47)*n(9)*n(15)+k(49)*n(9)*n(19)...
    +k(50)*n(11)*n(15)+k(54)*n(15)*n(27)+k(58)*n(24)*n(37)...
    +k(65)*n(1)*n(4)+k(75)*n(1)*n(23)+k(83)*n(4)*n(23)...
    +k(92)*n(9)*n(23)+k(96)*n(15)*n(23)+k(102)*n(24)*n(23)...
    +k(139)*n(1)*n(6)+k(140)*n(3)*n(4)-k(13)*n(1)*n(22)...
    -k(55)*n(15)*n(22)-k(79)*n(5)*n(22)-k(88)*n(8)*n(22)...
    -k(100)*n(21)*n(22)-k(101)*n(22)*n(27)...
    -k(130)*n(22)*n(43)-k(146)*n(6)*n(22);
    %Chemin CO+:
    dndt(23)= k(7)*n(2)*n(36)+k(15)*n(5)*n(9)+k(33)*n(4)*n(11)...
    +k(49)*n(9)*n(19)+k(50)*n(11)*n(15)+k(63)*n(2)*n(4)...
    +k(67)*n(1)*n(5)+k(79)*n(5)*n(22)+k(88)*n(8)*n(22)...
    +k(100)*n(21)*n(22)+k(101)*n(22)*n(27)-k(75)*n(1)*n(23)...
    -k(83)*n(4)*n(23)-k(92)*n(9)*n(23)...
    -k(96)*n(15)*n(23)-k(102)*n(24)*n(23);
    %Chemin NO:
    dndt(24)= k(20)*n(4)*n(32)+k(22)*2*n(4)*n(30)+k(23)*n(4)*n(20)...
    +k(24)*n(4)*n(36)+k(25)*n(4)*n(26)+k(34)*n(8)*n(15)...
    +k(35)*n(7)*n(15)+k(51)*n(15)*n(37)+k(106)*n(3)*n(25)...
    +k(114)*n(14)*n(25)+k(115)*n(12)*n(25)+k(116)*n(18)*n(25)...
    +k(147)*n(6)*n(7)-k(10)*n(1)*n(24)-k(14)*n(3)*n(24)-k(19)*n(4)*n(24)...
    -k(37)*n(7)*n(24)-k(52)*n(15)*n(24)-k(57)*2*n(24)*n(24)...
    -k(58)*n(24)*n(37)-k(72)*n(2)*n(24)-k(90)*n(8)*n(24)...
    -k(94)*n(11)*n(24)-k(98)*n(19)*n(24)-k(102)*n(24)*n(23)...
    -k(103)*n(24)*n(27)-k(131)*n(24)*n(43)-k(145)*n(6)*n(24);
    %Chemin NO+:
    dndt(25)= k(6)*n(2)*n(30)+k(16)*n(5)*n(26)+k(17)*n(5)*n(32)...
    +k(18)*n(5)*n(20)+k(30)*n(4)*n(21)+k(36)*n(7)*n(19)...
    +k(54)*n(15)*n(27)+k(72)*n(2)*n(24)+k(90)*n(8)*n(24)...
    +k(94)*n(11)*n(24)+k(98)*n(19)*n(24)+k(102)*n(24)*n(23)...
    +k(103)*n(24)*n(27)-k(106)*n(3)*n(25)-k(114)*n(14)*n(25)...
    -k(115)*n(12)*n(25)-k(116)*n(18)*n(25);
    %Chemin CN:
    dndt(26)= k(4)*n(1)*n(37)+k(6)*n(2)*n(30)+k(7)*n(2)*n(36)...
    +k(10)*n(1)*n(24)+k(11)*n(1)*n(20)+k(12)*n(1)*n(35)...
    +k(28)*n(4)*n(35)+k(41)*n(7)*n(13)+k(42)*n(7)*n(14)...
    +k(43)*n(7)*n(9)+k(45)*n(7)*n(11)+k(46)*2*n(7)*n(35)...
    +k(69)*n(1)*n(7)+k(73)*n(1)*n(27)+k(93)*n(9)*n(27)...
    +k(101)*n(22)*n(27)+k(103)*n(24)*n(27)...
    +k(112)*n(8)*n(28)+k(142)*n(3)*n(7)...
    +k(164)*n(7)*n(1)*n(1)-k(9)*n(1)*n(26)-k(16)*n(5)*n(26)...
    -k(25)*n(4)*n(26)-k(26)*n(4)*n(26)-k(39)*n(7)*n(26)...
    -k(56)*n(15)*n(26)-k(85)*n(6)*n(26)...
    -k(89)*n(8)*n(26)-k(99)*n(21)*n(26);
    %Chemin CN+:
    dndt(27)= k(62)*n(2)*n(7)+k(89)*n(8)*n(26)+k(99)*n(21)*n(26)...
    -k(40)*n(7)*n(27)-k(54)*n(15)*n(27)...
    -k(73)*n(1)*n(27)-k(93)*n(9)*n(27)...
    -k(101)*n(22)*n(27)-k(103)*n(24)*n(27);
    %Chemin CN-:
    dndt(28)= k(14)*n(3)*n(24)+k(44)*n(7)*n(12)+k(85)*n(6)*n(26)-k(112)*n(8)*n(28);
    %Chemin O3:
    dndt(29)= k(163)*n(1)*n(15)*n(4);
    %Chemin N2O:
    dndt(30)= k(38)*n(7)*n(32)+k(53)*n(15)*n(20)+k(58)*n(24)*n(37)...
    -k(6)*n(2)*n(30)-k(21)*n(4)*n(30)...
    -k(22)*n(4)*n(30)-k(80)*n(5)*n(30);
    %Chemin N2O+:
    dndt(31)= k(80)*n(5)*n(30);
    %Chemin NO2:
    dndt(32)= k(52)*n(15)*n(24)+k(145)*n(6)*n(24)-k(17)*n(5)*n(32)...
    -k(20)*n(4)*n(32)-k(38)*n(7)*n(32)...
    -k(97)*n(19)*n(32)-k(165)*n(1)*n(32)*n(4);
    %Chemin NO2+:
    dndt(33)= k(97)*n(19)*n(32);
    %Chemin CO2:
    dndt(34)= k(51)*n(15)*n(37)+k(55)*n(15)*n(22)...
    +k(141)*n(3)*n(15)+k(146)*n(6)*n(22);
    %Chemin C2N:
    dndt(35)= -k(12)*n(1)*n(35)-k(28)*n(4)*n(35)-k(46)*n(7)*n(35);
    %Chemin CNO:
    dndt(36)= -k(7)*n(2)*n(36)-k(24)*n(4)*n(36);
    %Chemin OCN
    dndt(37)= -k(4)*n(1)*n(37)-k(51)*n(15)*n(37)...
    -k(58)*n(24)*n(37)+k(56)*n(15)*n(26);
    %Chemin NO3:
    dndt(38)= k(165)*n(1)*n(32)*n(4);
    %Chemin C2O:
    dndt(39)= k(144)*n(4)*n(12);
    %Chemin C4:
    dndt(40)= k(64)*n(1)*n(13)+k(138)*n(1)*n(14)+k(148)*n(9)*n(12);
    %Chemin C5:
    dndt(41)= k(149)*n(12)*n(13);
    %Chemin e:
    dndt(42)= k(118)*n(1)*n(43)+k(119)*n(3)*n(43)+k(120)*n(6)*n(43)...
    +k(121)*n(12)*n(43)+k(123)*n(12)*n(43)+k(125)*n(14)*n(43)...
    +k(127)*n(18)*n(43)+k(136)*n(1)*n(3)+k(137)*n(1)*n(12)...
    +k(138)*n(1)*n(14)+k(139)*n(1)*n(6)+k(140)*n(3)*n(4)...
    +k(141)*n(3)*n(15)+k(142)*n(3)*n(7)+k(143)*n(4)*n(6)...
    +k(144)*n(4)*n(12)+k(145)*n(6)*n(24)+k(146)*n(6)*n(22)...
    +k(147)*n(6)*n(7)+k(148)*n(9)*n(12)+k(149)*n(12)*n(13)...
    +k(150)*n(9)*n(42)+k(151)*n(15)*n(42)+k(152)*n(15)*n(42)...
    +k(153)*n(10)*n(42)+k(156)*n(18)*n(15)+k(157)*n(18)*n(16)...
    +k(158)*n(18)*n(17)+k(159)*n(20)*n(42)...
    +2*k(160)*n(20)*n(42)-k(1)*n(2)*n(42)-k(2)*n(5)*n(42)...
    -k(3)*n(8)*n(42)-k(59)*n(11)*n(42)-k(60)*n(21)*n(42)...
    -k(132)*n(1)*n(42)-k(133)*n(4)*n(42)-k(134)*n(9)*n(42)...
    -k(135)*n(13)*n(42)-k(150)*n(9)*n(42)-k(151)*n(15)*n(42)...
    -k(152)*n(15)*n(42)-k(153)*n(10)*n(42)-k(159)*n(20)*n(42)...
    -k(160)*n(20)*n(42)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15);
    %Chemin hv:
    dndt(43)= -k(118)*n(1)*n(43)-k(119)*n(3)*n(43)-k(120)*n(6)*n(43)...
    -k(121)*n(9)*n(43)-k(122)*n(11)*n(43)-k(123)*n(12)*n(43)...
    -k(124)*n(13)*n(43)-k(125)*n(14)*n(43)-k(126)*n(15)*n(43)...
    -k(127)*n(18)*n(43)-k(128)*n(19)*n(43)-k(129)*n(20)*n(43)...
    -k(130)*n(22)*n(43)-k(131)*n(24)*n(43)+k(1)*n(2)*n(42)...
    +k(2)*n(5)*n(42)+k(3)*n(8)*n(42)+k(61)*n(2)*n(1)...
    +k(62)*n(2)*n(7)+k(63)*n(2)*n(4)+k(64)*n(1)*n(13)...
    +k(65)*n(1)*n(4)+k(66)*n(1)*n(9)+k(67)*n(1)*n(5)...
    +k(68)*n(1)*n(1)+k(69)*n(1)*n(7)+k(70)*n(4)*n(4)...
    +k(71)*n(8)*n(7)+k(132)*n(1)*n(42)+k(133)*n(4)*n(42)...
    +k(134)*n(9)*n(42)+k(135)*n(13)*n(42);
 end 

.

[Imene Yed]

thnk you it's working but I dont know why it still running and didnt get me results.

[Star Strider]

The system is ‘stiff’.

I changed the solver to ode15s in my edit just now, and it works.

[Imene Yed]

thnk you so much for your help.

[Star Strider]

As always, my pleasure!

[Imene Yed]

@Star Strider I'm trying to add this line expr_at_i(:)= n(:,1) + 2*n(:,2); I dont know where I should put it sorry I'm really new to this.

[Star Strider]

No worries!

It would seem logical that it would go after the ode15s call in ‘call_myequations’ as:

expr_at_i = n(:,1) + 2*n(:,2);

The (:) is not necessary (and will likely throw an error), because the right-hand-side will return a column vector.

If you want to plot it, it would be necessary to add it to the plot, either with a hold call and separate plot call afterwards:

figure
plot(t,n)
hold on
plot(t, expr_at_i)
hold off

or concatenate it to ‘n’ as:

n = [n expr_at_i];
figure
plot(t, n)

and then plot it as originally described.

That decision is yours.

.

[Imene Yed]

@Star Strider expr_at_i = n(:,1) + 2*n(:,2); yes I have put it after the call ode15s but it dosent seem to give me the results ,in my workspace I didnt find the values of expr_at_i

[Star Strider]

Just after it is calculated, add this line temporarily:

ViewVector = expr_at_i(1:10)

That will tell you something about what it is, and if it is being calculated correctly. (Note that there is no ending semicolon, so it should print to the Command Window.) Delete that line when you are satisfied that it is being created.

.

[Imene Yed]

sorry for disturbing you again but is it possible to add this function I mean to call it in mycall_equation script

  function res = func(t,n)
       
            %Conservation de la masse:
%C:
masse_C(:)= n(:,1) + 2*n(:,9)+2*n(:,10)+3*n(:,13)+...
    n(:,22)+n(:,26)+n(:,34)+2*n(:,35)+n(:,36)...
    +n(:,37)+2*n(:,39)+4*n(:,40)+5*n(:,41)
%O:
masse_O(:)= n(:,4)+2*n(:,15)+2*n(:,16)+2*n(:,17)+n(:,22)...
    +n(:,24)+ 3*n(:,29)+n(:,30)+2*n(:,32)+2*n(:,34)+...
    n(:,36)+n(:,37)+3*n(:,38)+ n(:,39);
 %N 
 masse_N(:)= n(:,7)+2*n(:,20)+n(:,24)+n(:,26)+2*n(:,30)+n(:,32)...
     +n(:,35)+n(:,36)+n(:,37)+n(:,38)

[Star Strider]

Yes, however it needs to be changed slightly —

function res = func(t,n)
       
            %Conservation de la masse:
%C:
masse_C = n(:,1) + 2*n(:,9)+2*n(:,10)+3*n(:,13)+...
    n(:,22)+n(:,26)+n(:,34)+2*n(:,35)+n(:,36)...
    +n(:,37)+2*n(:,39)+4*n(:,40)+5*n(:,41);
%O:
masse_O = n(:,4)+2*n(:,15)+2*n(:,16)+2*n(:,17)+n(:,22)...
    +n(:,24)+ 3*n(:,29)+n(:,30)+2*n(:,32)+2*n(:,34)+...
    n(:,36)+n(:,37)+3*n(:,38)+ n(:,39);
 %N 
masse_N = n(:,7)+2*n(:,20)+n(:,24)+n(:,26)+2*n(:,30)+n(:,32)...
     +n(:,35)+n(:,36)+n(:,37)+n(:,38);
 
res = [masse_C, masse_O,  masse_N];
end

Then call it as:

res = func(t,n);

The output has to be assigned in the function or its absence will throw an error.

That should work because all the components of ‘res’ are column vectors (or so it appears).

.

[Imene Yed]

function [t,n]= call_myequations()
    %condition intial 
    %IC=IC1=IC2=....=IC43=10^6cm^3
    k= parameters;
    IC =10^19*ones(1,43);
    tic,[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-3],IC);
    toc
    res = func(t,n);
     
    plot(t,n)
Unrecognized function or variable 'func'.

Error in call_myequations (line 8)
    res = func(t,n);

[Star Strider]

The ‘func’ function need to be placed at the end of the code, anywhere after ‘call_myequations’. It must either be in the same .m file or in a .m file all its own, named func.m.

[Imene Yed]

I still cant see the res in my workspace .I dont understand why there isnt any value of res thnk you

[Star Strider]

Try it without the ending semicolon:

res = func(t,n)

the columns of ‘res’ should appear in the Command Window.

[Imene Yed]

yes I saw it ,is it possible to make it in an excel file because I need this data

[Star Strider]

I would do something like this —

res = func(t,n)

res_table = array2table(res, 'VariableNames',{'masse_C', 'masse_O', 'masse_N'})

writetable(res_table, 'YourFileName.xlsx')

This will also display the table so that you can see that it was created, and what it looks like.

It should be written to the Excel file the same way, including the column headings.

.

[Imene Yed]

thnk you so much !!!!!!!!!!!!!!!!

[Star Strider]

As always, my pleasure!

[Imene Yed]

@Star Strider sorry for disturbing you again I'm trying to display the value of t and n since I needed it in my thesis because this script display only my graphes.

function [t,n]= call_myequations(k)

    %condition intial 
    %IC=IC1=IC2=....=IC43=10^6cm^3
    k= parameters;
    IC =10^18*ones(1,43);
    tic,[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-3],IC);
    toc
    res = func(t,n)
     

res_table = array2table(res, 'VariableNames',{'conservation_masse_C', 'conservation_masse_O', 'conservation_masse_N','conservation_charge_positive','conservation_charge_negative'})

writetable(res_table, 'myFileName.xlsx')
     plot(t,n)

[Star Strider]

I would do the same sort of thing for those —

tn_table = array2table([t n], 'VariableNames',{'conservation_masse_C', 'conservation_masse_O', 'conservation_masse_N','conservation_charge_positive','conservation_charge_negative'})
writetable(tn_table, 'myFileName.xlsx')

However it will be necessary to create a new vector of variable names for the ‘n’ variables, perhaps 'Chemin C', 'Chemin C+', (or 'Chemin_C', 'Chemin_Cplus', depending on the MATLAB release you are using) and so for the rest, otherwise the function will automatically created names for them, such as 'n_1', 'n_2', and so forth. The first variable name obviously will be 'Time'.

.

[Imene Yed]

I know their named it's not really a pb for me ,can I just display them without writing all their names in the variablesNames???

[Walter Roberson]

When you writetable() you can use 'writevariablenames', false

However if you are going to do that, then there is not a lot of point in going to table. You might as well leave res as an array and use writematrix()

[Star Strider]

@Imene Yed — As Walter noted, yes, you can. I suggested the table because that would make understanding the results for all 43 columns of ‘n’ (and additional columns of data derived from them) and referring to them easier for others who might look at the table or want to reference the data in it. That is a significant advantage of using table arrays.

[Imene Yed]

Thank for the explanations.can you please write to me the syntax that would be the most simple .I'm still trying to get used to all of this .I'm really new to programming

[Star Strider]

As always, my pleasure!

The syntax already presented is the simplest syntax. It would be necessary to put the variable names in a cell array, and then add them as appropriate arguments to the table function call when the table is created. If you do not want to use a table, see the documentation for writematrix to create the file.

[Imene Yed]

@Star Strider can you help me please I tried to add some initial conditions to my program here is what I got

unction [t,n]= call_myequations(k)

    %condition intial 
    %IC=IC1=IC2=....=IC43=10^6cm^3
    k= parameters;
    
    tic,[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-8],[1.0487*10^21 0 0 0 0 0 0 0 1.0487*10^21...
        0 0 0 1.0487*10^21 0 2*10^16 0 0 0 0 7*10^15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0])
    toc
    res = func(t,n)
     

res_table = array2table(res, 'VariableNames',{'conservation_masse_C', 'conservation_masse_O', 'conservation_masse_N','conservation_charge_positive','conservation_charge_negative'})

writetable(res_table, 'myFileName.xlsx')

errors

>> call_myequations

Error: File: call_myequations.m Line: 7 Column: 99

Invalid use of operator.

>> call_myequations

Error: File: call_myequations.m Line: 7 Column: 99

Invalid use of operator.

[Star Strider]

Put a space:

% tic,[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-8],[1.0487*10^21 0 0 0 0 0 0 0 1.0487*10^21 ...
%                                                                                              ↑ ← HERE

and it works!

(I commented it out because otherwise it would execute out-of-sequence here, and throw an error in the online Run feature.)

I thought that was the problem, so just to be certain I am testing it here —

[t,n]= call_myequations;

Elapsed time is 0.397776 seconds.

function [t,n]= call_myequations()
    %condition intial 
    %IC=IC1=IC2=....=IC43=10^6cm^3
    k= parameters;
    IC =10^19*ones(1,43);
    tic
    [t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-8],[1.0487*10^21 0 0 0 0 0 0 0 1.0487*10^21 ...
        0 0 0 1.0487*10^21 0 2*10^16 0 0 0 0 7*10^15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]);
    toc
    plot(t,n)
 end

And it was!

function k= parameters
     %constantes :::
    %temperature
    T=17000;
    %constante k: recombinaison radiative
    k(1)=2.36 *(1e-12)*(T/300)^(-0.29)*exp(17.60/T); 
    k(2) =3.24 *(1e-12) *(T/300)^(-0.66);
    k(3)=3.50 *(1e-12)*(T/300)^(-0.53)* exp(3.20/T);
    %constante k::Reconmbinaison 
    k(4)=1*10^(-10); 
    k(5)=4.54*10^(-10);
    k(6)=9.10*10^(-10);
    k(7)=8.98*10^(-9)*(T/300)^(-0.5);
    k(8)=5.56*10^(-11)*(T/300)^(0.41)*exp(26.90/T);
    k(9)=4.98*10^(-10)*exp(-18116/T);
    k(10)=6*10^(-11)*(T/300)^(-0.16);
    k(11)=8.69*10^(-11)*exp(-22600/T);
    k(12)=1*10^(-10);
    k(13)=2.94*10^(-11)*(T/300)^(0.5)*exp(-58025/T);
    k(14)=1*10^(-9);
    k(15)=4.80*10^(-10);
    k(16)=1.00*10^(-9)*(T/300)^(-0.5);
    k(17)=8.30*10^(-10);
    k(18)=2.42*10^(-12)*(T/300)^(-0.21)*exp(44/T);
    k(19)=1.18*10^(-11)*exp(-20413/T);
    k(20)=9.82*10^(-12)*(T/300)^(-0.21)*exp(-5.20/T);
    k(21)=1.66*10^(-10)*exp(-14100/T);
    k(22)=1.15*10^(-10)*exp(-13400/T);
    k(23)=2.51*10^(-10)*exp(-38602/T);
    k(24)=1.00*10^(-10);
    k(25)=5.37*10^(-11)*exp(-13800/T);
    k(26)=5.00*10^(-11)*exp(-200/T);
    k(27)=5.00*10^(-12)*exp(-900/T);
    k(28)=6.00*10^(-12);
    k(29)=2.00*10^(-10)*(T/300)^(-0.12);
    k(30)=1.30*10^(-10);
    k(31)=5.00*10^(-10);
    k(32)=2.90*10^(-10);
    k(33)=3.10*10^(-10);
    k(34)=3.66*10^(-11);
    k(35)=2.26*10^(-12)*(T/300)^(0.86)*exp(-3134.0/T);
    k(36)=1.80*10^(-10);
    k(37)=3.38*10^(-11)*(T/300)^(-0.17)*exp(2.80/T);
    k(38)=3*10^(-12);
    k(39)=1.00*10^(-10)*(T/300)^(0.4);
    k(40)=6.10*10^(-10);
    k(41)=1*10^(-13);
    k(42)=5*10^(-11);
    k(43)=5*10^(-11);
    k(44)=1.15*10^(-10);
    k(45)=4*10^(-11);
    k(46)=1*10^(-10);
    k(47)=1.50*10^(-11)*exp(-4300/T);
    k(48)=5.3137*10^(-10);
    k(49)=4.10*10^(-10);
    k(50)=8*10^(-10);
    k(51)=1.32*10^(-12);
    k(52)=2.8*10^(-12)*exp(23400/T);
    k(53)=1*10^(-10)*exp(-55200/T);
    k(54)=8.60*10^(-11);
    k(55)=5.99*10^(-12)*exp(-24075/T);
    k(56)=2.02*10^(-11)*(T/300)^(-0.19)*exp(31.90/T);
    k(57)=2.51*10^(-11)*exp(-30653/T);
    k(58)=3.47*10^(-11)*(T/300)^(-1.33)*exp(-242/T);
    
    %constante k::Recombinaison dissociative
    k(59)=3*10^(-7)*(T/300)^(-0.5);
    k(60)=1.70*10^(-7)*(T/300)^(-0.3);
     
    %constante k::Association radiative
    k(61)=4.01*10^(-18)*(T/300)^(0.17)*exp(-101.50/T);
    k(62)=1.08*10^(-18)*(T/300)^(0.07)*exp(-57.50/T);
    k(63)=3.14*10^(-18)*(T/300)^(-0.15)*exp(-68/T);
    k(64)=4.00*10^(-14)*(T/300)^(-1);
    k(65)=4.69*10^(-19)*(T/300)^(1.52)*exp(50.50/T);
    k(66)=3*10^(-16)*(T/300)^(-1);
    k(67)=5*10^(-10)*(T/300)^(-3.7)*exp(-800/T);
    k(68)=4.36*10^(-18)*(T/300)^(0.35)*exp(-161.30/T);
    k(69)=5.72*10^(-19)*(T/300)^(0.37)*exp(-51/T);
    k(70)=4.90*10^(-20)*(T/300)^(1.58);
    k(71)=3.71*10^(-18)*(T/300)^(0.24)*exp(-26.10/T);
    
    %constante k::Transfert de charge
    k(72)=7.05*10^(-10)*(T/300)^(-0.03)*exp(16.70/T);
    k(73)=1*10^(-10);
    k(74)=1.10*10^(-10);
    k(75)=1.10*10^(-10);
    k(76)=1.10*10^(-10);
    k(77)=5.20*10^(-11);
    k(78)=4.80*10^(-10);
    k(79)=4.90*10^(-12)*(T/300)^(0.5)*exp(-4580/T);
    k(80)=6.30*10^(-10);
    k(81)=1.90*10^(-11);
    k(82)=1.00*10^(-11);
    k(83)=1.40*10^(-10);
    k(84)=7.30*10^(-10)*exp(-890/T);
    k(85)=1.00*10^(-9);
    k(86)=3.11*10^(-10);
    k(87)=1.00*10^(-9);
    k(88)=8.25*10^(-10);
    k(89)=1.10*10^(-9)*(T/300)^(-0.50);
    k(90)=4.51*10^(-10);
    k(91)=1* 10^(-11);
    k(92)=8.40*10^(-10);
    k(93)=8.50^(-10);
    k(94)=3.40*10^(-10);
    k(95)=5*10^(-11);
    k(96)=1.20*10^(-10);
    k(97)=6.60*10^(-10);
    k(98)=4.60*10^(-10);
    k(99)=1.00*10^(-10)*(T/300)^(-0.5);
    k(100)=7.40*10^(-11);
    k(101)=6.30*10^(-10);
    k(102)=3.30*10^(-10);
    k(103)=5.70*10^(-10);
    
    %constante k::Neutralisation
    k(104)=7.51*10^(-8)*(T/300)^(-0.5);
    k(105)=7.51*10^(-8)*(T/300)^(-0.5);
    k(106)=7.51*10^(-8)*(T/300)^(-0.5);
    k(107)=7.51*10^(-8)*(T/300)^(-0.5);
    k(108)=7.51*10^(-8)*(T/300)^(-0.5);
    k(109)=7.51*10^(-8)*(T/300)^(-0.5);
    k(110)=7.51*10^(-8)*(T/300)^(-0.50);
    k(111)=7.51*10^(-8)*(T/300)^(-0.50);
    k(112)=7.51*10^(-8)*(T/300)^(-0.50);
    k(113)=7.51*10^(-8)*(T/300)^(-0.50);
    k(114)=7.51*10^(-8)*(T/300)^(-0.5);
    k(115)=7.51*10^(-8)*(T/300)^(-0.5);
    k(116)=7.51*10^(-8)*(T/300)^(-0.5);
    k(117)=7.51*10^(-8)*(T/300)^(-0.5);
    %constante k::Photoionisation
    k(118)=3.10*10^(-10)*exp(-3.3/T);
    k(119)=4.90*10^(-8)*exp(-0.50/T);
    k(120)=1.09*10^(-8)*exp(-0.5/T);
    k(121)=4.10*10^(-10)*exp(-3.80/T);
    k(122)=1*10^(-11)*exp(-1.70/T);
    k(123)=2.96*10^(-9)*exp(-2/T);
    k(124)=5*10^(-9)*exp(-2.1/T);
    k(125)=6.88*10^(-9)*exp(-1.5/T);
    k(126)=7.90*10^(-10)*exp(-2.1/T);
    k(127)=6.10*10^(-9)*exp(-0.5/T);
    k(128)=3.50*10^(-11)*exp(-2/T);
    k(129)=2.30*10^(-10)*exp(-3.90/T);
    k(130)=2.10*10^(-10)*exp(-3.50/T);
    k(131)=4.70*10^(-10)*exp(-2.10/T);
    %constante k::Attachement radiatif
    
    k(132)=2.25*10^(-15);
    k(133)=1.50*10^(-15);
    k(134)=2.00*10^(-15)*(T/300)^(-0.5);
    k(135)=1.70*10^(-14)*(T/300)^(-0.5);
    %constante k::Détachement associatif
    
    k(136)=5*10^(-10);
    k(137)=1*10^(-9);
    k(138)=1*10^(-9);
    k(139)=5*10^(-10);
    k(140)=5*10^(-10);
    k(141)=5*10^(-11);
    k(142)=5*10^(-10);
    k(143)=1.90*10^(-10);
    k(144)=2.90*10^(-10);
    k(145)=3.10*10^(-10)*(T/300)^(-0.83);
    k(146)=6.50*10^(-10);
    k(147)=2.20*10^(-10);
    k(148)=1.00*10^(-9);
    k(149)=1*10^(-9);
    
    %constante k::Excitation/ Désexcitation
    k(150)=2.0425*10^(-7)*T^(-0.3)*exp(-29489.55/T);
    k(151)=2.6*10^(-10)*T^(0.89)*exp(-11372.9/T);
    k(152)=1.2*10^(-10)*T^(0.59)*exp(-18916.15/T);
    k(153)=6.0278*10^(-7)*T^(-0.4)*exp(-19880.36/T);
    k(154)=1.6306*10^(-17);
    k(155)=3.6532*10^(-11);
    
    %constante k::Détachement électronique
    k(156)=7*10^(-10);
    k(157)=2*10^(-10);
    k(158)=3.6*10^(-10);
    
    %constante k::Dissociation par impact électronique
    k(159)= 6.2*10^(-10) * T^(1.75) * exp(-113032.7/T );
    
    %constante k::Ionisation par impact électronique
    k(160)= 7.1*10^(-9) * T^(0.76) * exp(-181038/T );
    
    %constante k::Attachement électronique cn(1)6/s
    k(161)=1*10^(-31);
    k(162)=1*10^(-31);
    
    %constante k::Processus de recombinaison cn(1)6/s
    k(163)=5.70*10^(-34)*T^(-2.6);
    k(164)=2.5919*10^(-80);
    k(165)=1.31*10^(-31)*T^(-1.5);
end 

function dndt = myequations(t,n,k)
    dndt = zeros(43,1);
    %Chemin C
    dndt(1)= k(1)*n(2)*n(42)+k(15)*n(5)*n(9)+k(16)*n(5)*n(26)+k(25)*n(4)*n(26)...
    +k(29)*n(4)*n(9)+k(33)*n(4)*n(11)+k(39)*n(7)*n(26)+k(40)*n(7)*n(27)...
    +k(43)*n(7)*n(9)+k(44)*n(7)*n(12)+k(48)*n(9)*n(9)+k(59)*n(11)*n(42)...
    +k(72)*n(2)*n(24)+k(104)*n(3)*n(5)+k(105)*n(3)*n(8)+k(106)*n(3)*n(25)...
    +k(117)*n(18)*n(2)+k(119)*n(3)*n(43)+k(122)*n(11)*n(43)+k(124)*n(13)*n(43)...
    +k(130)*n(22)*n(43)-k(4)*n(1)*n(37)-k(8)*n(1)*n(15)-k(9)*n(1)*n(26)...
    -k(10)*n(1)*n(24)-k(11)*n(1)*n(20)-k(12)*n(1)*n(35)-k(13)*n(1)*n(22)...
    -k(61)*n(2)*n(1)-k(64)*n(1)*n(13)-k(65)*n(1)*n(4)-k(66)*n(1)*n(9)...
    -k(67)*n(1)*n(5)-k(68)*n(1)*n(1)-k(69)*n(1)*n(7)-k(73)*n(1)*n(27)...
    -k(74)*n(1)*n(11)-k(75)*n(1)*n(23)-k(76)*n(1)*n(21)-k(77)*n(1)*n(19)...
    -k(118)*n(1)*n(43)-k(132)*n(1)*n(42)-k(136)*n(1)*n(3)-k(137)*n(1)*n(12)...
    -k(138)*n(1)*n(14)-k(139)*n(1)*n(6)-2*k(164)*n(7)*n(1)*n(1)...
    +k(164)*n(7)*n(1)^(2)-k(163)*n(1)*n(15)*n(4)+k(163)*n(1)*n(15)*n(4)...
    -k(165)*n(1)*n(32)*n(4)+k(165)*n(1)*n(32)*n(4);
    %Chemin C+:
    dndt(2)= -k(1)*n(2)*n(42)-k(5)*n(2)*n(15)-k(6)*n(2)*n(30)-k(7)*n(2)*n(36)...
    -k(61)*n(2)*n(1)-k(62)*n(2)*n(7)-k(63)*n(2)*n(4)-k(72)*n(2)*n(24)...
    -k(117)*n(18)*n(2)+k(45)*n(7)*n(11)+k(73)*n(1)*n(27)...
    +k(74)*n(1)*n(11)+k(75)*n(1)*n(23)+k(76)*n(1)*n(21)...
    +k(77)*n(1)*n(19)+k(118)*n(1)*n(43)+k(122)*n(11)*n(43);
    %Chemin C-:
    dndt(3)= k(32)*n(4)*n(12)+k(132)*n(1)*n(42)-k(14)*n(3)*n(24)-k(104)*n(3)*n(5)...
    -k(105)*n(3)*n(8)-k(106)*n(3)*n(25)-k(119)*n(3)*n(43)...
    -k(136)*n(1)*n(3)-k(140)*n(3)*n(4)-k(141)*n(3)*n(15)-k(142)*n(3)*n(7);
    %Chemin O:
    dndt(4)= k(2)*n(5)*n(42)+k(8)*n(1)*n(15)+k(10)*n(1)*n(24)+k(13)*n(1)*n(22)...
    +k(14)*n(3)*n(24)+k(35)*n(7)*n(15)+k(36)*n(7)*n(19)...
    +k(37)*n(7)*n(24)+k(38)*n(7)*n(32)+k(52)*n(15)*n(24)...
    +k(53)*n(15)*n(20)+k(55)*n(15)*n(22)+k(56)*n(15)*n(26)...
    +k(78)*n(5)*n(9)+k(79)*n(5)*n(22)+k(80)*n(5)*n(30)...
    +k(81)*n(5)*n(15)+k(84)*n(6)*n(15)+k(85)*n(6)*n(26)...
    +k(104)*n(3)*n(5)+2*k(107)*n(5)*n(6)+k(108)*n(5)*n(18)...
    +k(109)*n(6)*n(8)+k(120)*n(6)*n(43)+2*k(126)*n(15)*n(43)...
    +k(128)*n(19)*n(43)+k(130)*n(22)*n(43)+k(131)*n(24)*n(43)...
    +k(162)*n(42)*n(4)*n(15)-k(19)*n(4)*n(24)-k(20)*n(4)*n(32)...
    -k(21)*n(4)*n(30)-k(22)*n(4)*n(30)-k(23)*n(4)*n(20)...
    -k(24)*n(4)*n(36)-k(25)*n(4)*n(26)-k(26)*n(4)*n(26)...
    -k(27)*n(4)*n(13)-k(28)*n(4)*n(35)-k(29)*n(4)*n(9)...
    -k(30)*n(4)*n(21)-k(31)*n(4)*n(14)-k(32)*n(4)*n(12)...
    -k(33)*n(4)*n(11)-k(63)*n(2)*n(4)-k(65)*n(1)*n(4)...
    -2*k(70)*n(4)*n(4)-k(82)*n(4)*n(21)-k(83)*n(4)*n(23)...
    -k(133)*n(4)*n(42)-k(140)*n(3)*n(4)-k(143)*n(4)*n(6)...
    -k(144)*n(4)*n(12)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15)-k(163)*n(1)*n(15)*n(4)-k(165)*n(1)*n(32)*n(4);
    %Chemin O+:
    dndt(5)= +k(5)*n(2)*n(15)+k(34)*n(8)*n(15)+k(82)*n(4)*n(21)...
    +k(83)*n(4)*n(23)+k(128)*n(19)*n(43)-k(2)*n(5)*n(42)-k(15)*n(5)*n(9)...
    -k(16)*n(5)*n(26)-k(17)*n(5)*n(32)-k(18)*n(5)*n(20)-k(67)*n(1)*n(5)...
    -k(78)*n(5)*n(9)-k(79)*n(5)*n(22)-k(80)*n(5)*n(30)-k(81)*n(5)*n(15)...
    -k(104)*n(3)*n(5)-k(107)*n(5)*n(6)-k(108)*n(5)*n(18);
    %Chemin O-:
    dndt(6)= k(133)*n(4)*n(42)+k(161)*n(42)*n(4)*n(15)-k(84)*n(6)*n(15)...
    -k(85)*n(6)*n(26)-k(107)*n(5)*n(6)-k(109)*n(6)*n(8)-k(120)*n(6)*n(43)...
    -k(139)*n(1)*n(6)-k(143)*n(4)*n(6)-k(145)*n(6)*n(24)...
    -k(146)*n(6)*n(22)-k(147)*n(6)*n(7);
    %Chemin N:
    dndt(7)= k(3)*n(8)*n(42)+k(9)*n(1)*n(26)+k(11)*n(1)*n(20)...
    +k(18)*n(5)*n(20)+k(19)*n(4)*n(24)+k(23)*n(4)*n(20)...
    +k(26)*n(4)*n(26)+k(30)*n(4)*n(21)+2*k(60)*n(21)*n(42)...
    +k(86)*n(8)*n(15)+k(87)*n(8)*n(9)+k(88)*n(8)*n(22)...
    +k(89)*n(8)*n(26)+k(90)*n(8)*n(24)+k(105)*n(3)*n(8)...
    +k(109)*n(6)*n(8)+k(110)*n(8)*n(18)+k(111)*n(8)*n(12)...
    +k(112)*n(8)*n(28)+k(113)*n(8)*n(14)+2*k(129)*n(20)*n(43)...
    +k(131)*n(24)*n(43)+2*k(159)*n(20)*n(42)-k(35)*n(7)*n(15)...
    -k(36)*n(7)*n(19)-k(37)*n(7)*n(24)-k(38)*n(7)*n(32)...
    -k(39)*n(7)*n(26)-k(40)*n(7)*n(27)-k(41)*n(7)*n(13)...
    -k(42)*n(7)*n(14)-k(43)*n(7)*n(9)-k(44)*n(7)*n(12)...
    -k(45)*n(7)*n(11)-k(46)*n(7)*n(35)-k(62)*n(2)*n(7)...
    -k(69)*n(1)*n(7)-k(71)*n(8)*n(7)-k(91)*n(7)*n(21)...
    -k(142)*n(3)*n(7)-k(147)*n(6)*n(7)-k(164)*n(7)*n(1)*n(1);
    %Chemin N+:
    dndt(8)=k(91)*n(7)*n(21)-k(3)*n(8)*n(42)-k(34)*n(8)*n(15)...
    -k(71)*n(8)*n(7)-k(86)*n(8)*n(15)-k(87)*n(8)*n(9)...
    -k(88)*n(8)*n(22)-k(89)*n(8)*n(26)-k(90)*n(8)*n(24)...
    -k(105)*n(3)*n(8)-k(109)*n(6)*n(8)-k(110)*n(8)*n(18)...
    -k(111)*n(8)*n(12)-k(112)*n(8)*n(28)-k(113)*n(8)*n(14);
    %Chemin C2:
    dndt(9)= k(9)*n(1)*n(26)+k(12)*n(1)*n(35)+k(13)*n(1)*n(22)...
    +k(27)*n(4)*n(13)+k(41)*n(7)*n(13)+k(68)*n(1)*n(1)...
    +k(74)*n(1)*n(11)+k(94)*n(11)*n(24)+k(111)*n(8)*n(12)...
    +k(115)*n(12)*n(25)+k(121)*n(9)*n(43)+k(123)*n(12)*n(43)...
    +k(124)*n(13)*n(43)+k(136)*n(1)*n(3)+k(153)*n(10)*n(42)...
    +k(154)*n(10)^(2)+k(155)*n(10)*n(1)-k(15)*n(5)*n(9)...
    -k(29)*n(4)*n(9)-k(43)*n(7)*n(9)-k(47)*n(9)*n(15)...
    -2*k(48)*n(9)*n(9)-k(49)*n(9)*n(19)-k(66)*n(1)*n(9)...
    -k(78)*n(5)*n(9)-k(87)*n(8)*n(9)-k(92)*n(9)*n(23)...
    -k(93)*n(9)*n(27)-k(121)*n(9)*n(43)-k(134)*n(9)*n(42)...
    -k(148)*n(9)*n(12)-k(150)*n(9)*n(42);
    %Chemin C2(d)
    dndt(10)= k(150)*n(9)*n(42)-k(153)*n(10)*n(42)-k(154)*n(10)^(2)-k(155)*n(10)*n(1);
    %Chemin C2+:
    dndt(11)= k(61)*n(2)*n(1)+k(78)*n(5)*n(9)+k(87)*n(8)*n(9)...
    +k(92)*n(9)*n(23)+k(93)*n(9)*n(27)-k(33)*n(4)*n(11)...
    -k(45)*n(7)*n(11)-k(50)*n(11)*n(15)-k(59)*n(11)*n(42)...
    -k(74)*n(1)*n(11)-k(94)*n(11)*n(24)-k(122)*n(11)*n(43);
    %Chemin C2-:
    dndt(12)= -k(32)*n(4)*n(12)-k(44)*n(7)*n(12)-k(111)*n(8)*n(12)...
    -k(115)*n(12)*n(25)-k(123)*n(12)*n(43)-k(137)*n(1)*n(12)...
    -k(144)*n(4)*n(12)-k(148)*n(9)*n(12)-k(149)*n(12)*n(13)...
    +k(31)*n(4)*n(14)+k(42)*n(7)*n(14)+k(134)*n(9)*n(42);
    %Chemin C3:
    dndt(13)= k(48)*n(9)*n(9)+k(66)*n(1)*n(9)+k(113)*n(8)*n(14)...
    +k(114)*n(14)*n(25)+k(125)*n(14)*n(43)...
    +k(137)*n(1)*n(12)-k(27)*n(4)*n(13)-k(41)*n(7)*n(13)...
    -k(64)*n(1)*n(13)-k(124)*n(13)*n(43)...
    -k(135)*n(13)*n(42)-k(149)*n(12)*n(13);
    %Chemin C3-:
    dndt(14)= k(135)*n(13)*n(42)-k(31)*n(4)*n(14)...
    -k(42)*n(7)*n(14)-k(113)*n(8)*n(14)...
    -k(114)*n(14)*n(25)-k(125)*n(14)*n(43)-k(138)*n(1)*n(14);
    %Chemin O2:
    dndt(15)= k(17)*n(5)*n(32)+k(19)*n(4)*n(24)+k(20)*n(4)*n(32)...
    +k(21)*n(4)*n(30)+k(57)*n(24)*n(24)+k(70)*n(4)*n(4)...
    +k(77)*n(1)*n(19)+k(97)*n(19)*n(32)+k(98)*n(19)*n(24)...
    +k(108)*n(5)*n(18)+k(110)*n(8)*n(18)+k(116)*n(18)*n(25)...
    +k(117)*n(18)*n(2)+k(127)*n(18)*n(43)+k(143)*n(4)*n(6)...
    +2*k(156)*n(18)*n(15)+2*k(157)*n(18)*n(16)...
    +2*k(158)*n(18)*n(17)+k(161)*n(42)*n(4)*n(15)...
    -k(5)*n(2)*n(15)-k(8)*n(1)*n(15)-k(34)*n(8)*n(15)...
    -k(35)*n(7)*n(15)-k(47)*n(9)*n(15)-k(50)*n(11)*n(15)...
    -k(51)*n(15)*n(37)-k(52)*n(15)*n(24)-k(53)*n(15)*n(20)...
    -k(54)*n(15)*n(27)-k(55)*n(15)*n(22)-k(56)*n(15)*n(26)...
    -k(81)*n(5)*n(15)-k(84)*n(6)*n(15)-k(86)*n(8)*n(15)...
    -k(95)*n(15)*n(21)-k(96)*n(15)*n(23)-k(126)*n(15)*n(43)...
    -k(141)*n(3)*n(15)-k(151)*n(15)*n(42)-k(152)*n(15)*n(42)...
    -k(156)*n(18)*n(15)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15)-k(163)*n(1)*n(15)*n(4);
    %Chemin O2(a):
    dndt(16)= k(151)*n(42)*n(15)-k(157)*n(16)*n(18);
    %Chemin O2(b):
    dndt(17)= k(152)*n(42)*n(15)-k(158)*n(17)*n(18);
    %Chemin O2-:
    dndt(18)= k(84)*n(6)*n(15)+k(162)*n(42)*n(4)*n(15)...
    -k(108)*n(5)*n(18)-k(110)*n(8)*n(18)...
    -k(116)*n(18)*n(25)-k(117)*n(18)*n(2)-k(127)*n(18)*n(43)...
    -k(156)*n(18)*n(15)-k(157)*n(18)*n(16)-k(158)*n(18)*n(17);
    %Chemin O2+:
    dndt(19)= k(81)*n(5)*n(15)+k(86)*n(8)*n(15)+k(95)*n(15)*n(21)...
    +k(96)*n(15)*n(23)-k(36)*n(7)*n(19)-k(49)*n(9)*n(19)-k(77)*n(1)*n(19)...
    -k(97)*n(19)*n(32)-k(98)*n(19)*n(24)-k(128)*n(19)*n(43);
    %Chemin N2:
    dndt(20)= k(21)*n(4)*n(30)+k(37)*n(7)*n(24)+k(39)*n(7)*n(26)...
    +k(57)*n(24)*n(24)+k(76)*n(1)*n(21)+k(82)*n(4)*n(21)...
    +k(91)*n(7)*n(21)+k(95)*n(15)*n(21)+k(99)*n(21)*n(26)...
    +k(100)*n(21)*n(22)-k(11)*n(1)*n(20)-k(18)*n(5)*n(20)...
    -k(23)*n(4)*n(20)-k(53)*n(15)*n(20)...
    -k(129)*n(20)*n(43)-k(159)*n(20)*n(42)-k(160)*n(20)*n(42);
    %Chemin N2+:
    dndt(21)= k(40)*n(7)*n(27)+k(71)*n(8)*n(7)+k(160)*n(20)*n(42)...
    -k(30)*n(4)*n(21)-k(60)*n(21)*n(42)-k(76)*n(1)*n(21)...
    -k(82)*n(4)*n(21)-k(91)*n(7)*n(21)...
    -k(95)*n(15)*n(21)-k(99)*n(21)*n(26)-k(100)*n(21)*n(22);
    %Chemin CO:
    dndt(22)= k(4)*n(1)*n(37)+k(5)*n(2)*n(15)+k(8)*n(1)*n(15)...
    +k(24)*n(4)*n(36)+k(26)*n(4)*n(26)+k(27)*n(4)*n(13)...
    +k(28)*n(4)*n(35)+k(29)*n(4)*n(9)+k(31)*n(4)*n(14)...
    +k(32)*n(4)*n(12)+2*k(47)*n(9)*n(15)+k(49)*n(9)*n(19)...
    +k(50)*n(11)*n(15)+k(54)*n(15)*n(27)+k(58)*n(24)*n(37)...
    +k(65)*n(1)*n(4)+k(75)*n(1)*n(23)+k(83)*n(4)*n(23)...
    +k(92)*n(9)*n(23)+k(96)*n(15)*n(23)+k(102)*n(24)*n(23)...
    +k(139)*n(1)*n(6)+k(140)*n(3)*n(4)-k(13)*n(1)*n(22)...
    -k(55)*n(15)*n(22)-k(79)*n(5)*n(22)-k(88)*n(8)*n(22)...
    -k(100)*n(21)*n(22)-k(101)*n(22)*n(27)...
    -k(130)*n(22)*n(43)-k(146)*n(6)*n(22);
    %Chemin CO+:
    dndt(23)= k(7)*n(2)*n(36)+k(15)*n(5)*n(9)+k(33)*n(4)*n(11)...
    +k(49)*n(9)*n(19)+k(50)*n(11)*n(15)+k(63)*n(2)*n(4)...
    +k(67)*n(1)*n(5)+k(79)*n(5)*n(22)+k(88)*n(8)*n(22)...
    +k(100)*n(21)*n(22)+k(101)*n(22)*n(27)-k(75)*n(1)*n(23)...
    -k(83)*n(4)*n(23)-k(92)*n(9)*n(23)...
    -k(96)*n(15)*n(23)-k(102)*n(24)*n(23);
    %Chemin NO:
    dndt(24)= k(20)*n(4)*n(32)+k(22)*2*n(4)*n(30)+k(23)*n(4)*n(20)...
    +k(24)*n(4)*n(36)+k(25)*n(4)*n(26)+k(34)*n(8)*n(15)...
    +k(35)*n(7)*n(15)+k(51)*n(15)*n(37)+k(106)*n(3)*n(25)...
    +k(114)*n(14)*n(25)+k(115)*n(12)*n(25)+k(116)*n(18)*n(25)...
    +k(147)*n(6)*n(7)-k(10)*n(1)*n(24)-k(14)*n(3)*n(24)-k(19)*n(4)*n(24)...
    -k(37)*n(7)*n(24)-k(52)*n(15)*n(24)-k(57)*2*n(24)*n(24)...
    -k(58)*n(24)*n(37)-k(72)*n(2)*n(24)-k(90)*n(8)*n(24)...
    -k(94)*n(11)*n(24)-k(98)*n(19)*n(24)-k(102)*n(24)*n(23)...
    -k(103)*n(24)*n(27)-k(131)*n(24)*n(43)-k(145)*n(6)*n(24);
    %Chemin NO+:
    dndt(25)= k(6)*n(2)*n(30)+k(16)*n(5)*n(26)+k(17)*n(5)*n(32)...
    +k(18)*n(5)*n(20)+k(30)*n(4)*n(21)+k(36)*n(7)*n(19)...
    +k(54)*n(15)*n(27)+k(72)*n(2)*n(24)+k(90)*n(8)*n(24)...
    +k(94)*n(11)*n(24)+k(98)*n(19)*n(24)+k(102)*n(24)*n(23)...
    +k(103)*n(24)*n(27)-k(106)*n(3)*n(25)-k(114)*n(14)*n(25)...
    -k(115)*n(12)*n(25)-k(116)*n(18)*n(25);
    %Chemin CN:
    dndt(26)= k(4)*n(1)*n(37)+k(6)*n(2)*n(30)+k(7)*n(2)*n(36)...
    +k(10)*n(1)*n(24)+k(11)*n(1)*n(20)+k(12)*n(1)*n(35)...
    +k(28)*n(4)*n(35)+k(41)*n(7)*n(13)+k(42)*n(7)*n(14)...
    +k(43)*n(7)*n(9)+k(45)*n(7)*n(11)+k(46)*2*n(7)*n(35)...
    +k(69)*n(1)*n(7)+k(73)*n(1)*n(27)+k(93)*n(9)*n(27)...
    +k(101)*n(22)*n(27)+k(103)*n(24)*n(27)...
    +k(112)*n(8)*n(28)+k(142)*n(3)*n(7)...
    +k(164)*n(7)*n(1)*n(1)-k(9)*n(1)*n(26)-k(16)*n(5)*n(26)...
    -k(25)*n(4)*n(26)-k(26)*n(4)*n(26)-k(39)*n(7)*n(26)...
    -k(56)*n(15)*n(26)-k(85)*n(6)*n(26)...
    -k(89)*n(8)*n(26)-k(99)*n(21)*n(26);
    %Chemin CN+:
    dndt(27)= k(62)*n(2)*n(7)+k(89)*n(8)*n(26)+k(99)*n(21)*n(26)...
    -k(40)*n(7)*n(27)-k(54)*n(15)*n(27)...
    -k(73)*n(1)*n(27)-k(93)*n(9)*n(27)...
    -k(101)*n(22)*n(27)-k(103)*n(24)*n(27);
    %Chemin CN-:
    dndt(28)= k(14)*n(3)*n(24)+k(44)*n(7)*n(12)+k(85)*n(6)*n(26)-k(112)*n(8)*n(28);
    %Chemin O3:
    dndt(29)= k(163)*n(1)*n(15)*n(4);
    %Chemin N2O:
    dndt(30)= k(38)*n(7)*n(32)+k(53)*n(15)*n(20)+k(58)*n(24)*n(37)...
    -k(6)*n(2)*n(30)-k(21)*n(4)*n(30)...
    -k(22)*n(4)*n(30)-k(80)*n(5)*n(30);
    %Chemin N2O+:
    dndt(31)= k(80)*n(5)*n(30);
    %Chemin NO2:
    dndt(32)= k(52)*n(15)*n(24)+k(145)*n(6)*n(24)-k(17)*n(5)*n(32)...
    -k(20)*n(4)*n(32)-k(38)*n(7)*n(32)...
    -k(97)*n(19)*n(32)-k(165)*n(1)*n(32)*n(4);
    %Chemin NO2+:
    dndt(33)= k(97)*n(19)*n(32);
    %Chemin CO2:
    dndt(34)= k(51)*n(15)*n(37)+k(55)*n(15)*n(22)...
    +k(141)*n(3)*n(15)+k(146)*n(6)*n(22);
    %Chemin C2N:
    dndt(35)= -k(12)*n(1)*n(35)-k(28)*n(4)*n(35)-k(46)*n(7)*n(35);
    %Chemin CNO:
    dndt(36)= -k(7)*n(2)*n(36)-k(24)*n(4)*n(36);
    %Chemin OCN
    dndt(37)= -k(4)*n(1)*n(37)-k(51)*n(15)*n(37)...
    -k(58)*n(24)*n(37)+k(56)*n(15)*n(26);
    %Chemin NO3:
    dndt(38)= k(165)*n(1)*n(32)*n(4);
    %Chemin C2O:
    dndt(39)= k(144)*n(4)*n(12);
    %Chemin C4:
    dndt(40)= k(64)*n(1)*n(13)+k(138)*n(1)*n(14)+k(148)*n(9)*n(12);
    %Chemin C5:
    dndt(41)= k(149)*n(12)*n(13);
    %Chemin e:
    dndt(42)= k(118)*n(1)*n(43)+k(119)*n(3)*n(43)+k(120)*n(6)*n(43)...
    +k(121)*n(12)*n(43)+k(123)*n(12)*n(43)+k(125)*n(14)*n(43)...
    +k(127)*n(18)*n(43)+k(136)*n(1)*n(3)+k(137)*n(1)*n(12)...
    +k(138)*n(1)*n(14)+k(139)*n(1)*n(6)+k(140)*n(3)*n(4)...
    +k(141)*n(3)*n(15)+k(142)*n(3)*n(7)+k(143)*n(4)*n(6)...
    +k(144)*n(4)*n(12)+k(145)*n(6)*n(24)+k(146)*n(6)*n(22)...
    +k(147)*n(6)*n(7)+k(148)*n(9)*n(12)+k(149)*n(12)*n(13)...
    +k(150)*n(9)*n(42)+k(151)*n(15)*n(42)+k(152)*n(15)*n(42)...
    +k(153)*n(10)*n(42)+k(156)*n(18)*n(15)+k(157)*n(18)*n(16)...
    +k(158)*n(18)*n(17)+k(159)*n(20)*n(42)...
    +2*k(160)*n(20)*n(42)-k(1)*n(2)*n(42)-k(2)*n(5)*n(42)...
    -k(3)*n(8)*n(42)-k(59)*n(11)*n(42)-k(60)*n(21)*n(42)...
    -k(132)*n(1)*n(42)-k(133)*n(4)*n(42)-k(134)*n(9)*n(42)...
    -k(135)*n(13)*n(42)-k(150)*n(9)*n(42)-k(151)*n(15)*n(42)...
    -k(152)*n(15)*n(42)-k(153)*n(10)*n(42)-k(159)*n(20)*n(42)...
    -k(160)*n(20)*n(42)-k(161)*n(42)*n(4)*n(15)...
    -k(162)*n(42)*n(4)*n(15);
    %Chemin hv:
    dndt(43)= -k(118)*n(1)*n(43)-k(119)*n(3)*n(43)-k(120)*n(6)*n(43)...
    -k(121)*n(9)*n(43)-k(122)*n(11)*n(43)-k(123)*n(12)*n(43)...
    -k(124)*n(13)*n(43)-k(125)*n(14)*n(43)-k(126)*n(15)*n(43)...
    -k(127)*n(18)*n(43)-k(128)*n(19)*n(43)-k(129)*n(20)*n(43)...
    -k(130)*n(22)*n(43)-k(131)*n(24)*n(43)+k(1)*n(2)*n(42)...
    +k(2)*n(5)*n(42)+k(3)*n(8)*n(42)+k(61)*n(2)*n(1)...
    +k(62)*n(2)*n(7)+k(63)*n(2)*n(4)+k(64)*n(1)*n(13)...
    +k(65)*n(1)*n(4)+k(66)*n(1)*n(9)+k(67)*n(1)*n(5)...
    +k(68)*n(1)*n(1)+k(69)*n(1)*n(7)+k(70)*n(4)*n(4)...
    +k(71)*n(8)*n(7)+k(132)*n(1)*n(42)+k(133)*n(4)*n(42)...
    +k(134)*n(9)*n(42)+k(135)*n(13)*n(42);
end 

.

[Imene Yed]

Thank you so much you're the best !!!!!!!!!!!you helped me a lot with my thesis

[Star Strider]

As always, my pleasure!

Thank you! I know you will do well in your thesis defence!

[Imene Yed]

thank you it's so much kind !!!

[Star Strider]

As always, my pleasure!

[Imene Yed]

@Star Strider I need your help please ,I tried to display the values of n and t in a table because I need it in my thesis

 tn_table = array2table([t n], 'VariableNames',{'time', 'C', 'C+','C-'})
writetable(tn_table, 'c_file'.xlsx')  

I got this error

Error: File: call_myequations5.m Line: 57 Column: 30

Invalid use of operator.

[Star Strider]

If this line is throwing the error, remove the extra single quote:

writetable(tn_table, 'c_file'.xlsx')  
                            ↑ ← HERE

so the correct writetable call will be:

writetable(tn_table, 'c_file.xlsx')  

That should work.

[Imene Yed]

thank you so much it does work thank you so much.

[Star Strider]

As always, my pleasure!

[Imene Yed]

@Star Strider can you help me please ,I'm trying to call a script in another ,script(call_myequation5) is where I have calculated my n values,I want to calculate some properiety in another script called pression.m

I called the script call_myequation5 but it doesn't seem to recognize it.

call_myequation5;

%Calcul pression:
 calcul_Pression=15000*8.3144621*1e-6*(n(:,1)+n(:,2)+n(:,3)+n(:,4)+n(:,5)...
     +n(:,6)+n(:,7)+n(:,8)+n(:,9)+n(:,10)+n(:,11)+n(:,12)+n(:,13)+n(:,14)...
     +n(:,15))
 pressure=[t,calcul_Pression]
 plot(t,calcul_Pression)

[Star Strider]

I am slightly lost as to what you want to do.

The easiest way to get information from a script is to create the script as a function, and then return the necessary values it calculates as the output of the function. There are other ways (writing results to a file in one script, and then reading them in another scritp), however writing the scripts as functions instead is the easiest way to do this.

.

[Walter Roberson]

When you define a variable inside a script that you invoke, and then you use the variable in the outer script (or function), then there are situations in which MATLAB is permitted to ignore the variable you defined in the invoked script.

You have two choices:

1. Use functions instead (more robust!)
2. Before invoking a script, initialize all variables that are assigned to in the script and whose value you will need, as well as all variables that are assigned to in the script and whose value you do not need but which happen to be named the same as a MATLAB function. For example if the invoked script happens to write to a variable named min then even though you think you do not use that variable, assign to min before you call the script, or else you can end up with behaviour that is difficult to predict.

... You might have noticed that using functions is easier.

[Imene Yed]

I have fixed the problem thank you so much for your help@Walter Roberson ,@Star Strider

[Star Strider]

As always, our pleasure!