1) Added the case of electrode atoms 2) Generalized to two types of cations for the bulk liquid case
Now possible to calculate the cage function between one ion and the whole electrode--not very useful, though!