ger.m

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

function y=ger(tx,x0,x1)
 x2=x1-(fc(x1)-tx)*(x1-x0)/(fc(x1)-fc(x0));
 n=1;
 while (abs(x1-x0)>=1.0e-4)&(n<=100000000)
     x0=x1;
     x1=x2;
     x2=x1-(fc(x1)-tx)*(x1-x0)/(fc(x1)-fc(x0));
     n=n+1;
 end
 y=x2;