calling functions in other scripts

1 Ansicht (letzte 30 Tage)
Imene Yed
Imene Yed am 30 Mai 2021
Kommentiert: Imene Yed am 30 Mai 2021
Hi, I want call this function :
classdef functioncontainer
methods
function res = func(t,n)
%Conservation de la masse:
%C:
masse_C(:)= n(:,1) + 2*n(:,9)+2*n(:,10)+3*n(:,13)+...
n(:,22)+n(:,26)+n(:,34)+2*n(:,35)+n(:,36)...
+n(:,37)+2*n(:,39)+4*n(:,40)+5*n(:,41)
%O:
masse_O(:)= n(:,4)+2*n(:,15)+2*n(:,16)+2*n(:,17)+n(:,22)...
+n(:,24)+ 3*n(:,29)+n(:,30)+2*n(:,32)+2*n(:,34)+...
n(:,36)+n(:,37)+3*n(:,38)+ n(:,39);
%N
masse_N(:)= n(:,7)+2*n(:,20)+n(:,24)+n(:;26)+2*n(:,30)+n(:,32)...
+n(:,35)+n(:,36)+n(:,37)+n(:,38)
in this script
...............
tic,[t,n] = ode45(@(t,n)myequations(t,n,k1,k26,k57,k58,k59...),[0 10e-3],IC)
toc
plot(t,n)
myobj= functioncontainer
res1=myobj.func(t,n)
I don't know what is missing .
Elapsed time is 1.103608 seconds.
Error using functioncontainer
Error: File: functioncontainer.m Line: 15 Column: 44
Invalid expression. When calling a function or indexing a variable, use parentheses. Otherwise, check for
mismatched delimiters.
Error in finalODE (line 204)
myobj= functioncontainer
  7 Kommentare
5 ältere Kommentare anzeigen
Torsten
Torsten am 30 Mai 2021
Just a sidenote: I wonder why you want to use object-oriented programming when you just intend to check some mass balances. Don't you think it's kind of overkill ?
Imene Yed
Imene Yed am 30 Mai 2021
is there something else I can do ?please

Melden Sie sich an, um zu kommentieren.

Melden Sie sich an, um diese Frage zu beantworten.

Antworten (1)

Steven Lord
Steven Lord am 30 Mai 2021
This is not just a function. It is a class method.
classdef functioncontainer
methods
function res = func(t,n)
As written, either t or n must be a functioncontainer object. It's not really clear to me how you're trying to use this (your code excerpt refers to "myequations" which is not defined in that excerpt) as your ODE function.
But I concur with Torsten. In this case using a class seems unnecessary. Instead I would simply write a function like the two examples on this documentation page to evaluate the right hand side of your system of ODEs and call it in a similar manner as those examples do. I would also advise you to extract some or all of the components of the second input to variables with descriptive names at the start of your file as that may make your code look more like the equations in the textbook or paper that you're trying to implement.
function dydt = mysystem(t, y)
% Concentrations
O2conc = y(1);
CO2conc = y(2);
% etc.
% now work with O2conc, CO2conc, etc.
  1 Kommentar
