gamma = 1.4;
R = 287;
p1 = [6.78,2.79,9.8,11.93,7.81,8.83,5.81,7.84,5.29,4.83,9.77];
T1 = [23,23,22.5,23,22,22,22,22.5,24.5,22.5,24];
for K = 1:length(p1)
p4 = 29.85;
fun = @(p2) p2/p4*(1 -((gamma-1)*(p2/p1(K)-1))./(2*gamma*(2*gamma+(gamma+1)*(p2/p1(K)-1))).^.5).^(-2.*gamma./(gamma-1))-1;
result(K) = fzero(fun,10)/p1(K) ;
end
for i = 1:length(T1)
fun1 = @(T2) (result(K)*T1(i)/T2)*(((gamma+1)/(gamma-1))+result(K))/(1+(((gamma+1)/(gamma-1))*result(K)))-1;
Temp(i) = fzero(fun1,25);
end
for Q = 1:length(T1)
a(Q) = sqrt(gamma*R*T1(Q));
end
M_s = w_s(j)./a(Q); =>