Numerical routines for numerical calculation of Center of electrostatic potential, characteristic triangle point X(5626). Routines are mainly written in Wolfram Mathematica, but practical approximation is also implemented in C and Python.