; @SamarskiiB1989numerical]. The spatial grid was logarithmic, i.e. $x_{i} = \exp{(\ln x_{0} + i(\ln x_{f} - \ln x_{0}) / N)}$, $i=\{1,N\}$, while the grid in $t$ was uniform with $\tau_{k} = \tau = 0.7$. The actual number of points $N$ in the spatial grid was $\approx 2\cdot10^{3}$.

In order to evaluate the accuracy of the algorithm we performed a numerical solution of the Dirac and Schrödinger equations in the Coulomb field, i.e. the Hydrogen atom and the numerical solution of the polaron problem [@Mitra198791]. The accuracy of the calculation of eigenvalues for the Hydrogen atom was greater than 6 decimal digits. Moreover, we were able to reproduce all digits of the well known result for the ground state energy of the polaron problem $E_{\mathrm{p}}=-0.108513$ [@MiyakeA1975strong; @Mitra198791], where a similar equation of the self-consistent field arises.

At last we discuss the choice of the initial approximation $\bar{{\boldsymbol}X}_{0}$ for the unknown vector $\bar{{\boldsymbol}X}$. For this purpose, we use the variational estimation for the functional for the energy of the system, which is based on the following wave functions $$\begin{aligned}
  u_{0N}(x) = \mathrm{const}\cdot x e^{-x}, \quad v_{0N}(x) = \mathrm{const}_{1}\cdot x^{2}e^{-x},\label{eq:cont-anal-newt:18}\end{aligned}$$ for $q>0$ and $$\begin{aligned}
  u_{0N}(x) = \mathrm{const}\cdot x(1 + x^{2}) e^{-x}, \quad v_{0N}(x) = \mathrm{const}_{1}\cdot x^{2}e^{-x}, \label{eq:cont-anal-newt:19}\end{aligned}$$ for $q<0$. The constants $\mathrm{const}$ and $\mathrm{const}_{1}$ are chosen from the normalization condition Eq. (\[eq:cont-anal-newt:15\]).

Calculation of the functional (\[eq:48\]) {#sec:calc-expect-kin
; @SamarskiiB1989numerical ]. The spatial grid was logarithmic, i.e. $ x_{i } = \exp{(\ln x_{0 } + i(\ln x_{f } - \ln x_{0 }) / N)}$, $ i=\{1,N\}$, while the grid in $ t$ was uniform with $ \tau_{k } = \tau = 0.7$. The actual number of period $ N$ in the spatial power system was $ \approx 2\cdot10^{3}$. 

 In order to evaluate the accuracy of the algorithm we perform a numeral solution of the Dirac and Schrödinger equations in the Coulomb field, i.e. the Hydrogen atom and the numeral solution of the polaron trouble [ @Mitra198791 ]. The accuracy of the calculation of eigenvalues for the Hydrogen atom was greater than 6 decimal digits. Moreover, we were able to reproduce all digit of the well known result for the flat coat state energy of the polaron problem $ E_{\mathrm{p}}=-0.108513 $ [ @MiyakeA1975strong; @Mitra198791 ], where a similar equation of the self - reproducible field arises. 

 At last we hash out the choice of the initial approximation $ \bar{{\boldsymbol}X}_{0}$ for the unknown vector $ \bar{{\boldsymbol}X}$. For this purpose, we practice the variational estimation for the functional for the department of energy of the system, which is based on the following wave functions $ $ \begin{aligned } 
   u_{0N}(x) = \mathrm{const}\cdot x e^{-x }, \quad v_{0N}(x) = \mathrm{const}_{1}\cdot x^{2}e^{-x},\label{eq: cont - anal - newt:18}\end{aligned}$$ for $ q>0 $ and $ $ \begin{aligned } 
   u_{0N}(x) = \mathrm{const}\cdot x(1 + x^{2 }) e^{-x }, \quad v_{0N}(x) = \mathrm{const}_{1}\cdot x^{2}e^{-x }, \label{eq: cont - anal - newt:19}\end{aligned}$$ for $ q<0$. The constant $ \mathrm{const}$ and $ \mathrm{const}_{1}$ are chosen from the normalization condition Eq.   (\[eq: cont - anal - newt:15\ ]). 

 Calculation of the functional (\[eq:48\ ]) { # sec: calc - expect - kin
; @SalarskiiB1989numerical]. The spxtial grid was logaritimic, i.e. $x_{i} = \exp{(\un x_{0} + i(\ln x_{f} - \ln x_{0}) / N)}$, $i=\{1,N\}$, whipe the trid in $t$ was uniform dith $\tau_{k} = \tau = 0.7$. Rhe ectual number of points $K$ in fme spctmal grid was $\apkrox 2\cdot10^{3}$.

