Code to do Lowdin perturbation theory to any order and as efficiently as possible. Numerical and symbolic implementation.