clc;clear all;
%% Problem 5

%N2 Properties
R = 296.8; %N2 J/kg*K
Tv = 3521; %K
T1 = 300; %K
P1 = 101300; %Pa
T2 = 3000; %K

%% Calorically Perfect Gas
Cv_cpg = (5/2)*R;
Cp_cpg = Cv_cpg + R;
gamma_cpg = Cp_cpg/Cv_cpg;

% Since density is constant:
density_ratio = 1;

entropy_change_cpg = Cv_cpg*log(T2/T1) - R*log(density_ratio);
disp(['Assuming CPG, the change in entropy is: ', num2str(entropy_change_cpg), ' J/K', newline]);
%% Thermally Perfect Gas
Cv_tpg_1 = R*((5/2) + (exp(Tv/T1)*((Tv/T1)/(exp(Tv/T1)-1))^2));
Cv_tpg_2 = R*((5/2) + (exp(Tv/T2)*((Tv/T2)/(exp(Tv/T2)-1))^2));


entropy_change_tpg = (Cv_tpg_2-Cv_tpg_1)*log(T2/T1) - R*log(density_ratio);
disp(['Assuming tpg, the change in entropy is: ', num2str(entropy_change_tpg), ' J/K']);