The Landau-Zener-Stückelberg-Majorana (LZSM) problem models diabatic transitions between energy levels in quantum two-level systems with an avoided level-crossing. The diabatic transition is a consequence of quantum tunneling in energy space when the system's Hamiltonian is perturbed with a fast-acting bias. The probability of transition between the energy states for a linear bias is known as the LZSM transition probability. The objective of this work is to investigate the LZSM problem through analytical and numerical lenses. The LZSM transition probability is derived in two ways. The first approach is based on Majorana's solution using contour integrals. The second derivation follows Landau's quasi-classical treatment. The derivations demonstrate methods for transitions in the presence of time-dependent perturbations. The ubiquity of the two-level system is discussed and an application on qubits concerning LSZM interferometry is presented, with the latter arising after considering periodic biases. Lastly, a simulation of the two-level system is conducted using Trotter-decomposed time-evolution operators, perturbation theory, and vectorization. The simulated transition probabilities for linear and periodic biases are obtained for varied parameters. The results show that the simulation achieves an accurate and efficient emulation of the LZSM problem.