clc
P1=100*10^3
v1=1
v2=1/16
n=1.25
T1=300
P2=P1*((v1/v2)^n)
mprintf("P2=%fMPa\n",P2/(10^6))
T2=(T1*P2*v2)/(P1*v1)
mprintf("T2=%fK\n",T2)
R=8.314
W=(R*(T1-T2))/(n-1)
mprintf("W=%fkJ/mol\n",W/1000)
gama=1.4
q=((R*(T2-T1))/(gama-1))+W
mprintf("q=%fkJ/mol",q/1000)