i need some help

Below is a running shock tube illustration. - at: Sonic Speed in the driver section - P1 : Pressure (absolute) of the downstream of the shock - P2: Pressure (absolute) of the upstream of the shock - P3 : Pressure (absolute) of the upstream of the expansion wave - P4 : Pressure (absolute) of the downstream of the expansion wave - T1: Temperature of the downstream of the shock - T2: Temperature of the upstream of the shock - 1 : Density of the downstream of the shock - 2 : Density of the upstream of the shock - V1 : Flow velocity of the downstream of the shock - V2 : Flow velocity of the upstream of the shock Shock Relation: P1P2=+12M2+11=1++12(M21) Expansion relation: P4P3=(121aiV1)r12 Velocity in Region 2&3,V3=V2=+12a1(My,M1) This leads to: P4P1=[1+11a4a1(MM1)]12y Contact Surface (Slip Surface) Relation: P2=P1 Overall Relation from Driver and Driven Sections: P1P4=P3P4P1P1=P3P4P1P2 Or: P1Pt=(+12M2+11)[1+11a4a1(M,Mi1)]12 So Based on given P4/P1 ratio, you need to use Matlab to calculate the corresponding M6 then calculate a few more parameters in the same program using the following relations: 12=1++11(M21)Ms2T1T2=(a1a2)2=2+1M2(+12M2+11)(1+21M2)a1V2V1=+12(MM1) Use your program to solve initial driver section pressures at: P4=40psi, 60psi, and 80 psi respectively. Use input(prompt) function to allow users to enter the P4 pressure manually for calculations, so you don't need to change your program for different P4 values. You will need to submit the matlab code .m file and a brief report in word file explaining your algorithm used in the code and add results of your calculations for above pressures P4