the electron density along the x1-axis at x2 ¼ 0 and x3 ¼ h=2, are presented in Figs. 2–4, respectively. The accuracy of the present computational method is verified by results obtained using an alternative numerical method, namely, the local radial basis function collocation method (LRBFCM) [36]. Good agreement between the two sets of numerical results can be observed. The induced electric potential has the largest value at the center of beam, x1 ¼ l=2. The biggest difference between the FEM and LRBCFM results is observed at the ends of the beam. It can be due to a reduced support domain there. Nodes in a support domain are used for approximation of physical fields in the LRBCFM. In the FEM the approximation has a local character and there is no difference between the approximation at the center or at the end of the beam. This is similarly observed, in Fig. 4, for the variation of the density of electrons.