Initial conditions problem ODE15S
Show older comments
Hi ,here is my code if someone can help me to set initial conditions(densities IC) without changing the time so that I can have a graph with a meaning, I only have the density of IC1 which is 1.04e21 and knowing that I have a system of carbon plasma under air .
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
function [t,n]= call_myequations()
%Initial Conditions
IC1 = 1.04e21; % C
IC2 =0 ; % C+
IC3 =0 ; % C-
IC4 =0 ; % O
IC5 = 0; % O+
IC6 = 0; % O-
IC7 = 0; % N
IC8 = 0; % N+
IC9 = 0; % C2
IC10 = 0; % C2(d)
IC11 = 0; % C2+
IC12 = 0; % C2-
IC13 = 0; % C3
IC14 = 0; % C3-
IC15 = 07; % O2
IC16 = 0; % O2(a)
IC17 = 0; % O2(b)
IC18 = 0; % O2-
IC19 = 0; % O2+
IC20 = 0; % N2
IC21 = 0; % N2+
IC22 =0; % CO
IC23 = 0; % CO+
IC24 = 0; % NO
IC25 = 0; % NO+
IC26 = 0; % CN
IC27 = 0; % CN+
IC28 = 0; % CN-
IC29 = 0; % O3
IC30 = 0; % N2O
IC31 = 0; % N2O+
IC32 = 0; % NO2
IC33 = 0; % NO2+
IC34 = 0; % CO2
IC35 = 0; % C2N
IC36 = 0; % CNO
IC37 = 0; % OCN
IC38 = 0; % NO3
IC39 = 0; % C2O
IC40 = 0; % C4
IC41 = 0; % C5
IC42 = 0; % electron
IC43 = 0; % hv
k= parameters;
;
[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-6],[IC1 IC2 IC3 IC4 IC5 IC6 IC7 IC8 IC9 IC10...
IC11 IC12 IC13 IC14 IC15 IC16 IC17 IC18 IC19 IC20 IC21 IC22 IC23...
IC24 IC25 IC26 IC27 IC28 IC29 IC30 IC31 IC32 IC33 IC34 IC35 IC36 IC37...
IC38 IC39 IC40 IC41 IC42 IC43]);
title('chimie du plasma sous air ')
xlabel('Temps (s)')
ylabel('Densité (cm-3)')
I got this plot which has no meaning
Accepted Answer
Star Strider
on 10 Jun 2021
I will copy and run that (again) if necessary.
In the interim, see if actually calling the plot function works —
figure
plot(t,n)
title('chimie du plasma sous air ')
xlabel('Temps (s)')
ylabel('Densité (cm-3)')
.
8 Comments
Imene Yed
on 10 Jun 2021
Imene Yed
on 10 Jun 2021
this graph doesnt have a physical meaning I need to see at least most of the species
Star Strider
on 10 Jun 2021
That is going to be a bit difficult because of the number of variables (columns of ‘n’).
Providing that none of them are negative, see if plotting them as:
semilogy(t,n)
works to distinguish them.
speciesc = compose('%s','C','C+','C-', ,'electron','hv')
legend(speciesc, 'Location','eastoutside')
where all the comments in the ‘Initial Conditions’ vector are supplied to it as arguments. It will likely require the 'NumColumns' name-value pair (probably set to 4) in order to look acceptable.
.
Imene Yed
on 10 Jun 2021
Scaling the y-axis logarithmically is the best I can do. (I also limited the y-axis and took the absolute values of any elements in ‘n’ that were negative, since negative valuues will not plot on a logarithmic scale. I had to copy ‘parameters’ from you earlier post in order to run this.) This shows the data in detail. Change my code to remove the ‘XScale’ logarithmic scaling to get linear scaling.
[t,n] = call_myequations;
function [t,n]= call_myequations()
%Initial Conditions
IC1 = 1.04e21; % C
IC2 =0 ; % C+
IC3 =0 ; % C-
IC4 =0 ; % O
IC5 = 0; % O+
IC6 = 0; % O-
IC7 = 0; % N
IC8 = 0; % N+
IC9 = 0; % C2
IC10 = 0; % C2(d)
IC11 = 0; % C2+
IC12 = 0; % C2-
IC13 = 0; % C3
IC14 = 0; % C3-
IC15 = 07; % O2
IC16 = 0; % O2(a)
IC17 = 0; % O2(b)
IC18 = 0; % O2-
IC19 = 0; % O2+
IC20 = 0; % N2
IC21 = 0; % N2+
IC22 =0; % CO
IC23 = 0; % CO+
IC24 = 0; % NO
IC25 = 0; % NO+
IC26 = 0; % CN
IC27 = 0; % CN+
IC28 = 0; % CN-
IC29 = 0; % O3
IC30 = 0; % N2O
IC31 = 0; % N2O+
IC32 = 0; % NO2
IC33 = 0; % NO2+
IC34 = 0; % CO2
IC35 = 0; % C2N
IC36 = 0; % CNO
IC37 = 0; % OCN
IC38 = 0; % NO3
IC39 = 0; % C2O
IC40 = 0; % C4
IC41 = 0; % C5
IC42 = 0; % electron
IC43 = 0; % hv
k= parameters;
;
[t,n] = ode15s(@(t,n)myequations(t,n,k),[0 1e-6],[IC1 IC2 IC3 IC4 IC5 IC6 IC7 IC8 IC9 IC10...
IC11 IC12 IC13 IC14 IC15 IC16 IC17 IC18 IC19 IC20 IC21 IC22 IC23...
IC24 IC25 IC26 IC27 IC28 IC29 IC30 IC31 IC32 IC33 IC34 IC35 IC36 IC37...
IC38 IC39 IC40 IC41 IC42 IC43]);
n(n<0) = abs(n(n<0));
figure
plot(t,n)
title('chimie du plasma sous air ')
xlabel('Temps (s)')
ylabel('Densité (cm^{-3})')
grid
xlim([1E-22 1E-6])
ylim([1E-85 1E+25])
Ax = gca;
Ax.XScale = 'log';
Ax.YScale = 'log';
% pos = Ax.Position; % Use These To 'Stretch' The Plot Vertically
% Ax.Position = pos+[0 -200 0 200];
end
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
on 10 Jun 2021
Star Strider
on 10 Jun 2021
My pleasure!
More Answers (0)
Categories
Find more on Analysis and Applications 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!
