converging_section.m

% This is used to compress the subsonic flow to sonic and find the
% conditions at the choked condition. 

function [M_out,T2,P2,x,A2] = converging_section(M_in, T_in, P_in,A_in)
    % Defininng initial properties!
    gaama = 1.4;
    M2 = M_in;
    T_0 = T_in*(1 + .5*(gaama-1)*M_in^2);
    P_0 = P_in*(1 + .5*(gaama-1)*M_in^2)^(gaama/(gaama-1));

    % Variables the need to be deleted later!!
    A1 = A_in;
    i = 0; % Need to figure out how to get the x points. So this can be deleted. 

    % Defining arrays to store future state properties. 
    A2 = [A_in];
    T2 = [];
    P2 = [];
    M_out = [];
    while M2 <= 1.0
        M_out = [M_out,M2];
        i = i+1;
        M2 = M2 + .001;
        area_ratio = (M2/M_in)*((1+.5*(gaama-1)*M_in^2)/(1+.5*(gaama-1)*M2^2))^((gaama+1)/(2*(gaama-1)));
        A2 = [A2,A1/area_ratio];
        T2 = [T2,T_0/(1+.5*(gaama-1)*M2^2)];
        P2 = [P2,P_0/((1+.5*(gaama-1)*M2^2)^(gaama/(gaama-1)))];
    end
    M_out = 1.0;
    x = linspace(0,4*(A2(1)-A2(end)),length(A2));
    x = [x,x];
    A2 = [A2,-A2];
end