Issue with creating a plot from two for loops

Hi,

I have an issue with creating a plot in the same figure from two for loops.

This is my code which is quite simple but I am still unable to find the issue.

a=1.63279*10^-6;

b=264.9889;

for T=373:10:1273;

Mu = (a*T.^0.5)/(1+b./T)

plot(T, Mu)

end

n = 0.666;

T0 = 273;

Mu0 = 1.375*10^-5;

for T=373:10:1273;

MuT = Mu0.*(T./T0).^n

end

[Umar Naseef]

I should be creating 2 different graphs in the same plot from the same T values and the Mu and MuT values.

Answer 1

[Star Strider]

Neither of the loops are necessary. Use element-wise operations (here, exponentiation (.^ instead of ^) and division (./ instead of /)) to create the plots:

a=1.63279*10^-6;
b=264.9889;
T=373:10:1273;
Mu = (a*T.^0.5)./(1+b./T);
figure
plot(T, Mu)
grid

n = 0.666;
T0 = 273;
Mu0 = 1.375*10^-5;
T=373:10:1273;
MuT = Mu0.*(T./T0).^n;
figure
plot(T, MuT)
grid

See the documentation on Array vs. Matrix Operations for a full explanation.

.

[Umar Naseef]

Got it. Thank you.

[Star Strider]

As always, my pleasure!

Answer 2

[David Hill]

a=1.63279*10^-6;
b=264.9889;
T=373:10:1273;
Mu = (a*T.^0.5)./(1+b./T);
n = 0.666;
T0 = 273;
Mu0 = 1.375*10^-5;
MuT = Mu0.*(T./T0).^n;
plot(T, Mu,T,MuT);

[Umar Naseef]

This also works thank you.