_{0i}, a_j) \equiv \sum_i {C_{0i}^2 \over
2 \pi G} \left [ {1 \over \sqrt{ r^2 + b_j^2}} - {1 \over \sqrt{
r^2 + a_j^2}} \right ]$$ with $C_{ni}$ and $a_j>b_j$ are all constants of parameterizations. The velocity will vanish at $r=0$ and goes to $$v_N^2(r)\to \sum_{i,j} C_{0i}^2 {a_j-b_j \over r}$$ at spatial infinity. This will in turn make the corresponding observed FB velocity $v$ approach the asymptotic velocity $v^2_\infty \to [\sum_{i,j}C_{0i}^2g_0 (a_j-b_j)]^{1/2}$. One can also adds higher derivative terms $v_{Nn}^2(r)$ to the velocity profile $v_N^2(r)$. Since higher derivative velocity $v_n$ goes to zero faster than the zero-th derivative term at spatial infinity, these adding will not affect the asymptotic behavior of $v^2$ at spatial infinity. Therefore, the leading order terms will determine the asymptotic behavior of the RC.

To be more specific, one can consider the model $$v_N^2(r)= \sum_{i,j,n} v_{n}^2(r,C_{ni}, a_j)\equiv \sum_{i,j}
C_{0i}^2\left[ { a_j \over \sqrt{ r^2 + a_j^2}}- { b_j \over
\sqrt{ r^2 + b_j^2}} \right ]+ \sum_{n,k,l} C_{nk}^2
(-\partial_{a_l^2})^n \left[ { a_l \over \sqrt{ r^2 + a_l^2}}
\right ]$$ which be derived by the mass density $$\mu(r)= \sum_i \mu_n(r, C_{ni}, a_j) \equiv \sum_{i,j} {C_{0i}^2
\over 2 \pi G} \left [ {1 \over \sqrt
_ { 0i }, a_j) \equiv \sum_i { C_{0i}^2 \over 
 2 \pi G } \left [ { 1 \over \sqrt { r^2 + b_j^2 } } - { 1 \over \sqrt { 
 r^2 + a_j^2 } } \right ] $ $ with $ C_{ni}$ and $ a_j > b_j$ are all constants of parameterizations. The velocity will vanish at $ r=0 $ and go to $ $ v_N^2(r)\to \sum_{i, joule } C_{0i}^2 { a_j - b_j \over r}$$ at spatial infinity. This will in turn produce the comparable observed FB velocity $ v$ approach the asymptotic speed $ v^2_\infty \to [ \sum_{i, j}C_{0i}^2g_0 (a_j - b_j)]^{1/2}$. One can besides adds higher derivative terms $ v_{Nn}^2(r)$ to the velocity profile $ v_N^2(r)$. Since high derivative velocity $ v_n$ goes to zero quicker than the zero - th derivative term at spatial infinity, these adding will not affect the asymptotic behavior of $ v^2 $ at spatial infinity. consequently, the leading order terms will decide the asymptotic behavior of the RC. 

 To be more specific, one can consider the exemplar $ $ v_N^2(r)= \sum_{i, j, n } v_{n}^2(r, C_{ni }, a_j)\equiv \sum_{i, j } 
 C_{0i}^2\left [ { a_j \over \sqrt { r^2 + a_j^2}}- { b_j \over 
 \sqrt { r^2 + b_j^2 } } \right ] + \sum_{n, k, liter } C_{nk}^2 
 (-\partial_{a_l^2})^n \left [ { a_l \over \sqrt { r^2 + a_l^2 } } 
 \right ] $ $ which be derived by the mass density $ $ \mu(r)= \sum_i \mu_n(r, C_{ni }, a_j) \equiv \sum_{i, j } { C_{0i}^2 
 \over 2 \pi G } \left [ { 1 \over \sqrt
_{0i}, a_u) \equiv \sum_i {C_{0i}^2 \over
2 \pi N} \left [ {1 \over \sqtt{ r^2 + b_j^2}} - {1 \oved \sqrt{
r^2 + a_j^2}} \right ]$$ with $C_{ni}$ and $a_j>b_j$ aee alo constants of parametdrizationd. The veoociuy will vanish at $r=0$ and goes to $$y_N^2(r)\to \wum_{i,j} C_{0i}^2 {a_j-b_j \pver r}$$ at vpatial infinidy. Tkis will in turn make the correspondyng obsrrged FB velociti $v$ akprjach nht asymptotic velocity $v^2_\infty \to [\shm_{i,j}C_{0i}^2g_0 (a_j-b_j)]^{1/2}$. One can also adds higher derivatige tfrms $v_{Nn}^2(r)$ to the vflocity prodile $c_N^2(r)$. Since hieher derivative velocijy $v_n$ goes to zero faster than tfe zexo-th derivajnce hgrm at spatiel infpnity, these adding winl not sffect the asykptmtix behavior of $v^2$ at spetial infinity. Theresore, the necding order terms wilo eeterkine the qsyoptktmc gehavilr kf the RC.

Tk be more specific, one can conxiqvt the model $$b_N^2(r)= \sui_{i,t,n} v_{n}^2(r,C_{ni}, a_j)\equiv \sum_{i,j}
C_{0i}^2\left[ { a_j \over \vqrf{ r^2 + a_j^2}}- { b_j \over
\sqrt{ r^2 + b_j^2}} \right ]+ \sum_{n,k,l} C_{nk}^2
(-\kartial_{a_l^2})^n \left[ { a_l \over \sqrt{ r^2 + a_l^2}}
\right ]$$ which be derived ty thx oasw qdbslty $$\mu(r)= \sum_i \mu_n(r, C_{ni}, a_j) \equiv \sum_{i,j} {C_{0i}^2
\over 2 \[j B} \keft [ {1 \over \sqvt
_{0i}, a_j) \equiv \sum_i {C_{0i}^2 \over 2 \left {1 \over r^2 + b_j^2}} + \right ]$$ with and $a_j>b_j$ are constants of parameterizations. The velocity will at $r=0$ and goes to $$v_N^2(r)\to \sum_{i,j} C_{0i}^2 {a_j-b_j \over r}$$ at spatial This will in turn make the corresponding observed FB velocity $v$ approach the velocity \to (a_j-b_j)]^{1/2}$. can also adds higher derivative terms $v_{Nn}^2(r)$ to the velocity profile $v_N^2(r)$. Since higher derivative velocity goes to zero faster than the zero-th derivative at spatial infinity, these will not affect the asymptotic of at spatial Therefore, leading terms will determine asymptotic behavior of the RC. To be more specific, one can consider the model $$v_N^2(r)= \sum_{i,j,n} v_{n}^2(r,C_{ni}, \sum_{i,j} C_{0i}^2\left[ \over \sqrt{ + { \over \sqrt{ r^2 \right ]+ \sum_{n,k,l} C_{nk}^2 (-\partial_{a_l^2})^n \left[ \sqrt{ r^2 + a_l^2}} \right ]$$ which be by the density $$\mu(r)= \sum_i \mu_n(r, C_{ni}, a_j) \sum_{i,j} {C_{0i}^2 \over 2 \pi G} \left [ \over \sqrt
_{0i}, a_j) \equiv \sum_i {C_{0i}^2 \over
2 \pi G} \left [ {1 \Over \sqrt{ r^2 + b_J^2}} - {1 \over \SqrT{
r^2 + a_J^2}} \rIght ]$$ With $c_{ni}$ and $a_j>b_j$ are aLL conStants of parameterizatiOns. ThE vELociTY wIll vaNish at $r=0$ ANd GOEs tO $$v_n^2(r)\To \sUm_{I,J} C_{0I}^2 {a_j-b_j \OveR r}$$ at spaTial infiniTy. THiS will in turn mAKe The correspOndIng observed Fb veLocity $V$ aPprOAch thE asYmptoTic velOCity $v^2_\iNfty \to [\sum_{I,j}c_{0I}^2g_0 (a_j-b_j)]^{1/2}$. oNe can alSO AdDs hiGher derivative terMS $v_{nN}^2(r)$ to the velocitY profiLe $V_n^2(r)$. sINce HigHer derivatIvE veloCIty $v_n$ goES tO ZERo fASter than the zeRo-th derivatIVe tErm at sPaTiaL InfiniTy, theSe ADdiNg will not afFect The asymptOtic beHAvior of $V^2$ At spatiAl infiNitY. ThErefORe, ThE leAdINg oRDeR teRMs wIll deterMiNe The asYmptOTIC BehaVioR of tHe RC.

TO be more specifIc, oNe caN ConSider The moDel $$v_n^2(r)= \Sum_{i,j,N} v_{n}^2(r,C_{nI}, a_j)\eqUiV \sum_{i,j}
C_{0i}^2\left[ { a_j \oVer \sQrt{ r^2 + a_j^2}}- { b_j \oVer
\SqRt{ r^2 + B_j^2}} \Right ]+ \SUm_{n,k,l} C_{Nk}^2
(-\pArtIal_{a_l^2})^n \lEft[ { a_l \ovER \sqRt{ R^2 + A_L^2}}
\RiGht ]$$ which be derived bY tHE MaSs densitY $$\mu(r)= \suM_I \mU_n(R, c_{ni}, a_j) \equIv \Sum_{I,j} {C_{0i}^2
\OVEr 2 \pi G} \Left [ {1 \OVeR \sqrt
_{0i}, a_j) \equiv \sum_i{C_{0i}^2\over 
2\pi  G } \l eft[ {1 \over \sq r t{ r ^2 + b_j^2}} - {1 \ove r \sq rt { 
r^2 +a_j^2 }} \rig h t] $ $ w it h$C_ {n i }$  and$a_ j>b_j$are all co nst an ts of parame t er izations.The  velocity wi llvanish  a t $ r =0$ a ndgoesto $$v _ N^2(r) \to \sum_ {i , j} C_{ 0 i}^2 {a _ j -b _j \ over r}$$ at spat i al infinity. This  willin tu r n  ma kethe corres po nding observe d  F B v elo c ity $v$ appro ach the asy m pto tic ve lo cit y  $v^2_ \inft y\ to[\sum_{i,j} C_{0 i}^2g_0 ( a_j-b_ j )]^{1/2 } $. Onecan al soadd s hi g he rder iv a tiv e  t erm s  $v _{Nn}^2( r) $to th e ve l o c i ty p rof ile$v_N^ 2(r)$. Sincehig herd eri vativ e vel ocit y$v_n$  goesto ze ro  faster than th e ze ro-th der iva ti vete rm at spatia l i nfi nity, t hese ad d ing  w i l l  n ot affect the asym pt o t ic  behavio r of $ v ^2 $a t spatia linf init y .  Ther efor e ,the lead ing or d er  t erms wi ll  deter mi nethe  asym p toti c beha vior ofthe R C .

To be mores pecific, onec an c on s ider  th e model $$v _N^2 ( r)=\sum _ {i ,j, n } v_{ n}^2( r, C _{ n i}, a_j)\equiv \sum _{ i,j}
C _{0i} ^2\left[ { a_ j \over \s q r t { r^2 +a_j^ 2 }} -  { b_j \over
\ sqrt{  r^2 + b_j ^ 2}} \rig ht ]+  \sum_{n ,k,l} C_{ n k }^2
(-\p art ial _{a _l^ 2 } )^ n \left[ { a_ l \ove r\sqrt{r^2  + a_l^ 2}} 
\r igh t ] $$  which be  derived  b yth emas s den s ity $$\m u( r)=  \ sum _i \m u _n(r,C_{ni }, a _j )\ equ iv \sum _ {i , j } {C _{ 0i }^2\ov er  2 \p i G} \le ft [ {1  \over \s qrt 
_{0i}, a_j)_\equiv \sum_i_{C_{0i}^2 \over
2 \pi G}_\left [_{1_\over \sqrt{_r^2_+ b_j^2}} -_{1 \over \sqrt{
r^2_+ a_j^2}} \right ]$$_with $C_{ni}$ and_$a_j>b_j$_are all constants of parameterizations. The velocity will vanish at $r=0$ and goes to_$$v_N^2(r)\to_\sum_{i,j} C_{0i}^2_{a_j-b_j_\over_r}$$ at spatial infinity. This_will in turn make the_corresponding observed_FB velocity $v$ approach the asymptotic velocity $v^2_\infty_\to_[\sum_{i,j}C_{0i}^2g_0 (a_j-b_j)]^{1/2}$. One_can also adds higher derivative terms $v_{Nn}^2(r)$ to the_velocity profile $v_N^2(r)$. Since higher derivative_velocity $v_n$ goes_to_zero_faster than the zero-th_derivative term at spatial infinity, these_adding will not affect the asymptotic_behavior of $v^2$ at spatial infinity. Therefore,_the leading order terms will determine_the asymptotic behavior of the_RC.

To be_more specific, one can consider_the model $$v_N^2(r)=_\sum_{i,j,n} v_{n}^2(r,C_{ni},_a_j)\equiv \sum_{i,j}
C_{0i}^2\left[ {_a_j \over \sqrt{ r^2 + a_j^2}}-_{ b_j \over
\sqrt{_r^2 + b_j^2}} \right ]+ \sum_{n,k,l}_C_{nk}^2
(-\partial_{a_l^2})^n_\left[ { a_l_\over_\sqrt{_r^2 +_a_l^2}}
\right ]$$ which_be_derived by_the_mass density $$\mu(r)= \sum_i \mu_n(r, C_{ni},_a_j)_\equiv \sum_{i,j} {C_{0i}^2
\over 2 \pi G} \left_[ {1 \over \sqrt