Summary: This paper illustrates some of the advantages of using the boundary fitted coordinate technique in the numerical simulation of electron-optical systems. The boundary fitted coordinate technique can be used as an automatic mesh generator for the finiteelement solution of field distributions, while at the same time allowing for the possibility of high order interpolation on the mesh, typically required for the accurate plotting of electron trajectory paths. For the finite difference method of solving field distributions, the boundary fitted coordinate method can generate regular meshes, which fit curved shaped boundaries.