import numpy as np
import matplotlib.pyplot as plt

zp = 1
zm = -1

R, T = 8.314, 298

c_list = np.linspace(0,0.1,150)

mu_list = c_list

gamma = np.exp(-np.sqrt(mu_list/0.725))

plt.plot(c_list, R * T * np.log(gamma))
plt.ylabel(r'$\Delta_{mix}H$')
plt.show()

plt.plot(c_list, gamma)
plt.ylabel(r'$\gamma_{\pm}$')
plt.show()