This book presents a comprehensive mathematical approach for solving shastic magnetic field problems. It discusses variability in material properties and geometry, with an emphasis on the preservation of structural physical and mathematical properties. It especially addresses uncertainties in the computer simulation of magnetic fields originating from the manufacturing process. Uncertainties are quantified by approximating a shastic reformulation of the governing partial differential equation, demonstrating how statistics of physical quantities of interest, such as Fourier harmonics in accelerator magnets, can be used to achieve robust designs. The book covers a number of key methods and results such as: a shastic model of the geometry and material properties of magnetic devices based on measurement data; a detailed description of numerical algorithms based on sensitivities or on a higher-order collocation; an analysis of convergence and efficiency; and the application of the developed model and algorithms to uncertainty quantification in the complex magnet systems used in particle accelerators.