warning: Solving was not successful. The iterative integration loop exited at time t = 0.000000 before the endpoint at tend = 0.001000 was reached

8 views (last 30 days)
Imene Yed
Imene Yed on 28 May 2021
Commented: Imene Yed on 29 May 2021
Hi,I'm trying to solve a big system of ODE here is my code with ODE45
k1 k2 k3...
F=@(t,n)([k1*n(2)*n(42)+k15*n(5)*n(9)+k16*n(5)*n(26)+k25*n(4)*n(26)...
.+k29*n(4)*n(9)+k33*n(4)*n(11)+k39*n(7)*n(26)+k40*n(7)*n(27)...
.+k43*n(7)*n(9)+k44*n(7)*n(12)+k48*n(9)*n(9)+k59*n(11)*n(42)...
.+k72*n(2)*n(24)+k104*n(3)*n(5)+k105*n(3)*n(8)+k106*n(3)*n(25)...
.+k117*n(18)*n(2)+k119*n(3)*n(43)+k122*n(11)*n(43)+k124*n(13)*n(43)...
.+k130*n(22)*n(43)-k4*n(1)*n(37)-k8*n(1)*n(15)-k9*n(1)*n(26)...
.-k10*n(1)*n(24)-k11*n(1)*n(20)-k12*n(1)*n(35)-k13*n(1)*n(22)...
.-k61*n(2)*n(1)-k64*n(1)*n(13)-k65*n(1)*n(4)-k66*n(1)*n(9)...
.-k67*n(1)*n(5)-k68*n(1)*n(1)-k69*n(1)*n(7)-k73*n(1)*n(27)...
.-k74*n(1)*n(11)-k75*n(1)*n(23)-k76*n(1)*n(21)-k77*n(1)*n(19)...
.-k118*n(1)*n(43)-k132*n(1)*n(42)-k136*n(1)*n(3)-k137*n(1)*n(12)...
.-k138*n(1)*n(14)-k139*n(1)*n(6)-2*k164*n(7)*n(1)*n(1)...
.+k164*n(7)*n(1)^(2)-k163*n(1)*n(15)*n(4)+k163*n(1)*n(15)*n(4)...
.-k165*n(1)*n(32)*n(4)+k165*n(1)*n(32)*n(4);%Chemin réactionnel C
-k1*n(2)*n(42)-k5*n(2)*n(15)-k6*n(2)*n(30)-k7*n(2)*n(36)...
.-k61*n(2)*n(1)-k62*n(2)*n(7)-k63*n(2)*n(4)-k72*n(2)*n(24)...
.-k117*n(18)*n(2)+k45*n(7)*n(11)+k73*n(1)*n(27)...
.+k74*n(1)*n(11)+k75*n(1)*n(23)+k76*n(1)*n(21)...
.+k77*n(1)*n(19)+k118*n(1)*n(43)+k122*n(11)*n(43);%chemin reactionnel C+
k32*n(4)*n(12)+k132*n(1)*n(42)-k14*n(3)*n(24)-k104*n(3)*n(5)...
.-k105*n(3)*n(8)-k106*n(3)*n(25)-k119*n(3)*n(43)...
.-k136*n(1)*n(3)-k140*n(3)*n(4)-k141*n(3)*n(15)-k142*n(3)*n(7); %Chemin réactionnel C-
k2*n(5)*n(42)+k8*n(1)*n(15)+k10*n(1)*n(24)+k13*n(1)*n(22)...
.+k14*n(3)*n(24)+k35*n(7)*n(15)+k36*n(7)*n(19)...
.+k37*n(7)*n(24)+k38*n(7)*n(32)+k52*n(15)*n(24)...
.+k53*n(15)*n(20)+k55*n(15)*n(22)+k56*n(15)*n(26)...
.+k78*n(5)*n(9)+k79*n(5)*n(22)+k80*n(5)*n(30)...
.+k81*n(5)*n(15)+k84*n(6)*n(15)+k85*n(6)*n(26)...
.+k104*n(3)*n(5)+2*k107*n(5)*n(6)+k108*n(5)*n(18)...
.+k109*n(6)*n(8)+k120*n(6)*n(43)+2*k126*n(15)*n(43)...
.+k128*n(19)*n(43)+k130*n(22)*n(43)+k131*n(24)*n(43)...
.+k162*n(42)*n(4)*n(15)-k19*n(4)*n(24)-k20*n(4)*n(32)...
.-k21*n(4)*n(30)-k22*n(4)*n(30)-k23*n(4)*n(20)...
.-k24*n(4)*n(36)-k25*n(4)*n(26)-k26*n(4)*n(26)...
.-k27*n(4)*n(13)-k28*n(4)*n(35)-k29*n(4)*n(9)...
.-k30*n(4)*n(21)-k31*n(4)*n(14)-k32*n(4)*n(12)...
.-k33*n(4)*n(11)-k63*n(2)*n(4)-k65*n(1)*n(4)...
.-2*k70*n(4)*n(4)-k82*n(4)*n(21)-k83*n(4)*n(23)...
.-k133*n(4)*n(42)-k140*n(3)*n(4)-k143*n(4)*n(6)...
.-k144*n(4)*n(12)-k161*n(42)*n(4)*n(15)...
.-k162*n(42)*n(4)*n(15)-k163*n(1)*n(15)*n(4)-k165*n(1)*n(32)*n(4);%Chemin réactionnel O
k5*n(2)*n(15)+k34*n(8)*n(15)+k82*n(4)*n(21)...
.+k83*n(4)*n(23)+k128*n(19)*n(43)-k19*n(4)*n(24)...
.-k20*n(4)*n(32)-k21*n(4)*n(30)-k22*n(4)*n(30)...
.-k23*n(4)*n(20)-k24*n(4)*n(36)-k25*n(4)*n(26)...
.-k26*n(4)*n(26)-k27*n(4)*n(13)-k28*n(4)*n(35)...
.-k29*n(4)*n(9)-k30*n(4)*n(21)-k31*n(4)*n(14)...
.-k32*n(4)*n(12)-k33*n(4)*n(11)-k63*n(2)*n(4)...
.-k65*n(1)*n(4)-2*k70*n(4)*n(4)-k82*n(4)*n(21)...
.-k83*n(4)*n(23)-k133*n(4)*n(42)-k140*n(3)*n(4)...
.-k143*n(4)*n(6)-k144*n(4)*n(12)-k161*n(42)*n(4)*n(15)...
.-k162*n(42)*n(4)*n(15)-k163*n(1)*n(15)*n(4)-k165*n(1)*n(32)*n(4);%Chemin réactionnel O+
k133*n(4)*n(42)+k161*n(42)*n(4)*n(15)-k84*n(6)*n(15)...
.-k85*n(6)*n(26)-k107*n(5)*n(6)-k109*n(6)*n(8)-k120*n(6)*n(43)...
.-k139*n(1)*n(6)-k143*n(4)*n(6)-k145*n(6)*n(24)...
.-k146*n(6)*n(22)-k147*n(6)*n(7);% %Chemin réactionnel O-
k3*n(8)*n(42)+k9*n(1)*n(26)+k11*n(1)*n(20)...
.+k18*n(5)*n(20)+k19*n(4)*n(24)+k23*n(4)*n(20)...
.+k26*n(4)*n(26)+k30*n(4)*n(21)+2*k60*n(21)*n(42)...
.+k86*n(8)*n(15)+k87*n(8)*n(9)+k88*n(8)*n(22)...
.+k89*n(8)*n(26)+k90*n(8)*n(24)+k105*n(3)*n(8)...
.+k109*n(6)*n(8)+k110*n(8)*n(18)+k111*n(8)*n(12)...
.+k112*n(8)*n(28)+k113*n(8)*n(14)+2*k129*n(20)*n(43)...
.+k131*n(24)*n(43)+2*k159*n(20)*n(42)-k35*n(7)*n(15)...
.-k36*n(7)*n(19)-k37*n(7)*n(24)-k38*n(7)*n(32)...
.-k39*n(7)*n(26)-k40*n(7)*n(27)-k41*n(7)*n(13)...
.-k42*n(7)*n(14)-k43*n(7)*n(9)-k44*n(7)*n(12)...
.-k45*n(7)*n(11)-k46*n(7)*n(35)-k62*n(2)*n(7)...
.-k69*n(1)*n(7)-k71*n(8)*n(7)-k91*n(7)*n(21)...
.-k142*n(3)*n(7)-k147*n(6)*n(7)-k164*n(7)*n(1)*n(1);%Chemin réactionnel N
k91*n(7)*n(21)-k3*n(8)*n(42)-k34*n(8)*n(15)...
.-k71*n(8)*n(7)-k86*n(8)*n(15)-k87*n(8)*n(9)...
.-k88*n(8)*n(22)-k89*n(8)*n(26)-k90*n(8)*n(24)...
.-k105*n(3)*n(8)-k109*n(6)*n(8)-k110*n(8)*n(18)...
.-k111*n(8)*n(12)-k112*n(8)*n(28)-k113*n(8)*n(14);%Chemin réactionnel N+
k9*n(1)*n(26)+k12*n(1)*n(35)+k13*n(1)*n(22)...
.+k27*n(4)*n(13)+k41*n(7)*n(13)+k68*n(1)*n(1)...
.+k74*n(1)*n(11)+k94*n(11)*n(24)+k111*n(8)*n(12)...
.+k115*n(12)*n(25)+k121*n(9)*n(43)+k123*n(12)*n(43)...
.+k124*n(13)*n(43)+k136*n(1)*n(3)+k153*n(10)*n(42)...
.+k154*n(10)^(2)+k155*n(10)*n(1)-k15*n(5)*n(9)...
.-k29*n(4)*n(9)-k43*n(7)*n(9)-k47*n(9)*n(15)...
.-2*k48*n(9)*n(9)-k49*n(9)*n(19)-k66*n(1)*n(9)...
.-k78*n(5)*n(9)-k87*n(8)*n(9)-k92*n(9)*n(23)...
.-k93*n(9)*n(27)-k121*n(9)*n(43)-k134*n(9)*n(42)...
.-k148*n(9)*n(12)-k150*n(9)*n(42);%Chemin réactionnel C2
k150*n(9)*n(42)-k153*n(10)*n(42)-k154*n(10)^(2)-k155*n(10)*n(1);%C2(d)
k61*n(2)*n(1)+k78*n(5)*n(9)+k87*n(8)*n(9)...
.+k92*n(9)*n(23)+k93*n(9)*n(27)-k33*n(4)*n(11)...
.-k45*n(7)*n(11)-k50*n(11)*n(15)-k59*n(11)*n(42)...
.-k74*n(1)*n(11)-k94*n(11)*n(24)-k122*n(11)*n(43);%Chemin réactionnel C2+
-k32*n(4)*n(12)-k44*n(7)*n(12)-k111*n(8)*n(12)...
.-k115*n(12)*n(25)-k123*n(12)*n(43)-k137*n(1)*n(12)...
.-k144*n(4)*n(12)-k148*n(9)*n(12)-k149*n(12)*n(13)...
.+k31*n(4)*n(14)+k42*n(7)*n(14)+k134*n(9)*n(42);%Chemin réactionnel C2-
k48*n(9)*n(9)+k66*n(1)*n(9)+k113*n(8)*n(14)...
.+k114*n(14)*n(25)+k125*n(14)*n(43)...
.+k137*n(1)*n(12)-k27*n(4)*n(13)-k41*n(7)*n(13)...
.-k64*n(1)*n(13)-k124*n(13)*n(43)...
.-k135*n(13)*n(42)-k149*n(12)*n(13);%Chemin réactionnel C3
k135*n(13)*n(42)-k31*n(4)*n(14)...
.-k42*n(7)*n(14)-k113*n(8)*n(14)...
.-k114*n(14)*n(25)-k125*n(14)*n(43)-k138*n(1)*n(14);%Chemin réactionnel C3-
k17*n(5)*n(32)+k19*n(4)*n(24)+k20*n(4)*n(32)...
.+k21*n(4)*n(30)+k57*n(24)*n(24)+k70*n(4)*n(4)...
.+k77*n(1)*n(19)+k97*n(19)*n(32)+k98*n(19)*n(24)...
.+k108*n(5)*n(18)+k110*n(8)*n(18)+k116*n(18)*n(25)...
.+k117*n(18)*n(2)+k127*n(18)*n(43)+k143*n(4)*n(6)...
.+2*k156*n(18)*n(15)+2*k157*n(18)*n(16)...
.+2*k158*n(18)*n(17)+k161*n(42)*n(4)*n(15)...
.-k5*n(2)*n(15)-k8*n(1)*n(15)-k34*n(8)*n(15)...
.-k35*n(7)*n(15)-k47*n(9)*n(15)-k50*n(11)*n(15)...
.-k51*n(15)*n(37)-k52*n(15)*n(24)-k53*n(15)*n(20)...
.-k54*n(15)*n(27)-k55*n(15)*n(22)-k56*n(15)*n(26)...
.-k81*n(5)*n(15)-k84*n(6)*n(15)-k86*n(8)*n(15)...
.-k95*n(15)*n(21)-k96*n(15)*n(23)-k126*n(15)*n(43)...
.-k141*n(3)*n(15)-k151*n(15)*n(42)-k152*n(15)*n(42)...
.-k156*n(18)*n(15)-k161*n(42)*n(4)*n(15)...
.-k162*n(42)*n(4)*n(15)-k163*n(1)*n(15)*n(4);%Chemin réactionnel O2
k151*n(42)*n(15)-k157*n(16)*n(18);%Chemin réactionnel O2(a)
k152*n(42)*n(15)-k158*n(17)*n(18);%Chemin réactionnel O2(b)
k84*n(6)*n(15)+k162*n(42)*n(4)*n(15)...
.-k108*n(5)*n(18)-k110*n(8)*n(18)...
.-k116*n(18)*n(25)-k117*n(18)*n(2)-k127*n(18)*n(43)...
.-k156*n(18)*n(15)-k157*n(18)*n(16)-k158*n(18)*n(17);%Chemin réactionnel O2-
k81*n(5)*n(15)+k86*n(8)*n(15)+k95*n(15)*n(21)...
.+k96*n(15)*n(23)-k36*n(7)*n(19)-k49*n(9)*n(19)-k77*n(1)*n(19)...
.-k97*n(19)*n(32)-k98*n(19)*n(24)-k128*n(19)*n(43);%Chemin réactionnel O2+
k21*n(4)*n(30)+k37*n(7)*n(24)+k39*n(7)*n(26)...
.+k57*n(24)*n(24)+k76*n(1)*n(21)+k82*n(4)*n(21)...
.+k91*n(7)*n(21)+k95*n(15)*n(21)+k99*n(21)*n(26)...
.+k100*n(21)*n(22)-k11*n(1)*n(20)-k18*n(5)*n(20)...
.-k23*n(4)*n(20)-k53*n(15)*n(20)...
.-k129*n(20)*n(43)-k159*n(20)*n(42)-k160*n(20)*n(42);%Chemin réactionnel N2
k40*n(7)*n(27)+k71*n(8)*n(7)+k160*n(20)*n(42)...
.-k30*n(4)*n(21)-k60*n(21)*n(42)-k76*n(1)*n(21)...
.-k82*n(4)*n(21)-k91*n(7)*n(21)...
.-k95*n(15)*n(21)-k99*n(21)*n(26)-k100*n(21)*n(22);%Chemin réactionnel N2+
k4*n(1)*n(37)+k5*n(2)*n(15)+k8*n(1)*n(15)...
.+k24*n(4)*n(36)+k26*n(4)*n(26)+k27*n(4)*n(13)...
.+k28*n(4)*n(35)+k29*n(4)*n(9)+k31*n(4)*n(14)...
.+k32*n(4)*n(12)+2*k47*n(9)*n(15)+k49*n(9)*n(19)...
.+k50*n(11)*n(15)+k54*n(15)*n(27)+k58*n(24)*n(37)...
.+k65*n(1)*n(4)+k75*n(1)*n(23)+k83*n(4)*n(23)...
.+k92*n(9)*n(23)+k96*n(15)*n(23)+k102*n(24)*n(23)...
.+k139*n(1)*n(6)+k140*n(3)*n(4)-k13*n(1)*n(22)...
.-k55*n(15)*n(22)-k79*n(5)*n(22)-k88*n(8)*n(22)...
.-k100*n(21)*n(22)-k101*n(22)*n(27)...
.-k130*n(22)*n(43)-k146*n(6)*n(22);%Chemin réactionnel CO
k7*n(2)*n(36)+k15*n(5)*n(9)+k33*n(4)*n(11)...
.+k49*n(9)*n(19)+k50*n(11)*n(15)+k63*n(2)*n(4)...
.+k67*n(1)*n(5)+k79*n(5)*n(22)+k88*n(8)*n(22)...
.+k100*n(21)*n(22)+k101*n(22)*n(27)-k75*n(1)*n(23)...
.-k83*n(4)*n(23)-k92*n(9)*n(23)...
.-k96*n(15)*n(23)-k102*n(24)*n(23);%Chemin réactionnel CO+
k20*n(4)*n(32)+k22*2*n(4)*n(30)+k23*n(4)*n(20)...
.+k24*n(4)*n(36)+k25*n(4)*n(26)+k34*n(8)*n(15)...
.+k35*n(7)*n(15)+k51*n(15)*n(37)+k106*n(3)*n(25)...
.+k114*n(14)*n(25)+k115*n(12)*n(25)+k116*n(18)*n(25)...
.+k147*n(6)*n(7)-k10*n(1)*n(24)-k14*n(3)*n(24)-k19*n(4)*n(24)...
.-k37*n(7)*n(24)-k52*n(15)*n(24)-k57*2*n(24)*n(24)...
.-k58*n(24)*n(37)-k72*n(2)*n(24)-k90*n(8)*n(24)...
.-k94*n(11)*n(24)-k98*n(19)*n(24)-k102*n(24)*n(23)...
.-k103*n(24)*n(27)-k131*n(24)*n(43)-k145*n(6)*n(24);%Chemin réactionnel NO
k6*n(2)*n(30)+k16*n(5)*n(26)+k17*n(5)*n(32)...
.+k18*n(5)*n(20)+k30*n(4)*n(21)+k36*n(7)*n(19)...
.+k54*n(15)*n(27)+k72*n(2)*n(24)+k90*n(8)*n(24)...
.+k94*n(11)*n(24)+k98*n(19)*n(24)+k102*n(24)*n(23)...
.+k103*n(24)*n(27)-k106*n(3)*n(25)-k114*n(14)*n(25)...
.-k115*n(12)*n(25)-k116*n(18)*n(25);%Chemin réactionnel NO+
k4*n(1)*n(37)+k6*n(2)*n(30)+k7*n(2)*n(36)...
.+k10*n(1)*n(24)+k11*n(1)*n(20)+k12*n(1)*n(35)...
.+k28*n(4)*n(35)+k41*n(7)*n(13)+k42*n(7)*n(14)...
.+k43*n(7)*n(9)+k45*n(7)*n(11)+k46*2*n(7)*n(35)...
.+k69*n(1)*n(7)+k73*n(1)*n(27)+k93*n(9)*n(27)...
.+k101*n(22)*n(27)+k103*n(24)*n(27)...
.+k112*n(8)*n(28)+k142*n(3)*n(7)...
.+k164*n(7)*n(1)*n(1)-k9*n(1)*n(26)-k16*n(5)*n(26)...
.-k25*n(4)*n(26)-k26*n(4)*n(26)-k39*n(7)*n(26)...
.-k56*n(15)*n(26)-k85*n(6)*n(26)...
.-k89*n(8)*n(26)-k99*n(21)*n(26);%Chemin réactionnel CN
k62*n(2)*n(7)+k89*n(8)*n(26)+k99*n(21)*n(26)...
.-k40*n(7)*n(27)-k54*n(15)*n(27)...
.-k73*n(1)*n(27)-k93*n(9)*n(27)...
.-k101*n(22)*n(27)-k103*n(24)*n(27);%Chemin réactionnel CN+
k14*n(3)*n(24)+k44*n(7)*n(12)+k85*n(6)*n(26)-k112*n(8)*n(28);%CN-
k163*n(1)*n(15)*n(4);%Chemin réactionnel O3
k38*n(7)*n(32)+k53*n(15)*n(20)+k58*n(24)*n(37)...
.-k6*n(2)*n(30)-k21*n(4)*n(30)...
.-k22*n(4)*n(30)-k80*n(5)*n(30);%Chemin réactionnel N2O
k80*n(5)*n(30);%Chemin réactionnel N2O+
k52*n(15)*n(24)+k145*n(6)*n(24)-k17*n(5)*n(32)...
.-k20*n(4)*n(32)-k38*n(7)*n(32)...
.-k97*n(19)*n(32)-k165*n(1)*n(32)*n(4);%Chemin réactionnel NO2
k97*n(19)*n(32);%Chemin réactionnel NO2+
k51*n(15)*n(37)+k55*n(15)*n(22)...
.+k141*n(3)*n(15)+k146*n(6)*n(22);%Chemin réactionnel CO2
-k12*n(1)*n(35)-k28*n(4)*n(35)-k46*n(7)*n(35);%Chemin réactionnel C2N
-k7*n(2)*n(36)-k24*n(4)*n(36);%Chemin réactionnel CNO
-k4*n(1)*n(37)-k51*n(15)*n(37)...
.-k58*n(24)*n(37)+k56*n(15)*n(26);%Chemin réactionnel OCN
k165*n(1)*n(32)*n(4);%Chemin réactionnel NO3
k144*n(4)*n(12);%Chemin réactionnel C2O
k64*n(1)*n(13)+k138*n(1)*n(14)+k148*n(9)*n(12);%Chemin réactionnel C4
k149*n(12)*n(13);%Chemin réactionnel C5
k118*n(1)*n(43)+k119*n(3)*n(43)+k120*n(6)*n(43)...
.+k121*n(12)*n(43)+k123*n(12)*n(43)+k125*n(14)*n(43)...
.+k127*n(18)*n(43)+k136*n(1)*n(3)+k137*n(1)*n(12)...
.+k138*n(1)*n(14)+k139*n(1)*n(6)+k140*n(3)*n(4)...
.+k141*n(3)*n(15)+k142*n(3)*n(7)+k143*n(4)*n(6)...
.+k144*n(4)*n(12)+k145*n(6)*n(24)+k146*n(6)*n(22)...
.+k147*n(6)*n(7)+k148*n(9)*n(12)+k149*n(12)*n(13)...
.+k150*n(9)*n(42)+k151*n(15)*n(42)+k152*n(15)*n(42)...
.+k153*n(10)*n(42)+k156*n(18)*n(15)+k157*n(18)*n(16)...
.+k158*n(18)*n(17)+k159*n(20)*n(42)...
.+2*k160*n(20)*n(42)-k1*n(2)*n(42)-k2*n(5)*n(42)...
.-k3*n(8)*n(42)-k59*n(11)*n(42)-k60*n(21)*n(42)...
.-k132*n(1)*n(42)-k133*n(4)*n(42)-k134*n(9)*n(42)...
.-k135*n(13)*n(42)-k150*n(9)*n(42)-k151*n(15)*n(42)...
.-k152*n(15)*n(42)-k153*n(10)*n(42)-k159*n(20)*n(42)...
.-k160*n(20)*n(42)-k161*n(42)*n(4)*n(15)...
.-k162*n(42)*n(4)*n(15);%Chemin réactionnel e-
-k118*n(1)*n(43)-k119*n(3)*n(43)-k120*n(6)*n(43)...
.-k121*n(9)*n(43)-k122*n(11)*n(43)-k123*n(12)*n(43)...
.-k124*n(13)*n(43)-k125*n(14)*n(43)-k126*n(15)*n(43)...
.-k127*n(18)*n(43)-k128*n(19)*n(43)-k129*n(20)*n(43)...
.-k130*n(22)*n(43)-k131*n(24)*n(43)+k1*n(2)*n(42)...
.+k2*n(5)*n(42)+k3*n(8)*n(42)+k61*n(2)*n(1)...
.+k62*n(2)*n(7)+k63*n(2)*n(4)+k64*n(1)*n(13)...
.+k65*n(1)*n(4)+k66*n(1)*n(9)+k67*n(1)*n(5)...
.+k68*n(1)*n(1)+k69*n(1)*n(7)+k70*n(4)*n(4)...
.+k71*n(8)*n(7)+k132*n(1)*n(42)+k133*n(4)*n(42)...
.+k134*n(9)*n(42)+k135*n(13)*n(42)])
OPTs = OPTs = odeset('reltol', 1e-3, 'abstol', 1e-5)
tic,[Time, SOL]=ode45(F,[0 1e-3],IC, OPTs);
toc
plot(Time,SOL)
I got these errors
warning: Solving was not successful. The iterative integration loop exited
at time t = 0.000000 before the endpoint at tend = 0.001000 was reached. Thi
s may happen if the stepsize becomes too small. Try to reduce the value of '
InitialStep' and/or 'MaxStep' with the command 'odeset'.
warning: called from
integrate_adaptive at line 318 column 7
ode45 at line 241 column 12
imene14echkabilov at line 246 column 16

Accepted Answer

Jan
Jan on 28 May 2021
Edited: Jan on 28 May 2021
The first guess in chemical systems is usually: The problem is stiff. Then you need a solver for stiff problems as ODE15S instead of the non-stiff ODE45.
By the way, the code does not run, because there are no ".+" and ".-" operators in Matlab.
The coe would be much easier to debug, if you do not squeeze a huge equation in a function handle, but if you write it to a function. I assume it the the result of some matrix equations. Then the matrix notation would be easier the debug and it even would run faster, most likely.
  7 Comments

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!

Translated by