Scribd
Upload
24 views

Untitle 1 D

Uploaded by

arminalizadeh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
24 views

Untitle 1 D

Uploaded by

arminalizadeh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

clc

clear all
i=0;
% x0 = [0.998486651119084,0.00151334888091644,4.56391007714862e-
05,8.43472695199998e-08,7.98300159331115e-
07,2.06285840269399,0.000504411630754694,0.968065987872169,-1.11341342578616e-
22,1000];
x0 =
[2.99944633942996,0.000554431835037044,0.167189882489189,0.000243253252859082,3.127
60733121853e-07,4.33205114746987,9.58701388423067e-
05,1.83387846644230,9.32362151064267e-11,1];
for T=1300
i = i+1;
fun = @(x) calculatorr12(x,T);
options = optimoptions('fsolve');
options = optimoptions(options,'MaxFunctionEvaluations', 100000);
options = optimoptions(options,'MaxIterations', 10000);
options = optimoptions(options,'OptimalityTolerance', 1e-10);
options = optimoptions(options,'FunctionTolerance', 1e-10);
options = optimoptions(options,'StepTolerance', 1e-10);
[x0,fval] = fsolve(fun, x0,options);
x0 =abs(x0);
sol(i,:) = x0;
end
phai=[0.7953:0.01:2];
phai1 = phai;
phai = [];
phai = phai1(1:9);
figure
plot(phai,sol(1:end,end),'.')
xlabel('Phai')
ylabel('T(K)')

figure

plot(phai,sol(1:end,1),'.')
xlabel('Phai')
ylabel('CO2')

figure

plot(phai,sol(1:end,2),'.')
xlabel('Phai')
ylabel('CO')

figure

plot(phai,sol(1:end,3),'.')
xlabel('Phai')
ylabel('O2')

figure

plot(phai,sol(1:end,4),'.')
xlabel('Phai')
ylabel('O')
figure

plot(phai,sol(1:end,5),'.')
xlabel('Phai')
ylabel('H')

figure

plot(phai,sol(1:end,6),'.')
xlabel('Phai')
ylabel('OH')
figure

plot(phai,sol(1:end,7),'.')
xlabel('Phai')
ylabel('H2')
figure

plot(phai,sol(1:end,8),'.')
xlabel('Phai')
ylabel('H2O')

figure
plot(phai,sol(1:end,9),'.')
xlabel('Phai')
ylabel('C(g)')

You might also like