Numerical calculation of axisymmetric electrostatic modes for cold finite‐length non‐neutral plasmas