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