where F is a given vector, X is the vector of unknown degrees of freedom, and M is the given solution matrix. In all three cases, the entries of X allows to represent the (approximate) reaction potential W and can be used to compute the (approximate) electrostatic solvation energy that is of the form