We construct numerical solutions to the higher-dimensional Einstein-Maxwell theory. The solutions are based on embedding the four dimensional Bianchi type IX space in the theory. We find the solutions as superposition of two functions, which one of them can be found numerically. We show that the solutions in any dimensions, are almost regular everywhere, except a singular point. We find that the solutions interpolate between the two exact analytical solutions to the higher dimensional Einstein-Maxwell theory, which are based on Eguchi-Hanson type I and II geometries. Moreover, we construct the exact cosmological solutions to the theory, and study the properties of the solutions.