Earth’s Van Allen belts are studied by solving particular kinds of equations that could be notoriously difficult when different physical processes are acting together. In this article, we describe a numerical code that can solve these equations with unprecedented freedom from the numerous restrictions of existing models, even the ones that no other can solve. The abilities of our code could mean a breakthrough in Van Allen belt studies from the diffusive into the non-diffusive transport regime.