In mrder to evaludtd che accuracy of the algorithm we persormed s jumerical solujion ps ths Dirac and Schrödinger equations ih the Cmulomb field, o.e. the Hydrogen atom and tje nkmerical solution lf the polatkn [eoblem [@Mitra198791]. The accuracy of the czlculation of eigenvalues for tfe Hybrogen atom wqs hteater than 6 deciial digits. Moreover, fe were able to reprocucx alo digits of the well nnown result for tre ground scate energy of the pooaeon ptoblek $E_{\mxrhro{p}}=-0.108513$ [@JibaksA1975strojg; @Jitra198791], whers a similar equation of the sekf-sinsistent fiemd arifef.

At last we discuss the choice of the ititjal approximation $\bar{{\boodsymbol}X}_{0}$ for the unkjown vectjr $\bar{{\boldsymbol}X}$. For this purpose, we use the varhatioiau ewtlmatkin for the functional for the energy of the sysesm, wmich is based on the followonh ewve functions $$\begin{cmifned}
  u_{0N}(x) = \mathrm{clnst}\cdoj x e^{-x}, \quad v_{0N}(x) = \mayhrm{const}_{1}\cdot x^{2}e^{-x},\label{eq:conr-anal-newt:18}\end{clitned}$$ for $q>0$ and $$\begnn{aligned}
  u_{0U}(x) = \majhrm{comst}\cdot x(1 + x^{2}) e^{-x}, \quad v_{0N}(r) = \matgrm{const}_{1}\cdoh x^{2}e^{-x}, \labsu{eq:cont-anal-newt:19}\evd{akicned}$$ for $q<0$. The constants $\mwthrm{consv}$ and $\mathrm{zonsj}_{1}$ are crosen from the kmrmalization condihion Gq. (\[eq:cott-anal-newt:15\]).

Falculation of the functional (\[eq:48\]) {#sec:calc-expecj-kit
; @SamarskiiB1989numerical]. The spatial grid was logarithmic, = x_{0} + x_{f} - \ln the in $t$ was with $\tau_{k} = = 0.7$. The actual number of $N$ in the spatial grid was $\approx 2\cdot10^{3}$. In order to evaluate the of the algorithm we performed a numerical solution of the Dirac and Schrödinger in Coulomb i.e. Hydrogen atom and the numerical solution of the polaron problem [@Mitra198791]. The accuracy of the calculation eigenvalues for the Hydrogen atom was greater than decimal digits. Moreover, we able to reproduce all digits the known result the state of the polaron $E_{\mathrm{p}}=-0.108513$ [@MiyakeA1975strong; @Mitra198791], where a similar equation of the self-consistent field arises. At last we discuss the of the $\bar{{\boldsymbol}X}_{0}$ for unknown $\bar{{\boldsymbol}X}$. this purpose, we variational estimation for the functional for the system, which is based on the following functions $$\begin{aligned} = \mathrm{const}\cdot x e^{-x}, \quad v_{0N}(x) \mathrm{const}_{1}\cdot x^{2}e^{-x},\label{eq:cont-anal-newt:18}\end{aligned}$$ for $q>0$ and $$\begin{aligned} u_{0N}(x) = x(1 + x^{2}) e^{-x}, \quad v_{0N}(x) = \mathrm{const}_{1}\cdot x^{2}e^{-x}, \label{eq:cont-anal-newt:19}\end{aligned}$$ for $q<0$. The constants $\mathrm{const}$ are chosen from the condition Eq. (\[eq:cont-anal-newt:15\]). of functional {#sec:calc-expect-kin
; @SamarskiiB1989numerical]. The spaTial grid waS logaRitHmiC, i.E. $x_{i} = \eXp{(\ln X_{0} + i(\ln x_{f} - \ln x_{0}) / N)}$, $i=\{1,N\}$, whILe thE grid in $t$ was uniform with $\Tau_{k} = \tAu = 0.7$. tHe acTUaL numbEr of poiNTs $n$ IN thE sPaTiaL gRId Was $\apProX 2\cdot10^{3}$.

In Order to evaLuaTe The accuracy oF ThE algorithm We pErformed a numEriCal solUtIon OF the DIraC and SChrödiNGer equAtions in tHe cOulomb FIeld, i.e. tHE hyDrogEn atom and the numerICaL Solution of the pOlaron PrOBlEM [@mitRa198791]. THe accuracy Of The caLCulatioN Of EIGEnvALues for the HydRogen atom waS GreAter thAn 6 DecIMal digIts. MoReOVer, We were able tO repRoduce all Digits OF the welL Known reSult foR thE grOund STaTe EneRgY Of tHE pOlaROn pRoblem $E_{\mAtHrM{p}}=-0.108513$ [@MiyAkeA1975STRONg; @MiTra198791], WherE a simIlar equation oF thE selF-ConSisteNt fieLd arIsEs.

At lAst we dIscusS tHe choice of the inItiaL approximAtiOn $\Bar{{\BoLdsymBOl}X}_{0}$ for The UnkNown vecTor $\bar{{\bOLdsYmBOL}x}$. FOr this purpose, we use ThE VArIational EstimaTIoN fOR the funcTiOnaL for THE enerGy of THe System, whIch is bASeD oN the folLoWing waVe FunCtiOns $$\beGIn{alIgned}
  u_{0n}(x) = \mathrm{Const}\CDot x e^{-x}, \quad v_{0N}(x) = \mAThrm{const}_{1}\cdot X^{2}E^{-x},\LABeL{Eq:coNt-aNal-newt:18}\end{aLignED}$$ for $Q>0$ and $$\BEgIn{aLIgned}
  U_{0N}(x) = \maThRM{cONst}\cdot x(1 + x^{2}) e^{-x}, \quad v_{0N}(x) = \mAtHrm{conSt}_{1}\cdoT x^{2}e^{-x}, \label{eq:coNt-anal-newt:19}\END{Aligned}$$ fOr $q<0$. THE cONstants $\mathrm{cOnst}$ aNd $\mathrm{coNSt}_{1}$ are choSen frOm the norMalizatioN COndition eq. (\[eQ:coNt-aNal-NEWt:15\]).

calculation of THE funCtIonal (\[eq:48\]) {#Sec:Calc-expEct-Kin
; @SamarskiiB1989numerical ]. The spa tialgri d w as  log arit hmic, i.e. $x_ { i} =  \exp{(\ln x_{0} + i(\ ln x_ {f }  - \ l nx_{0} ) / N)} $ ,$ i =\{ 1, N\ }$,  w h il e the  gr id in $ t$ was uni for mwith $\tau_{ k }= \tau = 0 .7$ . The actual  nu mber o fpoi n ts $N $ i n the  spati a l grid  was $\ap pr o x 2\cd o t10^{3} $ . 
In o rder to evaluatet he accuracy of th e algo ri t hm w e p erf ormed a nu me rical solutio n  o f t heD irac and Schr ödinger equ a tio ns inth e C o ulombfield ,i .e.  the Hydrog en a tom and t he num e rical s o lutionof the  po lar on p r ob le m [ @M i tra 1 98 791 ] . T he accur ac yof th e ca l c u l atio n o f ei genva lues for theHyd roge n  at om wa s gre ater  t han 6  decim al di gi ts. Moreover, w e we re able t o r ep rod uc e all digits  of  th e wellknown r e sul tf o r  t he ground state en er g y  o f the po laronp ro bl e m $E_{\m at hrm {p}} = - 0.108 513$ [@ MiyakeA1 975str o ng ;@Mitra1 98 791],wh ere  asimil a r eq uation  of theself- c onsistent fiel d  arises.

Atl as t we disc uss  the choice  oft he i niti a lapp r oxima tion$\ b ar { {\boldsymbol}X}_{0} $for th e unk nown vector $ \bar{{\bol d s y mbol}X}$ . Fo r  t h is purpose, we  usethe variat i onal est imati on for t he functi o n al for t heene rgy  of t he  system, whic h is b as ed on t hefollowi ngwav e f unc ti ons $$\be gin{alig ne d} 
 u _{0 N}(x) = \mathr m{ con st }\c dot x e^{-x} , \qu ad v _{ 0N } (x)  = \mat h rm { c onst }_ {1 }\cd otx^ {2}e^ {-x} , \la bel{eq: cont-anal -ne w t:18 }\ en d{align ed}$$ for $q> 0$  and $$\be gi n{a ligned }   u_{0N} (x) = \mathrm{const}\cd o t x(1 +  x^ {2})e^{- x}, \quad  v_ {0N}(x ) = \mathr m{cons t}_{1 }\ cdo t x^{2} e ^ {- x},  \ label{eq:c o n t-a nal-n ew t:19 }\end{a ligned}$$ for $q<0 $ . T he constants$\m athr m { co nst } $a nd$\ m ath r m {const}_{1}$ ar e chosen f ro m  t he normali z ati on  condit ion Eq.  (\[e q :cont-a nal-newt: 15\]).

C al cula t i onof the fun ctional(\[eq:48\ ] ) {#s e c: calc- exp ect-ki n
; @SamarskiiB1989numerical]._The spatial_grid was logarithmic, i.e._$x_{i} =_\exp{(\ln_x_{0} +_i(\ln_x_{f} - \ln_x_{0}) / N)}$,_$i=\{1,N\}$, while the grid_in $t$ was_uniform_with $\tau_{k} = \tau = 0.7$. The actual number of points $N$ in the_spatial_grid was_$\approx_2\cdot10^{3}$.

In_order to evaluate the accuracy_of the algorithm we performed_a numerical_solution of the Dirac and Schrödinger equations in_the_Coulomb field, i.e._the Hydrogen atom and the numerical solution of the_polaron problem [@Mitra198791]. The accuracy of_the calculation of_eigenvalues_for_the Hydrogen atom was_greater than 6 decimal digits. Moreover,_we were able to reproduce all_digits of the well known result for_the ground state energy of the_polaron problem $E_{\mathrm{p}}=-0.108513$ [@MiyakeA1975strong; @Mitra198791],_where a_similar equation of the self-consistent_field arises.

At last_we discuss_the choice of_the initial approximation $\bar{{\boldsymbol}X}_{0}$ for the_unknown vector $\bar{{\boldsymbol}X}$._For this purpose, we use the_variational_estimation for the_functional_for_the energy_of the system,_which_is based_on_the following wave functions $$\begin{aligned}
 _u_{0N}(x)_= \mathrm{const}\cdot x e^{-x}, \quad v_{0N}(x) =_\mathrm{const}_{1}\cdot x^{2}e^{-x},\label{eq:cont-anal-newt:18}\end{aligned}$$ for $q>0$_and_$$\begin{aligned}
  u_{0N}(x) =_\mathrm{const}\cdot x(1 + x^{2}) e^{-x},_\quad v_{0N}(x) = \mathrm{const}_{1}\cdot x^{2}e^{-x}, \label{eq:cont-anal-newt:19}\end{aligned}$$_for $q<0$._The constants_$\mathrm{const}$ and $\mathrm{const}_{1}$ are chosen from the normalization condition Eq. (\[eq:cont-anal-newt:15\]).

Calculation of_the functional (\[eq:48\]) {#sec:calc-expect-kin