The thermodynamic properties of superconducting Nb are calculated numerically from the solution of the Eliashberg equations on the imaginary axis for several possible electron-phonon spectral densities ae(o)F(o). Comparison with experiment is made in order to see which spectrum gives the best agreement, and functional derivatives with respect to o~2 (o)F(o) are used to estimate how this agreement might be improved by small changes in a2 (o)F(o). Possible gap anisotropy effects are also considered with the help of a simple model anisotropy for the interaction.