  We describe the first convergent numerical method to determine static black
hole solutions (with S^3 horizon) in 5d compactified spacetime. We obtain a
family of solutions parametrized by the ratio of the black hole size and the
size of the compact extra dimension. The solutions satisfy the demanding
integrated first law. For small black holes our solutions approach the 5d
Schwarzschild solution and agree very well with new theoretical predictions for
the small corrections to thermodynamics and geometry. The existence of such
black holes is thus established. We report on thermodynamical (temperature,
entropy, mass and tension along the compact dimension) and geometrical
measurements. Most interestingly, for large masses (close to the
Gregory-Laflamme critical mass) the scheme destabilizes. We interpret this as
evidence for an approach to a physical tachyonic instability. Using
extrapolation we speculate that the system undergoes a first order phase
transition.
