Initial conditions problem ODE15S

2 Ansichten (letzte 30 Tage)
Imene Yed
Imene Yed am 10 Jun. 2021
Kommentiert: Star Strider am 10 Jun. 2021
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

Akzeptierte Antwort

Star Strider
Star Strider am 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 Kommentare
Imene Yed
Imene Yed am 10 Jun. 2021
Thank you so much
Star Strider
Star Strider am 10 Jun. 2021
My pleasure!

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Get Started with Curve Fitting Toolbox finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by