fc.m

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

function y=fc(Ma_cal);
global gama;
 y=sqrt((gama+1)/(gama-1))*atan(sqrt((gama-1)/(gama+1)*(Ma_cal^2-1)))-atan(sqrt((Ma_cal^2-1)));