Solution to Exercise 8.2

clc
disp('A gas is expanded in an engine cylinder, following the law PV^1.3=c')
disp('The initial pressure is 2550 kPa and the final pressure is 210 kPa.')
disp('If the volume at the end of expansion is 0.75 m3,')
disp('Compute the work done by the gas.')
disp(' ')                    % Display blank line
n=1.3;
P_i=2550;                    % Initial pressure
P_f=210;                     % Final pressure
V_f=.75;                     % Final volume
V_i=(P_f*(V_f^n)/P_i)^(1/n); % Initial volume
c=P_f*V_f^n;
v=V_i:.001:V_f;              % Creating a row vector for volume, v
p=c./(v.^n);                 % Computing pressure for volume
WorkDone=trapz(v,p)          % Integrating p*dv