function [P_t_inlet, P_inlet, T_i_t, P_outlet, T_o_t] = calcola_BC(NPR,M_i,P_a)

gamma = 1.4;

% Inlet;
P_t_inlet = NPR * P_a;   %[Pa]

P_inlet = P_t_inlet/( (1 + 0.5*(gamma - 1)*M_i^2)^(gamma/(gamma - 1)) );       %[Pa]
T_i_t = 300;                                                               %[K] ASSUNTO

% Outlet
P_outlet = P_a;  %[Pa]
T_o_t = 300;                                                               %[K]