shock capturing on flat plate
 

i am solving problem described in chapter 10 of CFD by john D.anderson jr. in TAUXY subroutine will T(x,y)be calculated at all nodes 1-1 to 70-70.please also guide me how case 1,2,3,4 will be executed. thanking in advance Muhammad zubair asif