Vectorize, Vectorize, Vectorize!
I recommend that you get in the habit of always using .* and ./ and .^ for everything (except perhaps when using scalar numeric constants as the first or last thing at the same artithmetic priority). Even if you are writing with variables you know are scalar now, if later you extend the variable to non-scalar, if you are in the habit of using the dot versions of the operators then you will not need to debug the failures later.
The exception, of course, being for cases where you are deliberately doing inner product (algebraic matrix multiplication) or deliberately doing matrix left divide \ or deliberately doing matrix right divide /
Even when using scalar numerics, keep in mind that 1/x with scalar x works the same as 1./x but that if you extend x to be non-scalar then 1/x suddenly will do something completely different. It is always safe to use / to divide by a scalar numeric value (unless the numeric value is being raised to a power, as power has higher priority)
It is safe to have expression * scalar_numeric . It is safe to have scalar_numeric * variable. However, expression * scalar_numeric * variable is not always safe, since expression * scalar_numeric might be non-scalar and non-scalar * variable could be a problem...
%Variables
g=5/3;
g1=(g-1)/2;
g2=(1-g)/2;
Ac=1.267*10^-4;
Pc=400:600; % Pc range from 400 to 600
Pcc= 400; %Variable check
T0=350;
P0=615;
Rg=207.85;
Tcc = T0*(Pcc/P0)^((g-1)/g)-1 %Check at Pcc=400
mc= sqrt((1/g1)*((P0/Pcc)^((g-1)/g)-1))%Check at Pcc=400
mdotc = (Pcc/(sqrt(Tcc))*mc*sqrt(g/Rg)*Ac)%Check at Pcc=400
Tc = T0.*(Pc./P0).^((g-1)./g)-1;
m = sqrt((1./g1).*((P0./Pc).^((g-1)./g)-1));
mdot = (Pc./(sqrt(Tc)).*m.*sqrt(g./Rg).*Ac);
plot(Pc,mdot)
xlabel('Background Pressure (Pc)')
ylabel('Mass flow rate (kg/s)')
