clc P1=200 //initial pressure in kPa V1=0.1 //initial volume in metre-cube P2=500 //final pressure in kPa V2=0.2 //final volume in metre-cube W=(P1+P2)*(V2-V1)/2 //work done,obtained after derivation in book mprintf("W=%ikJ",W)