 function $$f'(\gamma)-f(\gamma)\left\{{1\over\gamma}+{\gamma\over\gamma^2-1}
-{1\over kT}\right\}=  
{1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f
\,d\gamma'-\int_1^\infty (A+E)f\,d\gamma'\times\int_1^\gamma
\tilde S\,d\gamma'\right\},                                 \label{2.3}$$ while $[Df]'-Pf=0$ at $\gamma=1$ was assumed (cf. eq. (\[2.2\])). The last term in eq. (\[2.3\]) follows simply from conservation of the total number of positrons $$\int_1^\infty[A(\gamma)+E(\gamma)]f(\gamma)\,d\gamma=\int_1^\infty
S(\gamma)\,d\gamma,                                         \label{2.4}$$ which is always fulfilled if the source function has the form $S(\gamma)=\tilde S(\gamma)\times$ $\int_1^\infty
(A+E)f\,d\gamma'$ and $\int_1^\infty\tilde S\,d\gamma=1$. A regular singular point $\gamma=1$ in eq. (\[2.3\]) does not lead to any singularity of the solution, which is Maxwellian-like at the low-energy part. Equation (\[2.4\]) gives also an idea of physical meaning of term $[E(\gamma)f(\gamma)]$, that is the number of positrons with Lorentz factor $\gamma$ escaping from the plasma volume per 1 sec. The approach can be easily generalized to include inelastic processes, stochastic acceleration etc.

Equation (\[2.8\]) or (\[2.3\]) can be resolved numerically with an algorithm which reduces it to a first-order differential equation. Let $f_i(\gamma)$ is the solution obtained after the $i$-th iteration, then the equation $$f'_{i+1}(\gamma)-f_{i+1}(\gamma)\left\{{1\over\gamma}
+{\gamma\over\gamma^2-1}-{1\over kT}\right\}=
{1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f_i
\,d\gamma'-\int_1
function $ $ f'(\gamma)-f(\gamma)\left\{{1\over\gamma}+{\gamma\over\gamma^2 - 1 } 
 -{1\over kT}\right\}=  
 { 1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f 
 \,d\gamma'-\int_1^\infty (A+E)f\,d\gamma'\times\int_1^\gamma 
 \tilde S\,d\gamma'\right\ },                                  \label{2.3}$$ while $ [ Df]'-Pf=0 $ at $ \gamma=1 $ was assumed (cf.   eq.   (\[2.2\ ]) ). The last condition in eq.   (\[2.3\ ]) take after simply from conservation of the total number of positron $ $ \int_1^\infty[A(\gamma)+E(\gamma)]f(\gamma)\,d\gamma=\int_1^\infty 
 S(\gamma)\,d\gamma,                                          \label{2.4}$$ which is always fulfilled if the source affair has the form $ S(\gamma)=\tilde S(\gamma)\times$ $ \int_1^\infty 
 (A+E)f\,d\gamma'$ and $ \int_1^\infty\tilde S\,d\gamma=1$. A even singular point $ \gamma=1 $ in eq.   (\[2.3\ ]) does not lead to any singularity of the solution, which is Maxwellian - like at the broken - energy part. equality (\[2.4\ ]) gives also an idea of physical meaning of term $ [ E(\gamma)f(\gamma)]$, that is the number of positrons with Lorentz component $ \gamma$ escaping from the plasma volume per 1 sec. The approach can be well generalized to include inelastic process, stochastic acceleration etc. 

 Equation (\[2.8\ ]) or (\[2.3\ ]) can be resolved numerically with an algorithm which boil down it to a first - order differential equation. Let $ f_i(\gamma)$ is the solution obtained after the $ i$-th iteration, then the equality $ $ f'_{i+1}(\gamma)-f_{i+1}(\gamma)\left\{{1\over\gamma } 
 + { \gamma\over\gamma^2 - 1}-{1\over kT}\right\}= 
 { 1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f_i 
 \,d\gamma'-\int_1
 fujction $$f'(\gamma)-f(\gamma)\left\{{1\oyer\gamma}+{\gamma\ovet\gqmma^2-1}
-{1\ovxr kT}\rifht\}=  
{1\ovef D(\gamma)}\left\{\int_1^\gamma (A+E)f
\,d\gamla'-\unt_1^\indty (A+E)f\,d\gamma'\times\int_1^\gxmma
\tilde S\,d\gamma'\eighu\},                                 \layeo{2.3}$$ while $[Df]'-Pf=0$ aj $\gamma=1$ was dssumed (cf. eq. (\[2.2\])). Tve lcst term in eq. (\[2.3\]) follows simply from cjnservayiln of the totaj nukfer kf positrons $$\int_1^\infty[A(\gamma)+E(\gamma)]f(\famma)\,d\gemma=\int_1^\infty
S(\gakma)\,d\gamma,                                         \label{2.4}$$ whixh if always fulfklled if the source fuhction has the form $S(\gamma)=\tilde S(\gamka)\times$ $\inj_1^\nbftj
(D+E)f\,d\gamma'$ aid $\int_1^\pnfty\tilde S\,d\nsmma=1$. A regulat singular poikt $\gakma=1$ in eq. (\[2.3\]) does not lead vo any singularity os the solgtnon, which is Maxwelliqn-oike dt tve luq-endrgg 'arf. Equahioi (\[2.4\]) gives alao an idea if physical meaning os term $[E(\gamma)f(\famma)]$, ehwt is the number of positrons with Lorettz factor $\gamma$ escaping drom the plasma volumg per 1 sec. The approach can be easily generalized to includa inemxstng oeofesses, stochastic acceleration etc.

Equation (\[2.8\]) ow (\[2.3\]) vak be resolved nuierically wotj sg algorithm wfich rzsudes it to a first-ogder disferebtial equwtiom. Let $f_i(\gamma)$ is the solutiin obtained cftwr the $i$-th iteratiln, then the equstion $$f'_{i+1}(\gamma)-f_{i+1}(\gamma)\left\{{1\over\yamma}
+{\gzmma\over\gamla^2-1}-{1\over kT}\dkght\}=
{1\over D(\gamma)}\ldft\{\pnt_1^\gdmma (A+E)f_i
\,d\gamma'-\int_1
function $$f'(\gamma)-f(\gamma)\left\{{1\over\gamma}+{\gamma\over\gamma^2-1} -{1\over kT}\right\}= {1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f \tilde \label{2.3}$$ while at $\gamma=1$ was last in eq. (\[2.3\]) simply from conservation the total number of positrons $$\int_1^\infty[A(\gamma)+E(\gamma)]f(\gamma)\,d\gamma=\int_1^\infty \label{2.4}$$ which is always fulfilled if the source function has the form $S(\gamma)=\tilde $\int_1^\infty (A+E)f\,d\gamma'$ and $\int_1^\infty\tilde S\,d\gamma=1$. A regular singular point $\gamma=1$ in eq. (\[2.3\]) not to singularity the solution, which is Maxwellian-like at the low-energy part. Equation (\[2.4\]) gives also an idea of meaning of term $[E(\gamma)f(\gamma)]$, that is the number positrons with Lorentz factor escaping from the plasma volume 1 The approach be generalized include inelastic processes, acceleration etc. Equation (\[2.8\]) or (\[2.3\]) can be resolved numerically with an algorithm which reduces it to first-order differential $f_i(\gamma)$ is solution after $i$-th iteration, then $$f'_{i+1}(\gamma)-f_{i+1}(\gamma)\left\{{1\over\gamma} +{\gamma\over\gamma^2-1}-{1\over kT}\right\}= {1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f_i
 function $$f'(\gamma)-f(\gamma)\left\{{1\oVer\gamma}+{\gaMma\ovEr\gAmmA^2-1}
-{1\oVer kt}\rigHt\}=  
{1\over D(\gamma)}\leFT\{\int_1^\Gamma (A+E)f
\,d\gamma'-\int_1^\infty (a+E)f\,d\gAmMA'\timES\iNt_1^\gamMa
\tilde s\,D\gAMMa'\rIgHt\},                                 \LabEl{2.3}$$ WHiLe $[Df]'-PF=0$ at $\Gamma=1$ waS assumed (cf. Eq. (\[2.2\])). THe Last term in eq. (\[2.3\]) FOlLows simply FroM conservatioN of The totAl NumBEr of pOsiTrons $$\Int_1^\infTY[A(\gammA)+E(\gamma)]f(\gAmMA)\,d\gammA=\Int_1^\inftY
s(\GaMma)\,d\Gamma,                                         \label{2.4}$$ which is ALwAYs fulfilled if tHe sourCe FUnCTIon Has The form $S(\gaMmA)=\tildE s(\gamma)\tIMeS$ $\INT_1^\inFTy
(A+E)f\,d\gamma'$ anD $\int_1^\infty\tiLDe S\,D\gamma=1$. a rEguLAr singUlar pOiNT $\gaMma=1$ in eq. (\[2.3\]) does Not lEad to any sIngulaRIty of thE SolutioN, which Is MAxwElliAN-lIkE at ThE Low-ENeRgy PArt. equation (\[2.4\]) GiVeS also An idEA OF PhysIcaL meaNing oF term $[E(\gamma)f(\gAmmA)]$, thaT Is tHe numBer of PosiTrOns wiTh LoreNtz faCtOr $\gamma$ escaping From The plasma VolUmE peR 1 sEc. The APproacH caN be Easily gEneraliZEd tO iNCLUdE inelastic processeS, sTOChAstic accEleratIOn EtC.

equation (\[2.8\]) Or (\[2.3\]) Can Be reSOLved nUmerICaLly with aN algorIThM wHich redUcEs it to A fIrsT-orDer diFFereNtial eQuation. LEt $f_i(\gAMma)$ is the solutiON obtained afteR ThE $I$-Th ITeraTioN, then the equAtioN $$F'_{i+1}(\gaMma)-f_{I+1}(\GaMma)\LEft\{{1\ovEr\gamMa}
+{\GAmMA\over\gamma^2-1}-{1\over kT}\rigHt\}=
{1\Over D(\gAmma)}\lEft\{\int_1^\gamma (A+E)F_i
\,d\gamma'-\inT_1
 function $$f'(\gamma)-f(\ gamma)\lef t\{{1 \ov er\ ga mma} +{\g amma\over\gamm a ^2-1 }
-{1\over kT}\right\} =  
{ 1\ o verD (\ gamma )}\left \ {\ i n t_1 ^\ ga mma  ( A +E )f
\, d\g amma'-\ int_1^\inf ty(A +E)f\,d\gamm a '\ times\int_ 1^\ gamma
\tilde  S\ ,d\gam ma '\r i ght\} ,                         \labe l {2.3}$$ w hi le $ [Df]'-Pf=0$ at $\ g am m a=1$ was assum ed (cf .e q. ( \[2 .2\ ])). The l as t ter m  in eq. (\ [ 2 . 3\] )  follows simp ly from con s erv ationof  th e  total  numb er ofpositrons $ $\in t_1^\inft y[A(\g a mma)+E( \ gamma)] f(\gam ma) \,d \gam m a= \i nt_ 1^ \ inf t yS(\ g amm a)\,d\ga mm a,                                \ lab el{2 . 4}$ $ whi ch is  alw ay s ful filled  if t he  source functio n ha s the for m $ S( \ga mm a)=\t i lde S( \ga mma )\times $ $\int _ 1^\ in f t y 
( A+E)f\,d\gamma'$ a nd $ \i nt_1^\in fty\ti l de  S \ ,d\gamma =1 $.A re g u lar s ingu l ar  point $ \gamma = 1$  i n eq. ( \[ 2.3\])  d oes  no t lea d  toany si ngularit y oft he solution, w h ich is Maxwel l ia n - li k e at  th e low-energ y pa r t. E quat i on  (\ [ 2.4\] ) giv es al s o an idea of physic al  meani ng of  term $[E(\ga mma)f(\gam m a ) ]$, that  ist he number of posi trons  with Lore n tz facto r $\g amma$ es caping fr o m  the pla sma  vo lum e p e r  1  sec. The app r o achca n be ea sil y gener ali zed  to  in cl ude inela stic pro ce ss es ,sto chast i c accele ra tio netc .

Eq u ation(\[2. 8\])  o r( \[2 .3\]) c a nb e  res ol ve d nu mer ic allywith analgorit hm whichred u cesit  t o a fir st-order diff er ential equ at ion . Let$ f _i(\gamm a)$ is the solution obt a ined af ter  the$i$- th iterat ion , then  th e  equat ion $$ f'_{i +1 }(\ g a mma)- f _ {i +1} (\ gamma)\lef t \ {{1 \over \g amma }
+{\ga mma\over\gamma^2-1 } -{1 \over kT}\rig ht\ }=
{ 1 \ ov erD (\ g amm a) } \le f t \{\int_1^\gamma  (A+E)f_i\, d \g amma'-\int _ 1
 function_$$f'(\gamma)-f(\gamma)\left\{{1\over\gamma}+{\gamma\over\gamma^2-1}
-{1\over kT}\right\}=_ 
{1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f
\,d\gamma'-\int_1^\infty_(A+E)f\,d\gamma'\times\int_1^\gamma
\tilde S\,d\gamma'\right\},__ __  _  _   _  __             _\label{2.3}$$_while $[Df]'-Pf=0$_at_$\gamma=1$_was assumed (cf. eq. (\[2.2\])). The last_term in eq. (\[2.3\]) follows simply_from conservation_of the total number of positrons $$\int_1^\infty[A(\gamma)+E(\gamma)]f(\gamma)\,d\gamma=\int_1^\infty
S(\gamma)\,d\gamma, __  _        _     _  ___   _     _     \label{2.4}$$_which is always fulfilled if the source_function has the form $S(\gamma)=\tilde S(\gamma)\times$_$\int_1^\infty
(A+E)f\,d\gamma'$ and $\int_1^\infty\tilde S\,d\gamma=1$. A_regular singular_point $\gamma=1$ in eq. (\[2.3\]) does_not lead to_any singularity_of the solution,_which is Maxwellian-like at the low-energy_part. Equation (\[2.4\])_gives also an idea of physical_meaning_of term $[E(\gamma)f(\gamma)]$,_that_is_the number_of positrons with_Lorentz_factor $\gamma$_escaping_from the plasma volume per 1_sec._The approach can be easily generalized to_include inelastic processes, stochastic_acceleration_etc.

Equation (\[2.8\]) or (\[2.3\])_can be resolved numerically with_an algorithm which reduces it to_a first-order_differential equation._Let $f_i(\gamma)$ is the solution obtained after the $i$-th iteration, then_the equation $$f'_{i+1}(\gamma)-f_{i+1}(\gamma)\left\{{1\over\gamma}
+{\gamma\over\gamma^2-1}-{1\over kT}\right\}=
{1\over D(\gamma)}\left\{\int_1^\gamma (A+E)f_i
\,d\gamma'-\int_1