Предложен метод численного моделирования самосогласованной электронно-ионной оптики в источниках электронов на основе высоковольтного тлеющего разряда. Метод основан на использовании трубок токов и позволяет учесть влияние перезарядок и разброса электронов по скоростям на фокальные энергетические параметры формируемого электронного пучка. Приведены результаты моделирования.