• The spatial fractional quantum mechanics is applied for the electric screening effects. • The reduced form of the fractional Schrödinger equation for the electric screening potential is derived. • The equation is solved by applying the numerical simulation methods. • The wave function and the probability in the stationary state are found. In this article, we simulate the spatial form of the fractional Schrödinger equation for the electrical screening potential using Riemann-Liouville definition of the fractional derivatives and the numerical simulation methods. We find the wave function of systems described by the electrical screening interaction potential with a specific electrical permittivity. We find and apply the dimensionless formalism of the spatial fractional Schrödinger equation in case of the electrical screening interaction potential in the stationary state. We find the probabilities and the amplitude of the wave functions for multiple values of the spatial fractional parameter of the fractional Schrödinger equation. In every case of the spatial fractional parameter, we take multiple values of the dimensionless energy. The algorithm of this work is applied in case of the systems which obey the electrical screening interaction potential to be described by the fractional quantum mechanics such as the plasma systems, tokamak and some colloidal dispersion, all we need, the parameters of the system and applying the method.