Recently I wrote a paper "TBA Equations of 1D Hubbard Model and
High-Temperature Expansion".
Corresponding numerical program is written by Mathematica 4.1. We calculate
thermodynamic
quantities at 1/T=0.01,...0.10 at h=0, A=U/2 for U=4 and 8.
Neumerical details are here