[Pw_forum] help: free energy

[hailin yu]

Dear Stefano B.,
       Thank you for you reply. I have run the code (which you give me  in 
the attach) onto my date and the results are well good agreement with your 
paper (PRB,57,10421(1998)).
       But I am still confused about the code. In the 
Equ(1)(PRB,57,10421(1998)), there is "h" in the
ln[2sinh(h*omega/kB*T)] part. But I did't find "h" appear in the code. So 
can you tell me why "h"
did't appear in the code.
     Thank you very much!
regards!
     yuhailin

_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE! 
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/