mc_tm_1d_c.m

来自「这几个算例都是MacCormack格式的,喷管计算用时间步进法,外突钝角绕流用空」· M 代码 · 共 25 行

function time=MC_TM_1D_C(N,dt,t,tU2)
% MACCORMACK`S TECHNIQUE 
% Nonconservation Form
% Time Marching

tic;
initial_data(t,N,tU2);
for h=1:t
    preCal(h);
    predictor_step(N,h,dt);
    BoundaryConditions(h+1,N);
    preCal(h+1);
    corrector_step(N,h,dt);
    BoundaryConditions(h+1,N);
    Cal_rVT(h);
end
time=toc;