Initial conditions problem ODE15S

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

Answer 1

[Star Strider]

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

.

[Imene Yed]

@Star Strider it's working the problem is that I have to give values IC to my species in order to have a graph that I can study and dicuss .

[Imene Yed]

this graph doesnt have a physical meaning I need to see at least most of the species

[Star Strider]

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.

To identify them, use the legend function, speciffically something like:

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]

thnk you but I have to set first my initial conditions I only have that of C IC1 ,I need to get a graph first that Ihas a physical meaning (I can see the evolution of some species ).@Star Strider I dont know what values should I give to the other species in order to do my kinetic modelling since I couldnt do experiments because of corona virus .Can you help me please?

[Imene Yed]

for exmeple I need to see a behaviour that look like that

knowing that I'm working with densities (cm-3)

[Star Strider]

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]

Thank you so much

[Star Strider]

My pleasure!