The main objective of this research is the study of the magnetization of ferromagnetic spheres in the presence of external magnetic fields. The exact analytical solutions derived in this thesis are benchmark solutions, valuable in testing the correctness and accuracy of various approximate models and numerical methods. The total scalar magnetic potential outside the spheres, related to the magnetic field intensity, is obtained by the superposition of the potentials due to all spheres and the potential corresponding to the external field. The translational addition theorems for scalar Laplacian functions are used to solve boundary value by imposing exact boundary conditions. The scalar magnetic potential inside each sphere, related to the magnetic flux density, also satisfies the Laplace equation, which is solved by imposing the boundary conditions known from the solution of the outside field. Finally, the expressions derived are used to generate numerical results of controllable accuracy for field quantities.