SymBoltz represents cosmological models with symbolic equations grouped in physical components for the metric, gravity, and particle species. This is compiled to a numerical problem that splits equations into background and perturbation stages and generates fast code for ODE functions f and Jacobians J. The problem is then solved, and the result is stored in a solution object that gives access to any model variable S.