[x] Elaborate a small database with the critical temperature, critical pressure and acentric factor for pure components (species).
[x] Include automatically this data in the DB of the variable self
[x] Compute atraction and repulsion parameters for pure components
[x] Compute atraction and repulsion parameters for mixtures considering Van der Waals mixing rules and neglecting the Binary Interaction Parameter (BIP)
[x] Obtain molar volume for Van der Waals EoS
[x] Obtain molar volume for Peng-Robinson EoS
[x] Solve cubic EoS using a Newton-Raphson method
[x] Solve cubic EoS with the function roots
[x] Compute compressibility factor, Z
[x] Check results of Z for different pure components with the literature
[ ] Check results of Z for different mixtures with the literature
[x] Compute dimensionless chemical potential excess, muRT_excess
[x] Include muRT_excess in the the total dimensionless chemical potential of the system
[ ] Compute derivative of muRT_excess and include it in the matrix coefficients of the Newton-Raphson method
[ ] Include non-ideal contributions in the thermodynamic properties, i.e., enthalpy, internal energy, specific heat at constant pressure and at constant volume, entropy, adibatic index, sound velocity.