import numpy as np, matplotlib.pyplot as plt

pH2 = np.array([float(x) for x in '9 14.5 36 74 121 196'.split()])
H2_ads_ml = np.array([float(x) for x in '0.849 0.979 1.268 1.459 1.579 1.702'.split()])

H2_ads_mol = (H2_ads_ml*10**(-6)*10**5)/(8.315*298)
pH2 = pH2*133.3 #Pa

y_val = 1/H2_ads_mol
x_val = 1/np.sqrt(pH2)

coeff = np.polyfit(x_val,y_val,1)

plt.scatter(x_val,y_val)

reg_yval = coeff[0]*x_val + coeff[1]

print('a = ', coeff[0])
print('b = ', coeff[1])

print('C_t = ', 1/coeff[1])

plt.plot(x_val,reg_yval)
plt.savefig('Cads_mot_pH2')
plt.show()