这几个算例都是MacCormack格式的,喷管计算用时间步进法,外突钝角绕流用空间步进法.都写成模块形式,便于查错.

				主程序：NFlow.m
				计算采用MacCormack格式，时间步进法，给定初始条件，入口和出口边界条件，计算使用了无量纲处理。
				运行后可以生成几个全局变量
				        rou 密度
				        T   温度
				        V   速度
				        A   喷管型面
				        time 运行时间
				举例
				>> NFlow
				
				time =
				
				   30.3530
				
				>> plot(rou(:,31))
				>> plot(rou(1600,:))
				
				
				子模块说明：
				        function time=MC_TM_1D_NC(N,dt,t)
				        %MacCormack非守恒型无激波跨声速喷管主程序
				        %N是划分网格点数
				        %dt是时间步长
				        %t是计算的时间格点数 
				        %time返回计算所用时间
				
				        function time=MC_TM_1D_NC(N,dt,t,PN)
				        %MacCormack非守恒型无激波亚声速喷管主程序
				        %N是划分网格点数
				        %dt是时间步长
				        %t是计算的时间格点数 
				        %PN出口压力
				        %time返回计算所用时间
				
				        function time=MC_TM_1D_C(N,dt,t,tU2)
				        %MacCormack守恒型无激波跨声速喷管主程序
				        %N是划分网格点数
				        %dt是时间步长
				        %t是计算的时间格点数 
				        %tU2计算初始流量
				        %time返回计算所用时间
				
				        function time=MC_TM_1D_C(N,dt,t,pN,Cx)
				        %MacCormack守恒型有激波跨声速喷管主程序
				        %N是划分网格点数
				        %dt是时间步长
				        %t是计算的时间格点数 
				        %pN是出口压力
				        %Cx人工粘性系数
				        %time返回计算所用时间
				
				%子模块的变量对各个计算程序略有不同
				
				 
				        function initial_data(t,N)
				        %初始程序
				           function y=NF_Shape(N)
				           %面积函数
				           function init_condition(N)
				           %初始条件
				        
				        function preCal(h,dt)
				        %处理计算用到的计算变量
				
				        function predictor_step(N,h,dt);
				        %预估
				
				        function BoundaryConditions(h,N)
				        %边界条件
				
				        function corrector_step(N,h,dt);
				        %校正
				
				        function Cal_rVT(h);
				        %计算 rou V T