MATLAB code is needed! Please help!!

Using the Peng-Robinson EOS, do the following: (a) Construct a P-V diagram for methane. Use Isotherms in 10-degree increments starting at 75 degrees Kelvin until you are approximately 40 degrees above the critical point. Calculate the vapor pressure of methane as a function of temperature. Compare your computed results with the literature and describe the similarities and differences. On the P-V diagram, use the limits of stability derived in class to construct the metastable regime in your computer code. Use the requirement that fugacity's are equal in the liquid and vapor states to construct the coexistence lines. Compare the envelope of coexistence curves with the result in part (d). ln f^v (T, P)/P = 1/RT integral^V = Z^v RT/P _V = infinity (RT/V - P) DV - ln Z^v + (Z^V - 1) = (Z^V - 1) - ln (Z^v - B) - A/2 Squareroot B ln [Z^v + (1 + Squareroot 2) B/Z^v + (1 - Squareroot 2) B] ln f^L (T, P)/P = 1/RT integral^V = Z^L RT/P _V = infinity (RT/V - P) DV - ln Z^L + (Z^L - 1) = (Z^L - 1) - ln (Z^L - B) - A/2 Squareroot 2B ln [Z^L + (1 + Squareroot 2) B/Z^L + (1 - Squareroot 2) B]