# Error using call_myequations Too many input arguments.

2 views (last 30 days)
Imene Yed on 1 Jun 2021
Commented: Star Strider on 5 Jul 2021
Hi I divided my script into 4 script:
parameters.m
function k= parameters
%constantes :::
%temperature
T=17000;
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);
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);
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);

Star Strider on 1 Jun 2021
Edited: Star Strider on 1 Jun 2021
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;
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);
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);
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
.
Star Strider on 5 Jul 2021
As always, our pleasure!

### Categories

Find more on Systems of Nonlinear Equations in Help Center and File Exchange

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by