equation code not getting right answer

Hello, I'm trying to code this expression, this is my attempt but I'm not getting the same answer/value, not sure what I'm doing wrong. Any help/guidance will be appreciated.
% Given values
%AInfinity = 1;
%AInlet = 6.5;
MInlet = 0.5;
MInfinity = 2;
gamma = 1.4;
a = (1 + (gamma-1/2)* MInfinity^2)
b = (1 + (gamma-1/2)* MInlet^2)
b1 = (MInlet/MInfinity)
r = b1*(a/b)^(gamma + 1/(2*gamma -1))

 采纳的回答

% Given values
%AInfinity = 1;
%AInlet = 6.5;
MInlet = 0.5;
MInfinity = 2;
gamma = 1.4;
a = (1 + ((gamma-1)/2)* MInfinity^2) % for numerator
a = 1.8000
b = (1 + ((gamma-1)/2)* MInlet^2) % for numerator
b = 1.0500
b1 = (MInlet/MInfinity)
b1 = 0.2500
r = b1*(a/b)^((gamma + 1)/(2*(gamma -1))) % for num and den
r = 1.2595

3 个评论

use the parenthesis for num & den in the expression
Thank you, I just would like to know if these are well coded? Thanks
p02 = pa * (1 + eta_d * (gamma - 1) / 2 * Ma^2) ^ (gamma / (gamma - 1))
p2 = p02 * (1 + (gamma - 1) / 2 * M2^2) ^ (gamma / (gamma - 1))
I just would like to know if these are well coded?
Yes, the two expressions are correctly coded.

请先登录,再进行评论。

更多回答(0 个)

类别

帮助中心File Exchange 中查找有关 Gamma Functions 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by