主程序：PMEW.m
計算采用MacCormack格式，空間步進法，給定來流初始條件，壁面邊界條件采用Abbett`s boundary，采用貼體網格，加入人工粘性，計算未使用無量綱處理。（輸出為計算空間）
運行后可以生成幾個全局變量
   rou 密度
   p   壓力
   T   溫度
   u   x方向速度
   v   y方向速度
   iksi 計算的各個步長
舉例
>> PMEW

time =

    3.2550

>> mesh(u)


子模塊說明：
   function time=MC_SM_2D_C(M,N,dyita,C,Cy)
   %MacCormack守恒型外突鈍角繞流主程序
   %M,N是劃分網格數
   %dyita是yita所選取步長
   %C是計算ksi步長的Courant Number
   %Cy是人工粘性系數
  
   function initial_data(M,N)
   %初始程序
     function init_condition(irou,iT,iMa)
     %處理初始條件

   function cal_dksi(M,i,dyita,C); 
   %計算ksi步長

   function preCal(M,i,Cy)
   %處理計算用到的計算變量

   function predictor_step(M,i,dyita);
   %預估

   function Cal_rVT(i);
   %計算 rou V T

   function BoundaryConditions(M,i)
   %邊界條件

   function corrector_step(M,i,dyita);
   %校正