Simulation of a 1D Subsonic-Supersonic nozzle flow in MATLAB using MacCormack Method

Objective: 1D Subsonic-Supersonic Nozzle governing equations in non-conservation form: The governing equations in non-conservation form is given by `(∂ρ^')/(∂t^' )=-ρ^' (∂V^')/(∂x^' )-ρ^' V^' (∂(ln⁡〖A^' 〗 ))/(∂x^' )-V'(∂ρ')/(∂x^' )` `(∂V^')/(∂t^' )=-V^' (∂V^')/(∂x^'…