solution

The virial equation of state for fluids contains series of the virial coefficients from which the second virial coefficient can be calculated by the equation:

where B* is the normalized by the hard spheres volume the dimensionless second virial coefficient for Lennard-Jones-type potential function (see preceding problem); T* and r*—dimensionless temperature and intermolecular distance respectively, n—integer. Write and save the user-defined function with name Problem_4_16, and with vector T* and number n as input and B*-value as output function parameters. The function should calculate the B* for T*¼0.5, 1, 5, 10, 15, and 20 with n¼6. Use the quadgk command for calculating the semiinfinite integral of B*.