似乎已经完成了phase lag的校正, Oh Yeah,
Reversivle Heat capacity = best linear fitting slope 跟另外一个近似方法实质完全一样; 应该用Reversible Heat capacity = best linear fitting slope * cos( phase lag) 这样才行. 我尝试了一下, 果然如此. 结果得到了大大改善. 其中
