The causes of accuracy losses in electrostatic calculations of electron- and ion-optical systems with plane symmetry are analyzed. An algorithm is offered, that takes into account these causes, reducing their contribution to the computational error. The peculiarities of the surface charge density at the electrode ends and bends are also considered.