In the calculation of the compressibility factor for methane gas, the virial equation of state interrupted in the fourth term will be used.

=1+/+/2+/3 ;=-53.4cm3/mol;=2620cm6/mol2 ;=5000cm9/mol3

a) Plot the values for methane gas at 273 K as a function of pressure (range 0-3000 kPa).

Calculations must be performed via a computer program (eg Program1.py) written in Python. Perform the calculations based on pressure values of 100 kPa in the range of 0-3000 kPa. It is necessary to export the results to a file (such as result.txt) and open the file in Excel to produce a graph of with appropriate formatting.

b) With an acceptable error of 3%, determine the pressure at which the ideal gas equation can be used, based on the results in part (a).

c) Write a computer program (with Python, for example Program2.py) that can determine the state of b). It is expected that the program will ultimately yield a single result as follows. During this analysis, increase the pressure in steps of 1 kPa, starting from 0.

For methane, the ideal gas equation can be used up to the ABC kPa pressure value based on 3% margin of error.

d) The virial equation of state in the form of = 1 + /, truncated from the second term, can be used up to what pressure value (acceptable error limit of 3%) is qualified? It is expected that the calculation for this situation has also been performed in Program1.py.

The content of the result.txt file that program1.py will create can be as follows: P,kPa Z (4 term equation) Z (2 term equation)

------ ------------------------ --------------------- -----

0 1.0000 1.0000

100 a.bbbb a.bbbb ....

3000 y.zzzz y.zzzz.