Mathematica_scripts/coulomb_speedfunctions.nb