solution

The compressibility factor Z via the second virial coefficient has the approximate form:

here p is the dimensionless pressure, and B* is the dimensionless second virial coefficient of the Lennard-Jones potential. Write and save the user-defined function with name Problem_4_17. The function should have the input parameters only, i.e., vector T* and integer number n (in accordance to the 2n-n Lennard-Jones potential—see preceding problem), and generates plot with the obtained compressibility factor as function of the pressures p¼0, 0.01, …, 0.07 at three isotherms T*¼0.5, 1, and 2. The values of B* should be received using the Problem_4_16 function (see problem 16) with the three-isotherm T* values of and parameter n¼6 as the input of this function.