function [VIon]=IonEquilibriumPotential(positiveInR, positiveOutR)
%Calculate Equilibrium Potential
%Reverse for Green
if positiveInR==0
    VIon=58*log(positiveOutR);
elseif positiveOutR==0
    VIon=0;
else
    VIon=58*log(positiveOutR/positiveInR);
end