Week 7 NOzzle

Objective: close allclear allclc%initial conditionsn=31x=linspace(0,3,n);dx=x(2)-x(1);gamma=1.4;rho=1-0.3146*x;t=1-0.2314*x;v=(0.01+1.09*x).*t.^0.5;a=1+2.2*(x-1.5).^2;%time stepsnt=5;dt=0.001;%outer time loop for kt=1:nt%continuity equationfor j=2:n-1 %predictor method dvdx=(v(j+1)-v(j))/dx; drhodx=(rho(j+1)-rho(j))/dx; dtdx=(t(j+1)-t(j))/dx;…