initial.m

来自「求解激波管问题等产生的偏微分方程组的各种数值算法。」· M 代码 · 共 16 行

% the initial value of every variable
lenX=100;left=-2;right=2;
deltX=(right-left)/lenX;
deltT=0.01;endTime=5;
R=deltT/deltX;
% initial condition
rho=ones(1,lenX+1);pres=ones(1,lenX+1);
rho(1:lenX/2)=2*rho(1:lenX/2);pres(1:lenX/2)=2*pres(1:lenX/2);             % -2<=x<0
rho(lenX/2+1:lenX)=1*rho(lenX/2+1:lenX);pres(lenX/2+1:lenX)=1*pres(lenX/2+1:lenX);   %  0<=x<=2

vel=zeros(1,lenX+1);
m=rho.*vel;

gamma=1.4;
eng=pres./(gamma-1)+0.5*rho.*vel.*vel;