Present Value in the Case of General Interest Rate and General Rate of Inflation Početní základy finanční matematiky sylabus 5. přednášky (5.9) 24. 10. 2005 4 Úroková míra veličin, které nejsou v ustáleném stavu 4.1 Abstrakt Ve finanční matematice se používají dva pojmy: jednoduché a složené úročení. Úroková míra za jednotku času je v obou případech rovna f(x + l)-f(x) kde / je stavová funkce. Protože / nemusí být v bodě x + 1 vůbec definována (předpokládáme, že je definována na nějakém pravém okolí (x,x + h) bodu a;), počítá se v prvním případě úroková míra za jednotku času takto £=(1+f(x + h)-f(x) a ve druhém případě takto: = f f(x + h) - f(x)\ 1 v v m ) h a říká se, že je / úročena složené, pokud £ nezávisí na h a že je úročena jednoduše, pokud r\ nezávisí na h. / je úročena složené, právě tehdy, když je v ustáleném stavu. 4.1.0.1 Pokud je / v ustáleném stavu (tj. roste exponenciálně), nezávisí £ ani na x. 4.1.0.2 Pokud r) nezávisí na x, je / nulová konstanta. 4.1.0.3 Proto je číslo £ mnohem vhodnější míra než číslo t)av případech, které se dějí ze své a ne z naší vůle, jako je třeba inflace používáme číslo £. 4.1.0.4 Máme dvě různé míry růstu veličiny. Jejich nevýhodou je, že obě závisí nejen na bodě x, ve kterém růst měříme, ale i na délce času h. 4.1.0.5 V makroekonomii se jako míra růstu veličiny / používá hodnota výrazu f'(x) f(x + h)- f(x) 1 9 = —— = lim -—---- = hm r) f(x) h->o f(x) h h->0 4.1.0.6 Veličiny /, pro které je 6 konstantní jsou v ustáleném stavu, za povšimnutí stojí fakt, že se nepočítá limita £, ale r\. 4.1.0.7 Ve finanční matematice také se používá pojem spojité úročení namísto složené úročení a v tom případě je úroková míra jméno pro číslo 6, ne pro číslo £. Tato neobratnost ve vyjařování (nejde o jiné úročení, stav kapitálu pokaždé roste exponenciálně, nebo tvoří geometrickou posloupnopst, je-li čas diskrétní, ale o jiné číslo, které je nazváno úrokovou mírou) je ve většině učebnic finanční matematiky znejasňována ještě tím, že se předvádějí různé početní hříčky, kdy se z čísla £ dostane 6. (známá je například tak, kdy se za h postupně dosazují prvky posloupnosti (l/2n)^L0, číslům (1/2)™ se říká frekvence připisování úroků) 4.1.0.8 Cílem této poznámky je: 4.1.0.9 Najít limitní tvar (pro h —> 0) míry £, který lze použít namísto 6 4.1.0.10 Najít vzorec pro budoucí hodnotu kapitálu, v němž bude mít pojem úroková míra tentýž význam, jaký má při složeném úročení, totiž vzorec, do kterého bude možno dosadit dostatečně obecnou funkci, která bude vyjadřovat závislost úrokové míry na čase a který přejde do známého tvaru xt = xo(l + «)*, pokud bude úroková míra i ľit konstantní (poznamenejme, že vzorec xt = xoe-l používaný v učebnicích finanční matematiky tuto vlastnost nemá). 4.1.0.11 Na případech, jejichž chování je autonomní (nezávislé na naší vůli), jako je inflace, ukázat, že volba tohoto vzorce není libovolná, ale nutná: Totiž vzorec xt = xoe-l ln(1+& je limitním případem vzorce xt = xo(l + «)* tomto smyslu: pokud funkcionál, který přiřazuje úrokové míře budoucí hodnotu kapitálu má být spojitý, je budoucí hodnota při úročení libovolnou úrokovou mírou vynucena známými budoucími hodnotami při úroočení (po částech) konstantní úrokovou mírou, protože množina po částech konstantních funkcí je hustá v množině všech přípustných funkcí. 41 Present Value in the Case of General Interest Rate and General Rate of Inflation 4.2 Odvození hlavní formule 1. Introduction: 4.2.0.1 Nechťjsou aregované ceny nějaké funkce závislá na čase F = CPI: M—► M (1) Míra inflace za dobu (t±, í2) je relativní hodnota jejich přírůstku: L{tlM = LÍ{tlM) = EM_£M. (2) Takto definováno je l zobrazení l: M2 —> M (3) 4.2.0.2 1.1. Definition: (see. [2]) Steady state is a situation in which the various quantities splňují následující rovnici: D (z) (í) = k-x(t) (4) každá pro nějakou konstantu k. Tedy, každá funkce v ustáleném stavu závisí na dvou parametrech A a B a je tvaru x^A- eBx. (5) 4.2.0.3 Jediná quantity, kterou se nyní zabýváme j sou agregované ceny. 4.2.0.4 1.2. Lemma: Nechťjsou v čase (Ti, T2) ceny v ustáleném stavu. Pak Vk:Víi:Ví2:(íi,íi+k> c (Ti,T2) A (í2,í2 + k) c (Ti,r2) => ,h + k)) = i((í2,í2 + «» (6) i. e. i((í, í + k)) depends only on the length k of the interval ((t,t + n) but not on its beginning t. In a steady state we can define "i ' a; ^ ^(Tl^+o;)) (7) Při tom platí: r(r) = (l + i(í))í-l (8) takže známe-li hodnotu í v nějakém í € 1 můžeme rozšířit definiční obor i na celé ffi: ^ - ^ í(r) = (l + i(í))í-l (9) r Takto definované je íľ míra inflace agregovaných cen, které jsou v ustáleném stavu na ffi a které vyhovující rovnici F(r)=F(0)-(l + t(t))r (10) pro každé t G (To, To + ť) Tedy zejména, pokud zvolíme hodnotu t za jednotku času, máme: a naopak id) = + -l=|'Jn"+^-r"') . (12) Takto můžeme definovat míru inflace za jednotku času pokud jsou agregované ceny v ustáleném stavu na nějakém, libovolně malém, pravém okolí bodu iGl 4.2.0.5 V obecném případě, i pokud není stav ustálený, můžeme normovat míru inflace naměřenou na intervalu (t,t + n) s použitím vztahu (12). Definujeme: 42 Present Value in the Case of General Interest Rate and General Rate of Inflation 4.2.0.6 1.3. Definition: Míra inflace agregovaných cen F na intervalu (t,t + n) za jednotku času je: (F(t + K)-F(t)\1/K 4.2.0.7 1.4. Note: ľ závisí jen na hodnotách F v krajních bodech intervalu (t, t + n) a je stejná, jako by byla hodnota inflace v ustáleném stavu, ve kterém by agregované ceny v čase t byly F(t) a v čase t + k by byly F(t + n) tj. podle (5) by byly splněny rovnice: F(t) = AeBt (14) F{t + K)=AeB^ K ' které jsou ekvivaklentní s rovnicemi V F(t) ) F(t) J u 4.2.0.8 Nás ovšem zajímá míra inflace ľ = limK_>0+ ľ( —, n), v jednom časovém okamžiku t. Pokud v tomto okamžiku mají agregované ceny nespojitost prvního druhuje míra inflace l nenulová a je rovna ,m _ limT^t+ F(t) - limT^t- F(t) [) ~ lim_t- F(r) [ib) a míra inflace ľ za jednotku času bude rovna oo. 4.2.0.9 Pokud jsou v tomto okamžiku agregované ceny spojitě závislé na čase je míra inflace v tomto okamžiku 0. Pokud jsou navíc v tomto okamžiku agregované ceny diferencovatelné podle času, existuje limita (15) k^Oa platí: D(-F)(t)t lim A = e "(') F(t) k—*0 iimS = ^M F(t) (17) a míra inflace za jednotku času je lit) = l{l)AeB{t+1'-AeBt =e*-l = e*mi-l = cPOnoW) (18) a máme ln(i(í) + l) = D(lnoF)(í) (19) This is the main result of this part. Poslední rovnici můžeme psát také ve tvaru: B{F){t) =ln(i(í) + 1) -F{t) (20) 4.2.0.10 One can compare it with a equation B(F)(t) = ((t) ■ F(t) (2d která se obvykle používá k definici míry růstu of different quantities v makroekonomii a jejíž speciálním případem je (5). Podle toho, co jsme odvodili, by měla být makroekonomická míra růstu £ rovna ln(l + l). 4.2.0.11 Podotkněme ještě, že taylorův polynom stupmně 1 v bodě 0 funkce x i—> ln(l + a;) je stejný jako taylorův polynom identity, že tedy obě funkce mají dotyk řádu 1 (first order contact) a proto pro dostatečně malé hodnoty míry růstu (£ nebo ľ) vychází obě míry (£ a ln(l + t)) přibližně stejně. 4.2.0.12 2. Inverzní problém: 4.2.0.13 Nyní naopak předpokládejme, že známe míru inflace za jednotku času v každém bodě t G ffi. 4.2.0.14 (21) is differential equation for the state function /. The solution is f(x) = eUoln^+1^s) -C (22) 43 Present Value in the Case of General Interest Rate and General Rate of Inflation We determine the value of C substituting constant function for i,: í i—>■ H. It have to be fulfilled that: /(aO = (S + l)*./(0) (23) but f(x) = e(Ioln{E+1)ds) ■ C = C • (S + l)x (24) then: C = /(0) (25) And the state function has the form: /:^e(rin(E(s)+1H /(0) (26) where an i(t) is interest rate per unit of time in time t. 4.2.0.15 2.5. Example 4.2.0.1 6 Let us suppose, that we can check the rate of inflation per unit of time in every moment (it can be the result of measuring and some general theory). We know, that in time 0 the rate of unit time inflation was 0.1 and in time 1 the rate of unit time inflation was 0.2. 4.2.0.1 7 Předpokládejme, že uvnitř intervalu (0,1) se hodnoty inflace za jednotku času měnily podle jedné z následujících možností (skokem v jednom y obou krajních bodů, lineárně, kvadraticky): in this four following ways: _ , , í 0.1, if u < 1 _ , . M2,ni -is M,ni -is / O-1' if m < 0 ,ovs Li{^'-=\Q.2, ifU>l' jo+0-1, iz(u):=-+0.1, h(u) := ^ ifu>q (27) Míra inflace za dobu (0,1) byla v každém z těcto čtyř případů jiná. In general, we have i((0,l» = eUo^1+l^du) -1 (28) So, in our four cases we obtain: i = i! = 0.1: X(l) = eCío1"^-1) dM) - 1 = 0.1 (29) I = ij, = u ~ ^ u2 + 0.1: X(l) = e0"o Mi-i+«2/io) d") _ 1 = 0.132945354 ... (30) l=l3 = u^ l-u + o.i: X(l) = e0"o 1»(1-1+«/io) du) _ 1 = Q.Í49637533 ... (31) i = r4 = 0.2: X(l) = eCío1"^-2) d") - 1 = 0.2 (32) 4.2.0.1 8 3. Interest rate: Předpokládejme, že v každém okamžiku času t známe nominální hodnotu kapitálu f(t), jehož množství roste nebo klesá pouze úročením. Potom v souladu s obvyklou terminologií nazýváme relative increment of the capital during tne time tne interest rate xi: í((ÍI,a)).£!M_£(M m a tento vztah je stejný jako vztah (2), protože míra inflace je úroková míra, kterou se úročí ceny. 4.2.0.19 Je celkem obvyklé, že banka vyhlašuje úrokovou sazbu p. a., i když ji mění častěji, než jednou za rok. Je obvyklé, že střadatelé uloží peníze a chtějí je vybrat podle svých potřeb jindy, než po nějaké předem stanovené periodě času a banky jim to umožní. Ve všech těchto případech se stanovuje úrok jako odměna za dočasné přenechání práva disponovat s penězy zpravidla takto: Jsou-li li úrokové sazby p. a. konstantní na intervalech li = (U, íi+i) ;i = 0...n. 71-1 n(i+-fo(ti+i_ti) (34) i=0 (a úrok s touto mírou je připsán v okamžiku tn). Poznamenejme ještě, že jednotkou měření časuje časový interval na který je normovaná úroková míra (v našem případě je p. a., tedy jednotkou času je rok) a že banky obvykle zaokrouhlují časové údaje na celé dny, což může být jev pouze dočasný. 44 Present Value in the Case of General Interest Rate and General Rate of Inflation 4.2.0.20 V obchodních vztazích je obvyklé, že se úroková míra stanovuje jako po částech konstantní funkce a její hodnota se normuje, ale v obecných makroekonomických a ekonometrických úvahách, pokud například půjde o agregátní úrokovou sazbu, to již po částech konstantní funkce být nemusí. 4.2.0.21 V učebnicích finanční matematiky bývá definováno tzv. spojité úročení, kde je hodnota stavové funkce x v závislosti na úrokové míře £ dána vztahem: a:(í)=a;(0)ei'oí(s)ds (35) namísto vztahu x{t)=x{Q)e$lln(1+i(s))ds (36) jako v (22). Funkce 35 je řešením rovnice £>(*)(*)=£(*)■*(*) (37) namísto rovnice D{x){ť) =ln(l+ £(*)) -x{ť) jako v (21). Funkce (35) ale není zobecněním složeného úročení (33) protože v ustáleném stavu nedává výsledek x(t) =x(0) -(1+0*, (38) proto j e například nelze použít pro kvantitativní výpočty na datech v případě inflace a v příkladě (29) by nám vyšlo J0o.idu _ 1 = 0 105170918 v příkladě (30) by nám vyšlo e/o1i/io«2+o.id« _ 1 = 0.142630812, v příkladě (31) by nám vyšlo J\/iOu+o.idu _ ľ = Q 161834243 a v příkladě (32) by nám vyšlo e/o °'2dM _ i = 0.221402758 a první a poslední výsledekjsou očividně nesprávné. Lze tedy říci, že inflace se nechová jako tzv spojité úročení. 4.2.0.22 3.6. Theorem: The formula (38) je speciálním případem formule (22) pro konstantní úrokovou míru a formule (34) je speciálním případem formule (22) pro po částech konstantní úrokovou míru. 4.2.0.23 Proof: Let us suppose, that i(í) = / is constant: eCío ln(1+I) dM) = e(*'ln(i+i)) = eMi+i)') = (i + if (39) q. e. d. 4.2.0.24 Now let us suppose, that b is piecewise constant and that it has value Ii in every point of the interval Ii = (ti, íi+i); i = 0 ... n. Let \A be the characteristic function of the set A, then b(t) = X(t,,tl+1) ■ h = e(/oln(1+X(*i.*i+i)-í* e-io v K ' je spojitý v topologii stejnoměrné konvergence. Důležitos věty 38 plyne z následujícího lemmatu: 4.2.0.26 3.7. Lemma: pro každou spojitou funkci l definovanou na uzavřeném intervalu existuje posloupnost piecewise konstant funkcí £j takových, že ^ —> l 45 Present Value in the Case of General Interest Rate and General Rate of Inflation 4.2.0.27 Proof: Předpokládejme, že / je spojitá. Podle předpokladu je Dom(/) kompaktní. Zvolíme nějaké e kladné. Pro každé x G Dom(/) najdeme okolí 0(x) takové, aby F(0(x)) C Oe/2{f{x)). 0(x) tvoří pokrytí Dom(/) vybereme konečné podpokrytí fž. Definujeme S = min(yefi(Diam(č7)), kde (Diam(ř7)) je průměr množiny U. Rozdělíme Dom(/) na n disjunktních podintervalů (J?)j = 1™ délky S. Pro každý interval Jf vybereme nějaký bod xÍ7 který leží uvnitř něj a označíme íji = f(xi). Definujeme: x G Jf Ce{%) = Ví- Potom \/x G Dom(/): \f(x) — Cf-ix)\ < £■ A pro = jsou splněny předpoklady lemma. 4.2.0.28 3.8. Corollary: (38) říká, že vzorec (34) je speciálním případem vorce (22). Pro funkci, která má jen konečně mnho bodů nespojitosti můžeme provést aproximaci na každém intervalu, na kterém je spojitá jako v lemmatu (40). A tedy vzorec (22) je limitním případem vzorce (34). 4.2.0.29 Následuje několik aplikací vzorce (22): 4.2.0.30 4. Boundary interest rate: If the interest rate is constant and positive, the state function is increasing and constant. If it is positive, but very quickly decreasing, the state function will be concave. What interest rate makes the state function afine (polynomial of degree 1)? 4.2.0.31 State function is t » x{0)eSo lnd+«(s))^ (41) its derivative is t » a;(0) ln(l + £{t))etiln(1+«(s»ds (42) and the second derivative is x(0)eloln{1+^ds Um + (ln(l + m))2 + (ln(l + £(i)))2 £(*)) 1 »---;-77-,-- (43) We are looking for the interest rate, which makes the second derivative equal to zero. If a;(0) 7^ 0 and £ > 0, second derivative is equal to zero for such a £, which are the solutions of differential equation + (ln(l + £(í)))2 + (ln(l + £(í)))2 £(í) = 0 (44) i. e. we solve the differential equation |e(í) = -((in(i+e(í)))2(i+e(í)). (45) Its solution fulfilled the algebraical equation -(ln(l + e(í)))-1+í = C (46) hence it is any of the functions t = e(^) - 1 (47) for all constant C. We express C using the initial condition 1 ln(l+£(0)) and we have £(t) = e Vt+1»(1+«(°» / - 1 = (1 + ^(0)) (tm(i+e(o))+1) _ i (48) If we put this interest rate into the rule for interesting () we obtain: /' In ( (l+«(0)) Ci»(i+«(0))+i) x(t)=x(0)e\° V / / =a;(0)(tln(l + C(0)) + l) (49) and this really is an affine function. We can conclude: If the interest rate has in the time 0 value £(0) and if on the dependence on the time grows more quickly (decrease slowly) than function t ^ (1 +£(0))(^"1»(1+«(°»+1) -1 (50) 46 Present Value in the Case of General Interest Rate and General Rate of Inflation the state function is convex. If it decreases more quickly, the state function is concave. 4.2.0.32 5. Example with the real data of index of prices of non-regulated prices in Czech republic:. 4.2.0.33 We choose the year as a unit of time. We measured the index of non-regulated prices in every moment P= {1993+ j^}, where i is a natural number less or equal to 120. The values are relative. As a result of the measurement, CPI is a function cutting the points: (1993+^; 91.76), (1993+ \; 92.97), (1993+ \; 93.54), (1993+± 93.94), (1993+ ^; 94.32), (1993 + §; 94.62), (1993+ ^; 95.42), (1993+ f; 96.10), (1993 +1; 97.53), (1993+1 98.62), (1993 + i|; 99.17), (1994; 100.0), (1994 + ±; 100.76), (1994 + i; 101.12), (1994 + ±; 101.40), (1994 +1 101.88), (1994 + ^; 102.27), (1994 + §; 103.45), (1994 + ^; 103.81), (1994 + §; 104.62), (1994 + | 106.13), (1994+|; 107.45), (1994+^; 108.46), (1995; 109.23), (1995+^; 110.57), (1995+^; 111.62), (1995+| 111.97), (1995 + §; 112.72), (1995 + ^; 113.26), (1995 + \; 114.13), (1995 + ^; 113.47), (1995 + | 113.41), (1995+|; 114.41), (1995+f; 115.20), (1995+^; 116.09), (1996; 116.82), (1996+^; 118.97), (1996+J 119.62), (1996 + \; 120.43), (1996 + \; 121.15), (1996 + ^; 121.96), (1996 + \; 123.06), (1996 + ^ 123.18), (1996+f; 122.67), (1996+f; 123.04), (1996+|; 123.76), (1996+^; 124.44), (1997; 125.23), (1997+^ i; 126.71), (1997 + \; 126.85), (1997 + ±; 127.45), (1997 + ^ J_- no^ nnnv i 2. nnn? i 3. m nnnv i o 1 I 126.28), (1997 + ±; 126.71), (1997 + \; 126.85), (1997 + ±; 127.45')", (1997 + ^; 127.60), (1997 +'| 129.42), (1997 + ^; 129.67), (1997 + §; 130.77), (1997 + f; 131.58), (1997 + |; 132.33), (1997 + ±| 133.01), (1998; 133.75), (1998+^; 135.76), (1998+^; 136.71), (1998+^; 136.85), (1998+|; 137.12), (1998+^ 137.26), (1998 + \; 137.81), (1998 + ^; 137.53), (1998 + §; 137.12), (1998 + f; 137.25), (1998 + | 136.84), (1998+^; 136.43), (1999; 136.02), (1999+^; 136.70), (1999+±; 136.57), (1999+|; 136.29), (1999+| 136.84), (1999 + ^; 136.70), (1999 + \; 136.98), (1999 + ^; 136.98), (1999 + |; 137.11), (1999 + | 136.98), (1999+f; 136.98), (1999+±±; 137.39), (2000; 138.21), (2000+^; 139.04), (2000+±; 139.32), (2000+] 139.32), (2000 + i; 139.32), (2000 + ^; 139.74), (2000 + \; 140.71), (2000 + ^; 141.42), (2000 + | _5. 1 a 1 qq\ r U < tj are the moments, where CPI is measured, then we have to fulfil the equation: /t.+i ln(l + i(t))dt = lnoCPI(t„+i) -lnoCPI(i„) (54) and the accuracy of approximation should be measured as the measure of unreliableness of the equation. 4.2.0.36 More precisely: Let us choose the function: f. x i—> 0.04944210 + 0.009094682 cos (xn) + 0.01215932 sin (2 xn) + +0.01823419 sin (4 xtt) + 0.05263253 sin (1/3 xtt) - (55) -0.009667420 cos (1/3 aw) - 0.07212055 cos (1/4 aw) + 0.02834145 sin (1/2 aw) . The question is, how accurately this function corresponds to the measured values. The interesting thing is, that the function does not approximate the measured values, but some other values, which depend on the measured ones. And the attitude towards this two values is given by (54). So if we can compare the exactness of approximation with the same measure as we did in the previous cases, we have to compute the number vtf,ta(i+iW)di crimpy h v " "cpwJ (56) and this is the only relevant measure. In our case we obtain: 4.9914 We used the simplest function for the approximation of rate of inflation we considered yet, and it gives us far the best approximation we have obtained! 48 4.2.0.37 References: [I] Aczél, J.: Lectures on functional Equations and Their Applications, New York, Academic Press 1966 [2] Barro, Robert J., Sala-i-martin, Xavier: Economic Growth, MIT Press Cambridgr, Massachusetts, London, England 1999, ISBN 0-262-02459-4 [3] Dieudonné, J.: Treatise on analysis, Vol. Ill, Pure and Applied Mathematics, Vol. 10-111, Academic Press. New York-London, 1972 [4] Dupačová, J., Hurt, J., Stepán, J.: Stochastic Modeling in economics andFinance, Kluwert Academic Publishers, 2002, ISBN 1-4020-0840-6. [5] Studený V: Functional Equation of the Rate of Inflation, e-print archive of Coronell University, 2003, www.arxiv.org [6] 6. Conclusion: [7] Our sketch brings an open problem: what does the space of solution of (54) look like and what is the space of functions, the rate of inflation lies in. How can it help us to understand what kind of measure is this one which we called rate (of interest or inflation...). It shows us the possibility to approximate the measure of inflation except of the index of prices and it suggests, that the solution should be better. [8] 7. Apendix: Maple programs for computation. Let us suppose that in iota we have the function for approximation, in our case: [9] > iota:=.4944210e-l + .9094682e-2 * cos(x*Pi) + .1215932e-l * sin(2*x*Pi) + .1823419e-1 * sin(4*x*Pi)+ .5263253e-l * sin(l/3*x*Pi) - .9667420e-2 * cos(l/3*x*Pi) - .7212055e-1 * cos(l/4*x*Pi)+ .2834145e-l * sin(l/2*x*Pi); [9] Let us suppose, that the values of CPI are in the file Values and the corresponding moments of time, in which the values are measured, are in file Points. Then we can compute the preciseness of approximation using the following program: [10] > Mistake := 0; [II] > for i to NumberOfPoints do [12] > A := evalf(int(ln(l+iota),x = Points[1] .. Points[i+1])); [13] > A := Values[l]*exp(A) ; [14] > Mistake := Mistake+(A-Values[i+1] )~2 [15] > end do: [16] > print(Sum('Delta'~2,i)=Mistake); [1] References: [1] Aczél, J.: Lectures on functional Equations and Their Applications, New York, Academic Press 1966 [2] Dieudonné, J.: Treatise on analysis, Vol. Ill, Pure and Applied Mathematics, Vol. 10-111, Academic Press. New York-London, 1972 [3] Dupačová, J., Hurt, J., Stepán, J.: Stochastic Modeling in economics andFinance, Kluwert Academic Publishers, 2002, ISBN 1-4020-0840-6. [4] Studený V: Functional Equation of the Rate of Inflation, e-print archive of Coronell University, 2003, www.arxiv.org 49