Imene Yed
Imene Yed am 30 Mai 2021
here is my entire code
clear all,clc
%constantes :::
%temperature
T=17000;
%constante k: recombinaison radiative
k1=2.36 *(1e-12)*(T/300)^(-0.29)*exp(17.60/T);
k2 =3.24 *(1e-12) *(T/300)^(-0.66);
k3=3.50 *(1e-12)*(T/300)^(-0.53)* exp(3.20/T);
%constante k::Reconmbinaison
k4=1*10^(-10);
k5=4.54*10^(-10);
k6=9.10*10^(-10);
k7=8.98*10^(-9)*(T/300)^(-0.5);
k8=5.56*10^(-11)*(T/300)^(0.41)*exp(26.90/T);
k9=4.98*10^(-10)*exp(-18116/T);
k10=6*10^(-11)*(T/300)^(-0.16);
k11=8.69*10^(-11)*exp(-22600/T);
k12=1*10^(-10);
k13=2.94*10^(-11)*(T/300)^(0.5)*exp(-58025/T);
k14=1*10^(-9);
k15=4.80*10^(-10)
k16=1.00*10^(-9)*(T/300)^(-0.5);
k17=8.30*10^(-10);
k18=2.42*10^(-12)*(T/300)^(-0.21)*exp(44/T);
k19=1.18*10^(-11)*exp(-20413/T);
k20=9.82*10^(-12)*(T/300)^(-0.21)*exp(-5.20/T);
k21=1.66*10^(-10)*exp(-14100/T);
k22=1.15*10^(-10)*exp(-13400/T);
k23=2.51*10^(-10)*exp(-38602/T);
k24=1.00*10^(-10);
k25=5.37*10^(-11)*exp(-13800/T);
k26=5.00*10^(-11)*exp(-200/T);
k27=5.00*10^(-12)*exp(-900/T);
k28=6.00*10^(-12);
k29=2.00*10^(-10)*(T/300)^(-0.12);
k30=1.30*10^(-10);
k31=5.00*10^(-10);
k32=2.90*10^(-10);
k33=3.10*10^(-10);
k34=3.66*10^(-11);
k35=2.26*10^(-12)*(T/300)^(0.86)*exp(-3134.0/T);
k36=1.80*10^(-10);
k37=3.38*10^(-11)*(T/300)^(-0.17)*exp(2.80/T);
k38=3*10^(-12);
k39=1.00*10^(-10)*(T/300)^(0.4);
k40=6.10*10^(-10);
k41=1*10^(-13);
k42=5*10^(-11);
k43=5*10^(-11);
k44=1.15*10^(-10);
k45=4*10^(-11);
k46=1*10^(-10);
k47=1.50*10^(-11)*exp(-4300/T);
k48=5.3137*10^(-10);
k49=4.10*10^(-10);
k50=8*10^(-10);
k51=1.32*10^(-12);
k52=2.8*10^(-12)*exp(23400/T);
k53=1*10^(-10)*exp(-55200/T);
k54=8.60*10^(-11);
k55=5.99*10^(-12)*exp(-24075/T);
k56=2.02*10^(-11)*(T/300)^(-0.19)*exp(31.90/T);
k57=2.51*10^(-11)*exp(-30653/T);
k58=3.47*10^(-11)*(T/300)^(-1.33)*exp(-242/T);
%constante k::Recombinaison dissociative
k59=3*10^(-7)*(T/300)^(-0.5);
k60=1.70*10^(-7)*(T/300)^(-0.3);
%constante k::Association radiative
k61=4.01*10^(-18)*(T/300)^(0.17)*exp(-101.50/T);
k62=1.08*10^(-18)*(T/300)^(0.07)*exp(-57.50/T);
k63=3.14*10^(-18)*(T/300)^(-0.15)*exp(-68/T);
k64=4.00*10^(-14)*(T/300)^(-1);
k65=4.69*10^(-19)*(T/300)^(1.52)*exp(50.50/T);
k66=3*10^(-16)*(T/300)^(-1);
k67=5*10^(-10)*(T/300)^(-3.7)*exp(-800/T);
k68=4.36*10^(-18)*(T/300)^(0.35)*exp(-161.30/T);
k69=5.72*10^(-19)*(T/300)^(0.37)*exp(-51/T);
k70=4.90*10^(-20)*(T/300)^(1.58);
k71=3.71*10^(-18)*(T/300)^(0.24)*exp(-26.10/T);
%constante k::Transfert de charge
k72=7.05*10^(-10)*(T/300)^(-0.03)*exp(16.70/T);
k73=1*10^(-10);
k74=1.10*10^(-10);
k75=1.10*10^(-10);
k76=1.10*10^(-10);
k77=5.20*10^(-11);
k78=4.80*10^(-10);
k79=4.90*10^(-12)*(T/300)^(0.5)*exp(-4580/T);
k80=6.30*10^(-10);
k81=1.90*10^(-11);
k82=1.00*10^(-11);
k83=1.40*10^(-10);
k84=7.30*10^(-10)*exp(-890/T);
k85=1.00*10^(-9);
k86=3.11*10^(-10);
k87=1.00*10^(-9);
k88=8.25*10^(-10);
k89=1.10*10^(-9)*(T/300)^(-0.50);
k90=4.51*10^(-10);
k91=1* 10^(-11);
k92=8.40*10^(-10);
k93=8.50^(-10);
k94=3.40*10^(-10);
k95=5*10^(-11);
k96=1.20*10^(-10);
k97=6.60*10^(-10);
k98=4.60*10^(-10);
k99=1.00*10^(-10)*(T/300)^(-0.5);
k100=7.40*10^(-11);
k101=6.30*10^(-10);
k102=3.30*10^(-10);
k103=5.70*10^(-10);
%constante k::Neutralisation
k104=7.51*10^(-8)*(T/300)^(-0.5);
k105=7.51*10^(-8)*(T/300)^(-0.5);
k106=7.51*10^(-8)*(T/300)^(-0.5);
k107=7.51*10^(-8)*(T/300)^(-0.5);
k108=7.51*10^(-8)*(T/300)^(-0.5);
k109=7.51*10^(-8)*(T/300)^(-0.5);
k110=7.51*10^(-8)*(T/300)^(-0.50);
k111=7.51*10^(-8)*(T/300)^(-0.50);
k112=7.51*10^(-8)*(T/300)^(-0.50);
k113=7.51*10^(-8)*(T/300)^(-0.50);
k114=7.51*10^(-8)*(T/300)^(-0.5);
k115=7.51*10^(-8)*(T/300)^(-0.5);
k116=7.51*10^(-8)*(T/300)^(-0.5);
k117=7.51*10^(-8)*(T/300)^(-0.5);
%constante k::Photoionisation
k118=3.10*10^(-10)*exp(-3.3/T);
k119=4.90*10^(-8)*exp(-0.50/T);
k120=1.09*10^(-8)*exp(-0.5/T);
k121=4.10*10^(-10)*exp(-3.80/T);
k122=1*10^(-11)*exp(-1.70/T);
k123=2.96*10^(-9)*exp(-2/T);
k124=5*10^(-9)*exp(-2.1/T);
k125=6.88*10^(-9)*exp(-1.5/T);
k126=7.90*10^(-10)*exp(-2.1/T);
k127=6.10*10^(-9)*exp(-0.5/T);
k128=3.50*10^(-11)*exp(-2/T);
k129=2.30*10^(-10)*exp(-3.90/T);
k130=2.10*10^(-10)*exp(-3.50/T);
k131=4.70*10^(-10)*exp(-2.10/T);
%constante k::Attachement radiatif
k132=2.25*10^(-15);
k133=1.50*10^(-15);
k134=2.00*10^(-15)*(T/300)^(-0.5);
k135=1.70*10^(-14)*(T/300)^(-0.5);
%constante k::Détachement associatif
k136=5*10^(-10);
k137=1*10^(-9);
k138=1*10^(-9);
k139=5*10^(-10);
k140=5*10^(-10);
k141=5*10^(-11);
k142=5*10^(-10);
k143=1.90*10^(-10);
k144=2.90*10^(-10);
k145=3.10*10^(-10)*(T/300)^(-0.83);
k146=6.50*10^(-10);
k147=2.20*10^(-10);
k148=1.00*10^(-9);
k149=1*10^(-9);
%constante k::Excitation/ Désexcitation
k150=2.0425*10^(-7)*T^(-0.3)*exp(-29489.55/T);
k151=2.6*10^(-10)*T^(0.89)*exp(-11372.9/T);
k152=1.2*10^(-10)*T^(0.59)*exp(-18916.15/T);
k153=6.0278*10^(-7)*T^(-0.4)*exp(-19880.36/T);
k154=1.6306*10^(-17);
k155=3.6532*10^(-11);
%constante k::Détachement électronique
k156=7*10^(-10);
k157=2*10^(-10);
k158=3.6*10^(-10);
%constante k::Dissociation par impact électronique
k159= 6.2*10^(-10) * T^(1.75) * exp(-113032.7/T );
%constante k::Ionisation par impact électronique
k160= 7.1*10^(-9) * T^(0.76) * exp(-181038/T );
%constante k::Attachement électronique cn(1)6/s
k161=1*10^(-31);
k162=1*10^(-31);
%constante k::Processus de recombinaison cn(1)6/s
k163=5.70*10^(-34)*T^(-2.6);
k164=2.5919*10^(-80);
k165=1.31*10^(-31)*T^(-1.5);
%condition intial
%IC=IC1=IC2=....=IC43=10^6cm^3
IC =10^18*ones(1,43);
%opt=odeset('OutputFcn',@odeplot,'OutputSel',[1:43]);
tic,[t,n] = ode45(@(t,n)myequations(t,n,k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13,k14,k15,k16,k17,k18,k19,k20,k21,k22,k23,k24,k25,k26,k27,k28,k29,k30,k31,k32,k33,k34,k35,k36,k37,k38,k39,k40,k41,k42,k43,k44,k45,k46,k47,k48,k49,k50,k51,k52,k53,k54,k55,k56,k57,k58,k59,k60,k61,k62,k63,k64,k65,k66,k67,k68,k69,k70,k71,k72,k73,k74,k75,k76,k77,k78,k79,k80,k81,k82,k83,k84,k85,k86,k87,k88,k89,k90,k91,k92,k93,k94,k95,k96,k97,k98,k99,k100,k101,k102,k103,k104,k105,k106,k107,k108,k109,k110,k111,k112,k113,k114,k115,k116,k117,k118,k119,k120,k121,k122,k123,k124,k125,k126,k127,k128,k129,k130,k131,k132,k133,k134,k135,k136,k137,k138,k139,k140,k141,k142,k143,k144,k145,k146,k147,k148,k149,k150,k151,k152,k153,k154,k155,k156,k157,k158,k159,k160,k161,k162,k163,k164,k165),[0 10e-3],IC)
toc
myobj= functioncontainer
res1=myobj.func(t,n)
function dndt=myequations(t,n,k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13,k14,k15,k16,k17,k18,k19,k20,k21,k22,k23,k24,k25,k26,k27,k28,k29,k30,k31,k32,k33,k34,k35,k36,k37,k38,k39,k40,k41,k42,k43,k44,k45,k46,k47,k48,k49,k50,k51,k52,k53,k54,k55,k56,k57,k58,k59,k60,k61,k62,k63,k64,k65,k66,k67,k68,k69,k70,k71,k72,k73,k74,k75,k76,k77,k78,k79,k80,k81,k82,k83,k84,k85,k86,k87,k88,k89,k90,k91,k92,k93,k94,k95,k96,k97,k98,k99,k100,k101,k102,k103,k104,k105,k106,k107,k108,k109,k110,k111,k112,k113,k114,k115,k116,k117,k118,k119,k120,k121,k122,k123,k124,k125,k126,k127,k128,k129,k130,k131,k132,k133,k134,k135,k136,k137,k138,k139,k140,k141,k142,k143,k144,k145,k146,k147,k148,k149,k150,k151,k152,k153,k154,k155,k156,k157,k158,k159,k160,k161,k162,k163,k164,k165)
dndt=zeros(43,1)
%Chemin C
dndt(1)= 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 C+:
dndt(2)= -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 C-:
dndt(3)= 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 O:
dndt(4)= 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 O+:
dndt(5)= 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 O-:
dndt(6)= 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 N:
dndt(7)= 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 N+:
dndt(8)=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 C2:
dndt(9)= 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 C2(d)
dndt(10)= k150*n(9)*n(42)-k153*n(10)*n(42)-k154*n(10)^(2)-k155*n(10)*n(1);
%Chemin C2+:
dndt(11)= 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 C2-:
dndt(12)= -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 C3:
dndt(13)= 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 C3-:
dndt(14)= 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 O2:
dndt(15)= 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 O2(a):
dndt(16)= k151*n(42)*n(15)-k157*n(16)*n(18);
%Chemin O2(b):
dndt(17)= k152*n(42)*n(15)-k158*n(17)*n(18);
%Chemin O2-:
dndt(18)= 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 O2+:
dndt(19)= 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 N2:
dndt(20)= 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 N2+:
dndt(21)= 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 CO:
dndt(22)= 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 CO+:
dndt(23)= 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 NO:
dndt(24)= 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 NO+:
dndt(25)= 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 CN:
dndt(26)= 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 CN+:
dndt(27)= 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 CN-:
dndt(28)= k14*n(3)*n(24)+k44*n(7)*n(12)+k85*n(6)*n(26)-k112*n(8)*n(28);
%Chemin O3:
dndt(29)= k163*n(1)*n(15)*n(4);
%Chemin N2O:
dndt(30)= 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 N2O+:
dndt(31)= k80*n(5)*n(30);
%Chemin NO2:
dndt(32)= 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 NO2+:
dndt(33)= k97*n(19)*n(32);
%Chemin CO2:
dndt(34)= k51*n(15)*n(37)+k55*n(15)*n(22)...
+k141*n(3)*n(15)+k146*n(6)*n(22);
%Chemin C2N:
dndt(35)= -k12*n(1)*n(35)-k28*n(4)*n(35)-k46*n(7)*n(35);
%Chemin CNO:
dndt(36)= -k7*n(2)*n(36)-k24*n(4)*n(36);
%Chemin OCN
dndt(37)= -k4*n(1)*n(37)-k51*n(15)*n(37)...
-k58*n(24)*n(37)+k56*n(15)*n(26);
%Chemin NO3:
dndt(38)= k165*n(1)*n(32)*n(4);
%Chemin C2O:
dndt(39)= k144*n(4)*n(12);
%Chemin C4:
dndt(40)= k64*n(1)*n(13)+k138*n(1)*n(14)+k148*n(9)*n(12);
%Chemin C5:
dndt(41)= k149*n(12)*n(13);
%Chemin e:
dndt(42)= 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 hv:
dndt(43)= -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);
end

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Characters and Strings finden Sie in Help Center und File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by