mc_tm_1d_nc.m

来自「Matlab下实现的CFD算法」· M 代码 · 共 24 行

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

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