hx.m

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

function y=hx(x)
global OE OL theta;
%the shape of wall
if x<=OE
    y=40;
else
    y=40+(x-OE)*tan(theta);
end