采用lax格式用C语言模拟一维激波管问题

[1][i]=0;Р u[0][i]=u[1][i];u[1][i]=0;Р E[0][i]=E[1][i];E[1][i]=0;Р P[0][i]=P[1][i];P[1][i]=0;Р H[0][i]=H[1][i];H[1][i]=0;Р a[0][i]=a[1][i];a[1][i]=0;Р U1[0][i]=U1[1][i];U1[1][i]=0;Р U2[0][i]=U2[1][i];U2[1][i]=0;Р U3[0][i]=U3[1][i];U3[1][i]=0;Р F1[0][i]=F1[1][i];F1[1][i]=0;Р F2[0][i]=F2[1][i];F2[1][i]=0;Р F3[0][i]=F3[1][i];F3[1][i]=0;Р }Р}while(time<timeOut);Рdouble Ma[2][N];Р for(i=0;i<N;i++)Р {Р Ma[0][i]=u[0][i]/C[0][i];Р }Рprintf("X\t\t密度\t\t压力\t\t速度\t\t马赫数\n");Рfor(i=0;i<N-1;i++)Рprintf("%f\t%f\t%f\t%f\t%f\n",dx*i,D[0][i], P[0][i],u[0][i],Ma[0][i]);Рprintf("\n时间:%f\n",time);Рfp=fopen("date.txt","wb+");Рfprintf(fp,"X\t\t密度\t\t压力\t\t速度\t\t马赫数\r\n");Рfor(i=0;i<N-1;i++)Рfprintf(fp,"%f\t%f\t%f\t%f\t%f\r\n",dx*i,D[0][i], P[0][i],u[0][i],Ma[0][i]);Рfprintf(fp,"\r\n时间:%f\n",time);Рfclose(fp);Р}