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