%% We like it HOT HOT
% Written by me bb

clear all
close all
clc

%% Adiabatic wall heat something or other

T = 300; % ambient temp in K
y = 1.4; % spec heat of air
r = .82; % recovery factor
M = 3.0; % mach number

Taw = (1+(y-1)/2*r*M^2)*T; % wall temp in K

FTaw = Taw*9/5-459.67; % wall temp in F

disp(['temp in K: ', num2str(Taw)]);
disp(['temp in F: ', num2str(FTaw)]);