def write_heat_capacity_P_polyfit(self, filename='Cp-temperature_polyfit.dat', filename_ev='entropy-volume.dat', filename_cvv='Cv-volume.dat', filename_dsdvt='dsdv-temperature.dat'):
    wve = open(filename_ev, 'w')
    wvcv = open(filename_cvv, 'w')
    for i in range(1, self._len):
        t = self._temperatures[i]
        wve.write(('# temperature %20.15f\n' % t))
        wve.write(('# %20.15f %20.15f %20.15f %20.15f %20.15f\n' % tuple(self._volume_entropy_parameters[(i - 1)])))
        wvcv.write(('# temperature %20.15f\n' % t))
        wvcv.write(('# %20.15f %20.15f %20.15f %20.15f %20.15f\n' % tuple(self._volume_cv_parameters[(i - 1)])))
        for (ve, vcv) in zip(self._volume_entropy[(i - 1)], self._volume_cv[(i - 1)]):
            wve.write(('%20.15f %20.15f\n' % tuple(ve)))
            wvcv.write(('%20.15f %20.15f\n' % tuple(vcv)))
        wve.write('\n\n')
        wvcv.write('\n\n')
    wve.close()
    wvcv.close()
    w = open(filename, 'w')
    for i in range(self._len):
        w.write(('%20.15f %20.15f\n' % (self._temperatures[i], self._cp_polyfit[i])))
    w.close()
    w = open(filename_dsdvt, 'w')
    for i in range(self._len):
        w.write(('%20.15f %20.15f\n' % (self._temperatures[i], ((self._dsdv[i] * 1e+21) / Avogadro))))
    w.close()
def write_heat_capacity_P_polyfit(self, filename='Cp-temperature_polyfit.dat', filename_ev='entropy-volume.dat', filename_cvv='Cv-volume.dat', filename_dsdvt='dsdv-temperature.dat'):
    wve = open(filename_ev, 'w')
    wvcv = open(filename_cvv, 'w')
    for i in range(1, self._len):
        t = self._temperatures[i]
        wve.write(('# temperature %20.15f\n' % t))
        wve.write(('# %20.15f %20.15f %20.15f %20.15f %20.15f\n' % tuple(self._volume_entropy_parameters[(i - 1)])))
        wvcv.write(('# temperature %20.15f\n' % t))
        wvcv.write(('# %20.15f %20.15f %20.15f %20.15f %20.15f\n' % tuple(self._volume_cv_parameters[(i - 1)])))
        for (ve, vcv) in zip(self._volume_entropy[(i - 1)], self._volume_cv[(i - 1)]):
            wve.write(('%20.15f %20.15f\n' % tuple(ve)))
            wvcv.write(('%20.15f %20.15f\n' % tuple(vcv)))
        wve.write('\n\n')
        wvcv.write('\n\n')
    wve.close()
    wvcv.close()
    w = open(filename, 'w')
    for i in range(self._len):
        w.write(('%20.15f %20.15f\n' % (self._temperatures[i], self._cp_polyfit[i])))
    w.close()
    w = open(filename_dsdvt, 'w')
    for i in range(self._len):
        w.write(('%20.15f %20.15f\n' % (self._temperatures[i], ((self._dsdv[i] * 1e+21) / Avogadro))))
    w.close()