I am very unsure of why this does not work. I want to plug in a range of Xa into V equation and plot Xa vs V
k=1; %mol3L-3sec-1
Fa0=1; %molA/sec
v0=1; %L/sec
Ca0=Fa0/v0;
Xa=(0:.05:1);
V=(Fa0-(Ca0*v0*((1-Xa)/(1-(.5*Xa)))))./(-k*((Ca0*((1-Xa)/(1-(.5*Xa)))).^-2));
plot(Xa,V);
xlabel('Conversion');
ylabel('Volume');

 Réponse acceptée

Walter Roberson
Walter Roberson le 19 Fév 2021
k=1; %mol3L-3sec-1
Fa0=1; %molA/sec
v0=1; %L/sec
Ca0=Fa0./v0;
Xa=(0:.05:1);
V=(Fa0-(Ca0.*v0.*((1-Xa)./(1-(.5*Xa)))))./(-k.*((Ca0.*((1-Xa)./(1-(.5*Xa)))).^-2));
plot(Xa,V);
xlabel('Conversion');
ylabel('Volume');
In short: get in the habit of using ./ instead of / (except perhaps for simple fractions like 2/3 ) . The / function is not element-by-element division: it is "matrix right divide", in which A/B acts similar to A*pinv(B) where * is algebraic matrix multiplication.

1 commentaire

Kara Scully
Kara Scully le 19 Fév 2021
thank you so much! i honestly was not sure when to use ./ specifically if it had to be for like every step in the equation but now i know :)

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur 2-D and 3-D Plots dans Centre d'aide et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by