Error using call_myequations Too many input arguments.

1 Jun 2021
1 Answer
Answer Accepted
Updated 5 Jul 2021
3 Views (30 days)

0 votes

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);

Sign in to answer this question.

 Accepted Answer

Star Strider
Star Strider on 1 Jun 2021
Edited: Star Strider on 1 Jun 2021
Ran in:

1 vote

Ran in:
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
.

40 Comments
Show 38 older comments Hide 38 older comments

Imene Yed
Imene Yed on 1 Jun 2021
thnk you it's working but I dont know why it still running and didnt get me results.
Star Strider
Star Strider on 1 Jun 2021
The system is ‘stiff’.
I changed the solver to ode15s in my edit just now, and it works.
Imene Yed
Imene Yed on 1 Jun 2021
thnk you so much for your help.
Imene Yed
Imene Yed on 1 Jun 2021
@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
Star Strider on 1 Jun 2021
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
Imene Yed on 1 Jun 2021
@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
Star Strider on 1 Jun 2021
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
Imene Yed on 1 Jun 2021
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
Star Strider on 1 Jun 2021
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
Imene Yed on 1 Jun 2021
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
Star Strider on 1 Jun 2021
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
Imene Yed on 1 Jun 2021
I still cant see the res in my workspace .I dont understand why there isnt any value of res thnk you
Star Strider
Star Strider on 1 Jun 2021
Try it without the ending semicolon:
res = func(t,n)
the columns of ‘res’ should appear in the Command Window.
Imene Yed
Imene Yed on 1 Jun 2021
yes I saw it ,is it possible to make it in an excel file because I need this data
Star Strider
Star Strider on 1 Jun 2021
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
Imene Yed on 1 Jun 2021
thnk you so much !!!!!!!!!!!!!!!!
Imene Yed
Imene Yed on 2 Jun 2021
@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
Star Strider on 2 Jun 2021
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
Imene Yed on 2 Jun 2021
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
Walter Roberson on 2 Jun 2021
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
Star Strider on 2 Jun 2021
@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
Imene Yed on 3 Jun 2021
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
Star Strider on 3 Jun 2021
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
Imene Yed on 9 Jun 2021
@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
Star Strider on 9 Jun 2021
Ran in:
Ran in:
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
Imene Yed on 9 Jun 2021
Thank you so much you're the best !!!!!!!!!!!you helped me a lot with my thesis
Star Strider
Star Strider on 10 Jun 2021
As always, my pleasure!
Thank you! I know you will do well in your thesis defence!
Imene Yed
Imene Yed on 10 Jun 2021
thank you it's so much kind !!!
Imene Yed
Imene Yed on 1 Jul 2021
@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
Star Strider on 1 Jul 2021
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
Imene Yed on 1 Jul 2021
Edited: Imene Yed on 1 Jul 2021
thank you so much it does work thank you so much.
Imene Yed
Imene Yed on 5 Jul 2021
@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
Star Strider on 5 Jul 2021
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
Walter Roberson on 5 Jul 2021
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
Imene Yed on 5 Jul 2021
I have fixed the problem thank you so much for your help@Walter Roberson ,@Star Strider

Sign in to comment.

More Answers (0)

Categories

Find more on Surrogate Optimization in Help Center and File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!