|
|
| Line 1: |
Line 1: |
| {{Notability|date=June 2010}}
| | Surely the second option would be more beneficial for any website. Medical word press themes give you the latest medical designs. These templates are professionally designed and are also Adsense ready. If you are using videos on your site then this is the plugin to use. It's as simple as hiring a Wordpress plugin developer or learning how to create what is needed. <br><br>As you know today Word - Press has turn out to be a tremendously popular open source publishing and blogging display place. After all, Word - Press is free, many of the enhancements for Word - Press like themes and plugins are also free, and there is plenty of free information online about how to use Word - Press. There are number of web services that offer Word press development across the world. This is identical to doing a research as in depth above, nevertheless you can see various statistical details like the number of downloads and when the template was not long ago updated. W3C compliant HTML and a good open source powered by Word - Press CMS site is regarded as the prime minister. <br><br>The entrepreneurs can easily captivate their readers by using these versatile themes. After sending these details, your Word - Press blog will be setup within a few days. If you liked this information and you would certainly such as to get even more information regarding [http://www.twinurl.com/wordpressbackup91955 wordpress dropbox backup] kindly browse through our own web site. Use this section to change many formatting elements. Newer programs allow website owners and internet marketers to automatically and dynamically change words in their content to match the keywords entered by their web visitors in their search queries'a feat that they cannot easily achieve with older software. Premium vs Customised Word - Press Themes - Premium themes are a lot like customised themes but without the customised price and without the wait. <br><br>A built-in widget which allows you to embed quickly video from popular websites. This plugin allows a webmaster to create complex layouts without having to waste so much time with short codes. Thus it is difficult to outrank any one of these because of their different usages. IVF ,fertility,infertility expert,surrogacy specialist in India at Rotundaivf. Where from they are coming, which types of posts are getting top traffic and many more. <br><br>Every single module contains published data and guidelines, usually a lot more than 1 video, and when pertinent, incentive links and PDF files to assist you out. An ease of use which pertains to both internet site back-end and front-end users alike. As a result, it is really crucial to just take aid of some experience when searching for superior quality totally free Word - Press themes, Word - Press Premium Themes for your web site. ) Remote Login: With the process of PSD to Wordpress conversion comes the advantage of flexibility. Definitely when you wake up from the slumber, you can be sure that you will be lagging behind and getting on track would be a tall order. |
| [[File:Compound Interest with Varying Frequencies.svg|thumb|right|350px|The effect of earning 20% annual interest on an initial $1,000 investment at various compounding frequencies]]
| |
| | |
| Analogous to [[continuous compounding]], a continuous annuity<ref>{{Cite book
| |
| |ref = Dictionary
| |
| |last = James
| |
| |first = Robert C
| |
| |coauthors = James, Glen
| |
| |year = 1992
| |
| |title = Mathematics Dictionary
| |
| |url = http://books.google.com/books?id=UyIfgBIwLMQC
| |
| |publisher = Chapman and Hall
| |
| }} - Entry on [http://books.google.com/books?id=UyIfgBIwLMQC&lpg=PA86&dq=Continuous%20annuity&pg=PA86#v=onepage&q=Continuous%20annuity&f=false continuous annuity]</ref><ref>[[#Dictionary|Mathematics Dictionary]] p.86</ref> is an [[Annuity_(finance_theory)#Ordinary_annuity|ordinary annuity]] in which the payment interval is narrowed indefinitely. A (theoretical) '''continuous repayment mortgage''' is a mortgage loan paid by means of a continuous annuity.
| |
| | |
| Mortgages (i.e., mortgage loans) are generally settled over a period of years by a series of fixed regular payments commonly referred to as an [[Annuity (finance theory)|annuity]]. Each payment accumulates [[compound interest]] from time of deposit to the end of the mortgage timespan at which point the sum of the payments with their accumulated interest equals the value of the loan with interest compounded over the entire timespan. Given loan ''P''<sub>0</sub>, per period interest rate i, number of periods ''n'' and fixed per period payment ''x'', the end of term balancing equation is:
| |
| | |
| ::<math>P_0(1+i)^{n} = \sum_{k=1}^n x(1+i)^{n-k}=\frac{x[(1+i)^n - 1]}{i}</math>
| |
| | |
| Summation can be computed using the standard formula for summation of a [[geometric sequence]].
| |
| | |
| In a (theoretical) continuous-repayment mortgage the payment interval is narrowed indefinitely until the discrete interval process becomes continuous and the fixed interval payments become—in effect—a literal cash "flow" at a fixed annual rate. In this case, given loan ''P''<sub>0</sub>, annual interest rate ''r'', loan timespan ''T'' (years) and annual rate ''M''<sub>''a''</sub>, the [[infinitesimal]] cash flow elements ''M''<sub>''a''</sub>''δt'' accumulate [[continuously compounded interest]] from time t to the end of the loan timespan at which point the balancing equation is:
| |
| | |
| ::<math>P_0e^{rT}=\int\limits_{0}^{T} M_ae^{r(T-t)}\, dt=\frac{M_a(e^{rT}-1)}{r}.</math>
| |
| | |
| Summation of the cash flow elements and accumulated interest is effected by integration as shown. It is assumed that compounding interval and payment interval are equal—i.e., compounding of interest always occurs at the same time as payment is deducted.<ref>Strictly speaking compounding occurs momentarily before payment is deducted so that interest is calculated on the balance as it was before deduction of the period payment.</ref>
| |
| | |
| Within the timespan of the loan the time continuous mortgage balance function obeys a first order [[linear differential equation]] (LDE)<ref>[[#CFP|Beckwith]] p. 116: ''"Technically speaking, the underlying equation is known as an ordinary, linear, first order, inhomogenous, scalar differential equation with a boundary condition."''</ref> and an alternative derivation thereof may be obtained by solving the LDE using the method of [[Laplace transforms]].
| |
| | |
| Application of the equation yields a number of results relevant to the financial process which it describes. Although this article focuses primarily on mortgages, the methods employed are relevant to any situation in which payment or saving is effected by a regular stream of fixed interval payments (annuity).
| |
| | |
| ==Derivation of time-continuous equation==
| |
| The classical formula for the present value of a series of ''n'' fixed monthly payments amount ''x'' invested at a monthly interest rate ''i''% is:
| |
| | |
| :<math>P_v(n) = \frac{x(1 - (1 + i)^{-n})}{i}</math>
| |
| | |
| The formula may be re-arranged to determine the monthly payment ''x'' on a loan of amount ''P''<sub>0</sub> taken out for a period of ''n'' months at a monthly interest rate of ''i''%:
| |
| | |
| :<math>x = \frac{P_0\cdot i}{1 - (1 + i)^{-n}}</math>
| |
| | |
| We begin with a small adjustment of the formula: replace ''i'' with ''r''/''N'' where ''r'' is the annual interest rate and ''N'' is the annual frequency of compounding periods (''N'' = 12 for monthly payments). Also replace ''n'' with ''NT'' where ''T'' is the total loan period in years. In this more general form of the equation we are calculating ''x''(''N'') as the fixed payment corresponding to frequency ''N''. For example if ''N'' = 365, ''x'' corresponds to a daily fixed payment. As ''N'' increases, ''x''(''N'') decreases but the product ''N''·''x''(''N'') approaches a limiting value as will be shown:
| |
| | |
| :<math>x(N) = \frac{P_0\cdot r}{N(1 - (1 + \frac{r}{N})^{-NT})}</math>
| |
| | |
| :<math>N\cdot x(N) = \frac{P_0\cdot r}{1 - (1 + \frac{r}{N})^{-NT}}</math>
| |
| | |
| Note that ''N''·''x''(''N'') is simply the amount paid per year – in effect an annual repayment rate ''M''<sub>''a''</sub>.
| |
| | |
| It is well established that:
| |
| | |
| :<math>\lim_{N\to\infty}\left(1+\frac{r}{N}\right)^{Nt}=e^{rt}</math> <ref>[[#CFP|Beckwith]] p.115</ref><ref>[[#A&T|Munem and Foulis]] p.273</ref>
| |
| | |
| Applying the same principle to the formula for annual repayment, we can determine a limiting value:
| |
| | |
| :<math>M_a=\lim_{N\to\infty}N\cdot x(N)=\lim_{N\to\infty}\frac{P_0\cdot r}{1 - (1 + \frac{r}{N})^{-NT}}=\frac{P_0\cdot r}{1 - e^{-rT}}. </math> <ref>[[#CFP|Beckwith:]] Equation (29) p. 123.</ref>
| |
| | |
| At this point in the orthodox formula for present value, the latter is more properly represented as a function of annual compounding frequency ''N'' and time ''t'':
| |
| | |
| :<math>P_v(N,t)=\frac{N\cdot x(N)(1 - (1 + \frac{r}{N})^{-Nt})}{r}</math>
| |
| | |
| Applying the limiting expression developed above we may write present value as a purely time dependent function:
| |
| | |
| :<math>P_v(t)=\lim_{N\to\infty}P_v(N,t)=\frac{M_a}{r}(1-e^{-rt})</math><ref>See also: {{Cite book
| |
| |ref = USMAG
| |
| |last = Wisdom
| |
| |first = John C
| |
| |coauthors = Hasselback, James R.
| |
| |year = 2008
| |
| |title = U.S. Master Accounting Guide 2008.
| |
| |url = http://books.google.com/books?id=l1_FbS0Gj9EC
| |
| |publisher = C C H Inc 2008.
| |
| }} ps. 470–471</ref>
| |
| | |
| [[File:mortgage balance.jpg|thumb|300px|right|Figure 1]]Noting that the balance due ''P''(''t'') on a loan ''t'' years after its inception is simply the present value of the contributions for the remaining period (i.e. ''T'' − ''t''), we determine:
| |
| | |
| :<math>P(t) = \frac{M_a}{r}(1 - e^{-r(T-t)})=\frac{P_0(1 - e^{-r(T-t)})}{1 - e^{-rT}}</math> <ref>[[#CFP|Beckwith:]] Equation (31) p. 124.</ref>
| |
| | |
| The graph(s) in the diagram are a comparison of balance due on a mortgage (1 million for 20 years @ ''r'' = 10%) calculated firstly according to the above time continuous model and secondly using the Excel PV function. As may be seen the curves are virtually indistinguishable – calculations effected using the model differ from those effected using the Excel PV function by a mere 0.3% (max). The data from which the graph(s) were derived can be viewed [[Media:Mortgage Balance model data.pdf|here.]]
| |
| | |
| ==Comparison with similar physical systems==
| |
| | |
| Define the "reverse time" variable ''z'' = ''T'' − ''t''. (''t'' = 0, ''z'' = ''T'' and ''t'' = ''T'', ''z'' = 0). Then:
| |
| [[File:RC-circuit2.png|right|200px|thumb| Plotted on a time axis normalized to system time constant ('''τ''' = 1/''r'' years and '''τ''' = ''RC'' seconds respectively) the mortgage balance function in a CRM (green) is a mirror image of the step response curve for an RC circuit (blue).The vertical axis is normalized to system asymptote i.e. perpetuity value M<sub>a</sub>/r for the CRM and applied voltage V<sub>0</sub> for the RC circuit.]]
| |
| | |
| : <math>P(z) = \frac{M_a}{r}(1 - e^{-rz}).</math>
| |
| | |
| This may be recognized as a solution to the "reverse time" differential equation:
| |
| | |
| : <math>\frac{M_a}{r} = \frac{1}{r}\frac{dP(z)}{dz} + P(z).</math>
| |
| | |
| Electrical/electronic engineers and physicists will be familiar with an equation of this nature: it is an exact analogue of the type of differential equation which governs (for example) the charging of a capacitor in an RC circuit.
| |
| | |
| : <math>V_0 = RC\frac{dV(t)}{dt}+V(t).</math>
| |
| | |
| The key characteristics of such equations are explained in detail at [[RC_circuit#Time-domain_considerations|RC circuits]]. For home owners with mortgages the important parameter to keep in mind is the [[time constant]] of the equation which is simply the reciprocal of the annual interest rate ''r''. So (for example) the time constant when the interest rate is 10% is 10 years and the period of a home loan should be determined – within the bounds of affordability – as a minimum multiple of this if the objective is to minimise interest paid on the loan.
| |
| | |
| ==Mortgage difference and differential equation==
| |
| | |
| The conventional [[difference equation]] for a mortgage loan is relatively straightforward to derive - balance due in each successive period is the previous balance plus per period interest less the per period fixed payment.
| |
| | |
| Given an '''annual''' interest rate ''r'' and a borrower with an '''annual''' payment capability ''M''<sub>''N''</sub> (divided into N equal payments made at time intervals Δ''t'' where Δ''t'' = 1/''N'' years), we may write:
| |
| | |
| :<math>\begin{align}
| |
| P_{t+\Delta t} & = P_t+(rP_t-M_N)\Delta t\\[12pt]
| |
| \dfrac{P_{t+\Delta t}-P_t}{\Delta t} & = rP_t-M_N
| |
| \end{align}</math>
| |
| | |
| If ''N'' is increased indefinitely so that Δ''t'' → 0, we obtain the continuous time differential equation:
| |
| | |
| :<math>{\operatorname{d}P(t)\over\operatorname{d}t}=rP(t)-M_a </math> <ref>[[#CFP|Beckwith:]] Equation (25) p. 123</ref><ref>[[#GIT|Hackman:]] Equation (2) p.1</ref>
| |
| | |
| Note that for there to be a continually diminishing mortgage balance, the following inequality must hold:
| |
| | |
| :<math>P_0 \leqslant \frac{M_a}{r} </math><ref>Where equality holds, the mortgage becomes a [[perpetuity]].</ref>
| |
| | |
| ''P''<sub>0</sub> is the same as ''P''(0) – the original loan amount or loan balance at time ''t'' = 0.
| |
| | |
| ===Solving the difference equation===
| |
| We begin by re-writing the difference equation in recursive form:
| |
| | |
| :<math>P_{t+\Delta t}=P_t(1+r\Delta t)-M_N\Delta t\; </math>
| |
| | |
| Using the notation ''P''<sub>''n''</sub> to indicate the mortgage balance after ''n'' periods, we may apply the recursion relation iteratively to determine ''P''<sub>1</sub> and ''P''<sub>2</sub>:
| |
| | |
| :<math>P_1 = P_0(1+r\Delta t)-M_N\Delta t\;</math>
| |
| | |
| :<math>\begin{align}
| |
| P_2 &= [P_0(1+r\Delta t)-M_N\Delta t](1+r\Delta t)-M_N\Delta t\\
| |
| &= P_0(1+r\Delta t)^2 - M_N\Delta t(1+r\Delta t)-M_N\Delta t
| |
| \end{align}</math>
| |
| | |
| It can already be seen that the terms containing ''M''<sub>''N''</sub> form a geometric series with common ratio 1 + ''r''Δ ''t''. This enables us to write a general expression for ''P''<sub>''n''</sub>:
| |
| | |
| :<math>\begin{align}
| |
| P_n&=P_0(1+r\Delta t)^n-\sum_{k=1}^{n} M_N\Delta t(1+r\Delta t)^{n-k} \\
| |
| &=P_0(1+r\Delta t)^n-\dfrac{M_N\Delta t[(1+r\Delta t)^n - 1]}{r\Delta t}
| |
| \end{align}</math>
| |
| | |
| Finally noting that ''r'' Δ ''t'' = ''i'' the per-period interest rate and <math>M_N\Delta t=x</math> the per period payment, the expression may be written in conventional form:
| |
| | |
| :<math>P_n=P_0(1+i)^n-\dfrac{x[(1+i)^n - 1]}{i}</math>
| |
| | |
| If the loan timespan is m periods, then ''P''<sub>''m''</sub> = 0 and we obtain the standard present value formula: | |
| | |
| :<math>P_0=\dfrac{x[1-(1+i)^{-m}]}{i}</math>
| |
| | |
| ===Solving the differential equation===
| |
| :<math>{\operatorname{d}P(t)\over\operatorname{d}t}=rP(t)-M_a </math>
| |
| | |
| One method of solving the equation is to obtain the [[Laplace transform]] ''P''(''s''):
| |
| | |
| :<math>P(s)=\frac{M_a}{s(r-s)}=\frac{M_a}{r} \times \frac{(-r)}{s(s-r)}.</math>
| |
| | |
| Using a [[Laplace_transform#Table_of_selected_Laplace_transforms|table of Laplace transforms]] and their time domain equivalents, ''P''(''t'') may be determined:
| |
| | |
| :<math>P(t)=\frac{M_a}{r}(1-e^{rt}).</math>
| |
| | |
| In order to fit this solution to the particular start and end points of the mortgage function we need to introduce a time shift of ''T'' years (''T'' = loan period) to ensure the function reaches zero at the end of the loan period:
| |
| | |
| :<math>
| |
| \begin{align}
| |
| & P(t) = \frac{M_a}{r}(1-e^{r(t-T)}) \\[8pt]
| |
| & P_0 = \frac{M_a}{r}(1-e^{-rT})\Rightarrow\frac{M_a}{r}=\frac{P_0}{1-e^{-rT}} \\[8pt]
| |
| \Rightarrow & P(t) = \frac{P_0(1-e^{-r(T-t)})}{1-e^{-rT}}
| |
| \end{align}
| |
| </math>
| |
| | |
| Note that both the original solution and "time-shifted" version satisfy the original differential equation whence both are derived.
| |
| | |
| Similar to the expression derived above for ''P''<sub>''n''</sub> in the difference equation, the expression for ''P''(''t'') may be written in the following algebraically equivalent form:
| |
| | |
| :<math>P(t) = P_0e^{rt} - \frac{M_a}{r}(e^{rt}-1)</math> | |
| | |
| ==Calculation of accumulated interest and principal payments==
| |
| | |
| Re-arranging the original differential equation we obtain:
| |
| | |
| :<math>M_a = rP(t) - \frac{dP(t)}{dt}</math>
| |
| | |
| Integrating both sides of the equation yields:
| |
| | |
| :<math>M_a.t = \int_0^t rP(t)\,dt \, - \int_0^t \frac{dP(t)}{dt}\,dt \, </math>
| |
| | |
| The first integral on the right hand side determines the accumulated interest payments from time of inception to time t whilst the second determines the accumulated principal payments over the same period. The sum of these interest and principal payments must equal the cumulative fixed payments at time ''t'' i.e. ''M''<sub>''a''</sub>''t''. Evaluating the first integral on the right we obtain an expression for ''I''(''t''), the interest paid:
| |
| | |
| :<math>I(t)=M_at-\frac{M_a(e^{rt}-1)}{re^{rT}}</math>
| |
| | |
| Unsurprisingly the second integral evaluates to ''P''<sub>0</sub> − ''P''(''t'') and therefore:
| |
| | |
| :<math>I(t)=M_at-P_0+P(t) \,</math>
| |
| | |
| The reader may easily verify that this expression is algebraically identical to the one above.
| |
| | |
| ==Loan cost factor==
| |
| | |
| The cost of a loan is simply the annual rate multiplied by loan period:
| |
| | |
| :<math> C = M_aT= \frac{P_0 rT}{1-e^{-rT}} </math>
| |
| | |
| Let ''s'' = ''rT''. Then we may define loan cost factor ''C''(''s'') such that ''C'' = ''P''<sub>0</sub>''C''(s) i.e.: ''C''(''s'') is the cost per unit of currency loaned.
| |
| | |
| :<math> C(s) = \frac{rT}{1-e^{-rT}}=\frac{s}{1-e^{-s}} </math>
| |
| | |
| The function ''C''(''s'') is characterised by having a limiting value of 1 when ''s'' is close to zero since for small values of ''s'', exp(−''s'') ≈ 1 − ''s'' and the denominator simplifies to ''s''. Also when ''s'' is very large, exp(−''s'') is small so ''C''(''s'') ≈ ''s'' and thus loan cost ''C'' ≈ ''P''<sub>0</sub>''rT'' (''rT'' >> 0).
| |
| | |
| By way of example, consider a loan of 1000000 at 10% repaid over 20 years. Then ''s'' = 0.1 × 20 = 2.
| |
| | |
| :<math> C = 1000000 \times \frac{2}{1-e^{-2}} \approx 2.313\times10^6 </math>
| |
| | |
| The product rT is an easily obtained but important parameter in determining loan cost according to the equation C=P<sub>0</sub>xC(s). This is best illustrated by plotting the cost factor function for s values in domain [0;5]. The linear behaviour of the function for higher values of ''s'' is clear.
| |
| | |
| ==Equivalent simple interest cost factor==
| |
| For a fixed term loan of t years, we may compare the above loan cost factor against an equivalent simple interest cost factor ''1+s<sub>e</sub>'' where ''s<sub>e</sub>=r<sub>e</sub>t'' and ''r<sub>e</sub>'' is the equivalent simple interest rate:
| |
| | |
| :<math>\frac{s}{1-e^{-s}}=1+s_e </math>
| |
| | |
| It is straightforward to determine ''s<sub>e</sub>'' in terms of s. Dividing by loan time period t will then give the equivalent simple interest rate. More challenging is the reverse determination of s given ''s<sub>e</sub>''.
| |
| | |
| In his book ''Problem Solving with True Basic'',<ref>[[#BDH|Hahn]] p. 247</ref> Dr B.D. Hahn has a short section on certain 'hire purchase' schemes in which ''interest is calculated in advance in one lump sum, which is added to the capital amount, the sum being equally divided over the repayment period. The buyer, however, is often under the impression that the interest is calculated on a reducing balance.''
| |
| | |
| The above example is adapted from the one given in Dr Hahn's book in which he employs the Newton-Raphson algorithm to solve the same problem albeit for a discrete interval (i.e. monthly) repayment loan over the same time period (3 years). As with many similar examples the discrete interval problem and its solution is closely approximated by calculations based on the continuous repayment model - Dr Hahn's solution for interest rate is 40.8% as compared to the 41.6% calculated above.
| |
| | |
| ==Period of a loan==
| |
| | |
| If a borrower can afford an annual repayment rate ''M''<sub>''a''</sub>, then we can re-arrange the formula for calculating ''M''<sub>''a''</sub> to obtain an expression for the time period ''T'' of a given loan ''P''<sub>0</sub>:
| |
| | |
| :<math>
| |
| \begin{align}
| |
| & M_a = \frac{P_0 r}{1-e^{-rT}} \\[8pt]
| |
| \Rightarrow & T = \frac{1}{r}\ln\frac{M_a}{M_a-P_0 r} = -\frac{1}{r}\ln\left(1 - \frac{P_0 r}{M_a} \right)
| |
| \end{align}
| |
| </math>
| |
| | |
| ==Minimum payment ratio==
| |
| | |
| The minimum payment ratio of a loan is the ratio of minimum possible payment rate to actual payment rate. The minimum possible payment rate is that which just covers the loan interest – a borrower would in theory pay this amount forever because there is never any decrease in loan capital. We will use the letter ''k'' to denote minimum payment ratio:
| |
| | |
| :<math> k = \frac{M_\min}{M_a} = \frac{P_0r}{M_a} </math>
| |
| | |
| Now we may consider a small re-arrangement of the equation for loan period ''T'':
| |
| | |
| :<math> T = -\frac{1}{r}\ln\left(1-\frac{P_0r}{M_a}\right) </math>
| |
| :<math> rT = s(k) = -\ln(1-k) \; </math>
| |
| | |
| Plotting ''s''(''k'') against ''k'' gives a very graphic demonstration of why it is a good idea to keep the ''k'' value well below the asymptote at ''k'' = 1 since in the vicinity thereof, ''s''(''k'') increases sharply and therefore so does loan cost which is in turn an increasing function of parameter ''s'' (''rT'' product).
| |
| | |
| =="Half-life" of a loan==
| |
| | |
| A useful parameter of the mortgage model is the "half-life" of the loan which is the time it takes for the balance on the loan to reach half its original value. To determine the "half-life" we may write:
| |
| | |
| :<math>\frac{P(t)}{P_0}=\frac{1}{2}=\frac{1-e^{-r(T-t)}}{1-e^{-rT}}</math>
| |
| | |
| Solving for ''t'' we obtain:
| |
| | |
| :<math>t_{\frac{1}{2}}=\frac{1}{r}\ln\left(\frac{1+e^{rT}}{2}\right)</math> <ref>[[#CFP|Beckwith:]] Equation (23) p. 122. Beckwith uses this formula in relation to a sinking fund but notes (p.124) that the formula is identical for an amortization process.</ref>
| |
| | |
| For example applying the formula to some test data (loan of 1 million at 10% for 20 years) we obtain the half-life as 14.34 years. If in practice the loan is being repaid via monthly instalments, the decimal portion can be converted to months and rounded so this answer would equate to 172 months.
| |
| | |
| ==Calculation of interest rate==
| |
| | |
| In the discrete time interval model, calculation of a mortgage based interest rate given the remaining parameters has not been possible using analytic methods. Implementations such as the Excel "rate" function employ a numerical "trial and improvement" method to determine interest rate. At first glance this would also seem to be the case for the continuous repayment model. Given:
| |
| | |
| :<math>P_0=\frac{M_a}{r}(1-e^{-rT})</math>
| |
| | |
| we may write:
| |
| | |
| :<math>rP_0=M_a(1-e^{-rT})\,</math>
| |
| :<math>\Rightarrow rP_0-M_a+M_ae^{-rT}=0</math>
| |
| | |
| [[File:Interest rate plot.jpg|thumb|300px|right|Figure 1]]In order to visualise the above as a function of ''r'' (for which we wish to determine zeroes), it will be helpful to select numerical values of ''P''<sub>0</sub>, ''M''<sub>''a''</sub> and ''T'' as 10000, 6000 and 3 respectively and plot as shown at right. It will be noted that the function has a minimum value which can be determined by differentiation:
| |
| | |
| :<math>f'(r)=0\,</math>
| |
| :<math>\Rightarrow P_0-M_aTe^{-rT}=0</math>
| |
| :<math>\Rightarrow r=\frac{1}{T}\ln{\frac{M_aT}{P_0}}</math>
| |
| | |
| Since the function is approximately parabolic between the roots at ''r'' = 0 and the sought value, we may estimate the required root as:
| |
| | |
| :<math>r\approx\frac{2}{T}\ln{\frac{M_aT}{P_0}}</math>
| |
| | |
| Using this as a starting point, increasingly accurate values for the root may be determined by repeated iterations of the [[Newton–Raphson|Newton–Raphson algorithm]]:<ref>[[#CFP|Beckwith:]] (p.125):''"In the determination of rates of interest for given continuous payment schedules, it is frequently necessary to determine the roots of transcendental functions."''. Beckwith details two methods: Successive substitution and Newton–Raphson. (ps. 126–127).</ref>
| |
| | |
| :<math>r_{1} = r_0 - \frac{f(r_0)}{f'(r_0)}.\,\!</math>
| |
| | |
| Some experimentation on [http://www58.wolframalpha.com Wolfram Alpha] reveals that an [http://www.wolframalpha.com/input/?i=solve+for+r+in+P%3Dm%2Fr*%281-exp%28-r*T%29%29 exact analytical solution] employing the [[Lambert W function|Lambert-W]] or "product log" function can be obtained. Setting ''s'' = ''M''<sub>''a''</sub>''T''/''P''<sub>0</sub> we obtain:
| |
| | |
| :<math>r=\frac{1}{T}\left (W(-se^{-s})+s\right )</math>
| |
| | |
| In the region of interest ''W''(−''se''<sup>−''s''</sup>) is a bi-valued function. The first value is just −''s'' which yields the trivial solution ''r'' = 0. The second value evaluated within the context of the above formula will provide the required interest rate.
| |
| | |
| The following table shows calculation of an initial estimate of interest rate followed by a few iterations of the Newton–Raphson algorithm. There is rapid convergence to a solution accurate to several decimal places as may be corroborated against the [[Media:interest rate calc example.pdf|analytical solution]] using the Lambert ''W'' or "productlog" function on Wolfram Alpha.
| |
| | |
| {| border="1" cellpadding="2"
| |
| | align="center" style="background:#f0f0f0;" width="100"|'''Loan''' (''P'')
| |
| | align="center" style="background:#f0f0f0;" width="100"|'''Period''' (''T'')
| |
| | align="center" style="background:#f0f0f0;" width="100"|'''Annual payment rate''' (''Ma'')
| |
| | align="center" style="background:#f0f0f0;" width="100"|'''Initial estimate:''' 2 ln(''MaT''/''P'')/''T''
| |
| |-
| |
| | 10000||3||6000||39.185778%
| |
| |}
| |
| | |
| '''''Newton–Raphson iterations'''''
| |
| | |
| {| border="1" cellpadding="2"
| |
| | align="center" style="background:#f0f0f0;" width="100"|''n''
| |
| | align="center" style="background:#f0f0f0;" width="100"|''r''(''n'')
| |
| | align="center" style="background:#f0f0f0;" width="100"|''f''[''r''(''n'')]
| |
| | align="center" style="background:#f0f0f0;" width="100"|''f''<nowiki>'</nowiki>[''r''(''n'')]
| |
| |-
| |
| | 0||39.185778%||−229.57||4444.44
| |
| |-
| |
| | 1||44.351111%||21.13||5241.95
| |
| |-
| |
| | 2||43.948044%||0.12||5184.06
| |
| |-
| |
| | 3||43.945798%||0||5183.74
| |
| |}
| |
| | |
| ==Present value and future value formulae==
| |
| | |
| Corresponding to the standard formula for the present value of a series of fixed monthly payments, we have already established a time continuous analogue:
| |
| | |
| :<math>P_v(t)=\frac{M_a}{r}(1-e^{-rt}).</math>
| |
| | |
| In similar fashion, a future value formula can be determined:
| |
| | |
| :<math>F_v(t)=\frac{M_a}{r}(e^{rt}-1).</math> <ref>See also: {{Cite book
| |
| |ref = TOF
| |
| |last = King
| |
| |first = George
| |
| |year = 1898
| |
| |title = The Theory of Finance. Being a Short Treatise on the Doctrine of Interest and Annuities-Certain
| |
| |url = http://books.google.com/books?id=ZOWCNpT2cL4C
| |
| |publisher = London: Charles and Edwin Layton. Reprinted March 2010 Nabu Press
| |
| |ISBN = 1-146-31870-7
| |
| }} p. 22. Older actuarial textbooks refer to "interest convertible momently" and "payments momently" when discussing continuous annuities.</ref>
| |
| | |
| In this case the annual rate ''M''<sub>''a''</sub> is determined from a specified (future) savings or sinking fund target ''P''<sub>''T''</sub> as follows.
| |
| | |
| :<math>M_a=\lim_{N\to\infty}N\cdot x(N)=\lim_{N\to\infty}\frac{P_T\cdot r}{(1 + \frac{r}{N})^{NT}-1}=\frac{P_T\cdot r}{e^{rT}-1}. </math> <ref>[[#CFP|Beckwith:]] Equation (19) p. 121.</ref>
| |
| | |
| It will be noted that as might be expected:
| |
| | |
| :<math>F_v(t)=P_v(t) \times e^{rt}. \, </math>
| |
| | |
| Another way to calculate balance due ''P''(''t'') on a continuous-repayment loan is to subtract the future value (at time ''t'') of the payment stream from the future value of the loan (also at time ''t''):
| |
| | |
| :<math>P(t)=P_0 e^{rt}-\frac{M_a}{r}(e^{rt}-1).</math> <ref>[[#CFP|Beckwith:]] Equation (27) p. 123.</ref>
| |
| | |
| ==Example==
| |
| | |
| The following example from a school text book<ref>[[#MJG|Glencross]] p. 67</ref> will illustrate the conceptual difference between a savings annuity based on discrete time intervals (per month in this case) and one based on continuous payment employing the above future value formula:
| |
| | |
| ''On his 30th birthday, an investor decides he wants to accumulate R500000 by his 40th birthday. Starting in one month's time he decides to make equal monthly payments into an account that pays interest at 12% per annum compounded monthly. What monthly payments will he have to make?''
| |
| | |
| For the sake of brevity, we will solve the "discrete interval" problem using the Excel PMT function:
| |
| | |
| :<math>x(12) = PMT(1%, 120, 500000) = 2173.55</math>
| |
| | |
| The amount paid annually would therefore be 26082.57.
| |
| | |
| For a theoretical continuous payment savings annuity we can only calculate an annual ''rate'' of payment:
| |
| | |
| :<math>M_a=\frac{500000 \times 12%}{e^{0.12\cdot 10}-1}=25860.77</math>
| |
| | |
| At this point there is a temptation to simply divide by 12 to obtain a monthly payment. However this would contradict the primary assumption upon which the "continuous payment" model is based: namely that the annual payment ''rate'' is defined as:
| |
| | |
| :<math>M_a=\lim_{N\to\infty}N\cdot x(N) \,</math>
| |
| | |
| Since it is of course impossible for an investor to make an infinitely small payment infinite times per annum, a bank or other lending institution wishing to offer "continuous payment" annuities or mortgages would in practice have to choose a large but finite value of ''N'' (annual frequency of payments) such that the continuous time formula will always be correct to within some minimal pre-specified error margin. For example hourly fixed payments (calculated using the conventional formula) in this example would accumulate to an annual payment of 25861.07 and the error would be < 0.02%. If the error margin is acceptable, the hourly payment rate can be more simply determined by dividing ''M''<sub>''a''</sub> by 365×24. The (hypothetical) lending institution would then need to ensure its computational resources are sufficient to implement (when required) hourly deductions from customer accounts. In short cash "flow" for continuous payment annuities is to be understood in the very literal sense of the word.
| |
| | |
| : ''"Monies paid into a fund in the financial world are paid at discrete – usually equally spaced – points in calendar time. In the continuous process the payment is made continuously, as one might pour fluid from one container into another, where the rate of payment is the fundamental quantity".''<ref>[[#CFP|Beckwith]] p. 114.</ref> | |
| | |
| The following table shows how as ''N'' (annual compounding frequency) increases, the ''annual'' payment approaches the limiting value of ''M''<sub>''a''</sub>, the annual payment ''rate''. The difference (error) between annual payment and the limiting value is calculated and expressed as a percentage of the limiting value.
| |
| | |
| {| border="1" cellpadding="2"
| |
| | align="center" style="background:#f0f0f0;"|'''Compounding Period'''
| |
| | align="center" style="background:#f0f0f0;"|'''Frequency (N)'''
| |
| | align="center" style="background:#f0f0f0;"|'''Per period interest rate'''
| |
| | align="center" style="background:#f0f0f0;"|'''Per period payment x(N)'''
| |
| | align="center" style="background:#f0f0f0;"|'''Annual Payment'''
| |
| | align="center" style="background:#f0f0f0;"|'''% Error'''
| |
| |-
| |
| | Bi-annual||2||6.000000%|| 13,592.28|| 27,184.56||5.118918%
| |
| |-
| |
| | Quarterly||4||3.000000%|| 6,631.19|| 26,524.76||2.567558%
| |
| |-
| |
| | Monthly||12||1.000000%|| 2,173.55|| 26,082.57||0.857683%
| |
| |-
| |
| | Daily||365||0.032877%|| 70.87|| 25,868.07||0.028227%
| |
| |-
| |
| | Hourly||8760||0.001370%|| 2.95|| 25,861.07||0.001176%
| |
| |}
| |
| | |
| <ref>Further worked examples and problems with solutions can be found in Professor Hackman's course notes. See Reference Section.</ref><ref>Beckwith (pages 128–129) provides more complex examples involving interest rate calculation. The interested reader may verify the calculations by entering the resultant transcendental equations on Wolfram Alpha. Note: the line of working before eqn (38) in Beckwith's article is missing a pair of brackets</ref>
| |
| | |
| It will be apparent from the above that the concept of a "continuous repayment" mortgage is a somewhat theoretical construct. Whether it has practical value or not is a question that would need to be carefully considered by economists and actuaries. In particular the meaning of an annual repayment ''rate'' must be clearly understood as illustrated in the above example.
| |
| | |
| However the "continuous payment" model does provide some meaningful insights into the behaviour of the discrete mortgage balance function – in particular that it is largely governed by a [[time constant]] equal to the reciprocal of r the nominal annual interest rate. And if a mortgage were to be paid off via fixed daily amounts, then balance due calculations effected using the model would – in general – be accurate to within a small fraction of a percent. Finally the model demonstrates that it is to the modest advantage of the mortgage holder to increase frequency of payment where practically possible.
| |
| | |
| ==Summary of formulae and online calculators==
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=b8acb9f3726fce014595187a18ed05b9 Annual payment rate] (mortgage loan): <math>M_a=\lim_{N\to\infty}N\cdot x(N)=\frac{P_0\cdot r}{1 - e^{-rT}}</math>
| |
| | |
| Annual payment rate (sinking fund): <math>M_a=\lim_{N\to\infty}N\cdot x(N)=\frac{P_T\cdot r}{e^{rT}-1}</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=a99957c3b83e7d54e27d6f22d47b16f7 Future value:] <math>F_v(t) = \frac{M_a}{r}(e^{rt}-1)</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=7507dafaedac784a18c852536bbd3c88 Present value:] <math>P_v(t) = \frac{M_a}{r}(1 - e^{-rt})</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=695553c0c5ba292284a3e1e64da6ed47 Loan balance:] <math>P(t) = \frac{M_a}{r}(1 - e^{-r(T-t)})</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=cdad431591373e47eba5ab6c984c2b55 Loan period:] <math>T=-\frac{1}{r}\ln\left(1-\frac{P_0 r}{M_a}\right)</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=8e887cf8e64ab8e7173701a979476567 Half-life of loan:] <math>t_{\frac{1}{2}}=\frac{1}{r}\ln\left(\frac{1+e^{rT}}{2}\right)</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=65982e894fcde21153454b2ea4b2a8a Interest rate:] <math>r\approx\frac{2}{T}\ln{\frac{M_aT}{P_0}}</math> <math>r=\frac{1}{T}\left (W(-se^{-s})+s\right )\text{ with }s=\frac{M_at}{P_0}</math>
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=6c11cb78b7bbb5c22d5f5271b5494381 Universal mortgage calculator]. Given any three of four variables, this calculates the fourth (unknown) value.
| |
| | |
| [http://developer.wolframalpha.com/widgets/gallery/view.jsp?id=a2b56600097fa2ea81bfc9774f1baaf Mortgage graph]. This illustrates the characteristic curve of mortgage balance vs time over a given loan timespan. Loan amount and loan interest rate (''p''/''a'') may also be specified. A discrete interval loan will have a very similar characteristic.
| |
| | |
| ==Notes==
| |
| {{Reflist}}
| |
| | |
| == References ==
| |
| * {{Cite journal
| |
| |ref = CFP
| |
| |last = Beckwith
| |
| |first = R.E.
| |
| |date = June 1968
| |
| |title = Continuous Financial Processes
| |
| |jstor = 2329786
| |
| |journal = The Journal of Financial and Quantitative Analysis
| |
| |volume = 3
| |
| |issue = 2
| |
| |pages = 113–133}}
| |
| * {{Cite web
| |
| |ref = GIT
| |
| |author = Georgia Institute of Technology
| |
| |coauthors = Hackman, Steve
| |
| |date =
| |
| |title = Financial Engineering: ISyE 4803A Course Notes
| |
| |url = http://www.ontario-mortgage-rates.ca/FixedRateMortgages.pdf
| |
| |publisher = Georgia Institute of Technology
| |
| |accessdate = 2009-04-27
| |
| }}
| |
| * {{Cite book
| |
| |ref = MJG
| |
| |last = Glencross
| |
| |first = M.J.
| |
| |year = 2007
| |
| |title = Mathematics: Grade 12 OBE
| |
| |publisher = Effective Teaching Publishers, Cape Town, RSA.
| |
| |ISBN = 978-1-920116-36-1
| |
| }}
| |
| * {{Cite book
| |
| |ref = A&T
| |
| |last = Munem
| |
| |first = M.A.
| |
| |coauthors = Foulis D.J.
| |
| |year = 1986
| |
| |title = Algebra and Trigonometry with Applications.
| |
| |url = http://books.google.co.za/books?id=1cOyXr0W6XIC
| |
| |publisher = Worth Publishers, USA.
| |
| |ISBN = 0-87901-281-1
| |
| }}
| |
| * {{Cite book
| |
| |ref = BDH
| |
| |last = Hahn
| |
| |first = Brian D.
| |
| |year = 1989
| |
| |title = Problem Solving with True Basic.
| |
| |url = http://books.google.co.za/books?id=OQDhAAAACAAJ&dq
| |
| |publisher = Juta & Company Limited, Cape Town, South Africa.
| |
| |ISBN = 0-7021-2282-3
| |
| }}
| |
| | |
| == Bibliography ==
| |
| * Kreyszig, Erwin, ''Advanced Engineering Mathematics'' (1998, Wiley Publishers, USA), ISBN 0-471-15496-2.
| |
| | |
| {{DEFAULTSORT:Continuous-Repayment Mortgage}}
| |
| [[Category:Mortgage]]
| |
| [[Category:Mathematical finance]]
| |