计算流体力学作业-平板上的超声速流动的数值求解

function [rho1u1v1e1T1p1mu1k1]=BC（classImaxJmaxurTrTwprRcvcpmuoTouvTp）
switch class
    case 1     %前缘处或（1，1）点的边界条件
        u1=0;v1=0;
        T1=Tr;p1=pr;
    case 2     %左边界条件
        u1=ur*ones（1Jmax-1）;v1=zeros（1Jmax-1）;
        T1=Tr*ones（1Jmax-1）;
        p1=pr*ones（1Jmax-1）;
    case 3     %上边界条件
        u1=ur*ones（Imax-11）;v1=zeros（Imax-11）;
        T1=Tr*ones（Imax-11）;
        p1=pr*ones（Imax-11）;
    case 4     %下边界/壁面条件
        u1=zeros（Imax-11）;
        v1=zeros（Imax-11）;
        T1=Tw*ones（Imax-11）;
        p1=2*p（2:Imax2）-p（2:Imax3）;
    case 5      %右边界条件
        u1=2*u（Imax-12:Jmax-1）-u（Imax-22:Jmax-1）;
        v1=2*v（Imax-12:Jmax-1）-v（Imax-22:Jmax-1）;
        T1=2*T（Imax-12:Jmax-1）-T（Imax-22:Jmax-1）;
        p1