TIcurve function

Thermodynamic integration plot