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