PlasmaFunctions: Add ion charge density function
Fix incorrect handling of xi=phi/kt < 0
Rename RadialIntegral to IntegrateOverRadius