This is a problem for MATLAB program

Please, explain step by step with code.

Isentropic flow relationships are used for studying the airflow around an airfoil during high-speed flight. The flow relationships, along with the ideal gas law, let you compare one point on a streamline to another point- see the figure below Ideal gas law (for any point): P=RT (R-1716 ft-lb/slug"R) Isentropic flow equation: 2--(T,, 3.5 T. Energy balance. Cp7, + 0.5y? = G,T2 + 0.5y? must be in ft/s) (C-6006 ft-lb/slug"RL(Note: velocities Assume that point 1 has the following conditions: pressure (P1) 2116 lb/ft2, velocity (Vi) 500 miles/hr and density 0.002377 slug/ft3. Find the temperature (T2) and velocity (V2) at point 2 for pressure (P2) values of 1000 to 2000 lb/h2in increments of 10. Plot the following: P,V5%, T2 vs. P2, T2 vs. V. Create a modular program as follows: 1. 2. Create a script that initializes the point 1 values and any other needed constants. Create and call (from the main script) the following functions: a. Ideal gas law.m this function calculates the temperature at point 1. b. sentropic.m-this function uses the isentropic relationship to calculate the temperature at point 2. c. Velocity.m- this function uses the energy balance equation to calculate the velocity at point 2. Plot vel vs press.m-this function creates the three plots. The plots should be in the format shown below d. Your final output should be a plot like the one shown. Note: The main script is used only to initialize variables and call the functions. Do not do calculations in the main script. Isentropic flow relationships are used for studying the airflow around an airfoil during high-speed flight. The flow relationships, along with the ideal gas law, let you compare one point on a streamline to another point- see the figure below Ideal gas law (for any point): P=RT (R-1716 ft-lb/slug"R) Isentropic flow equation: 2--(T,, 3.5 T. Energy balance. Cp7, + 0.5y? = G,T2 + 0.5y? must be in ft/s) (C-6006 ft-lb/slug"RL(Note: velocities Assume that point 1 has the following conditions: pressure (P1) 2116 lb/ft2, velocity (Vi) 500 miles/hr and density 0.002377 slug/ft3. Find the temperature (T2) and velocity (V2) at point 2 for pressure (P2) values of 1000 to 2000 lb/h2in increments of 10. Plot the following: P,V5%, T2 vs. P2, T2 vs. V. Create a modular program as follows: 1. 2. Create a script that initializes the point 1 values and any other needed constants. Create and call (from the main script) the following functions: a. Ideal gas law.m this function calculates the temperature at point 1. b. sentropic.m-this function uses the isentropic relationship to calculate the temperature at point 2. c. Velocity.m- this function uses the energy balance equation to calculate the velocity at point 2. Plot vel vs press.m-this function creates the three plots. The plots should be in the format shown below d. Your final output should be a plot like the one shown. Note: The main script is used only to initialize variables and call the functions. Do not do calculations in the main script