Continuous-repayment mortgage

src: www.tsb.co.uk

Analogous to continuous compounding, a continuous annuity is an 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. 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₀, per period interest rate i, number of periods n and fixed per period payment x, the end of term balancing equation is:

${\displaystyle P_{0}(1+i)^{n}=\sum _{k=1}^{n}x(1+i)^{n-k}={\frac {x[(1+i)^{n}-1]}{i}}}$

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₀, annual interest rate r, loan timespan T (years) and annual rate M_a, the infinitesimal cash flow elements M_a?t accumulate continuously compounded interest from time t to the end of the loan timespan at which point the balancing equation is:

${\displaystyle P_{0}e^{rT}=\int \limits _{0}^{T}M_{a}e^{r(T-t)}\,dt={\frac {M_{a}(e^{rT}-1)}{r}}.}$

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.

Within the timespan of the loan the time continuous mortgage balance function obeys a first order linear differential equation (LDE) 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).

Video Continuous-repayment mortgage

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:

${\displaystyle P_{v}(n)={\frac {x(1-(1+i)^{-n})}{i}}}$

The formula may be re-arranged to determine the monthly payment x on a loan of amount P₀ taken out for a period of n months at a monthly interest rate of i%:

${\displaystyle x={\frac {P_{0}\cdot i}{1-(1+i)^{-n}}}}$

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:

${\displaystyle x(N)={\frac {P_{0}\cdot r}{N(1-(1+{\frac {r}{N}})^{-NT})}}}$

${\displaystyle N\cdot x(N)={\frac {P_{0}\cdot r}{1-(1+{\frac {r}{N}})^{-NT}}}}$

Note that N·x(N) is simply the amount paid per year - in effect an annual repayment rate M_a.

It is well established that:

${\displaystyle \lim _{N\to \infty }\left(1+{\frac {r}{N}}\right)^{Nt}=e^{rt}}$

Applying the same principle to the formula for annual repayment, we can determine a limiting value:

${\displaystyle 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}}}.}$

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:

${\displaystyle P_{v}(N,t)={\frac {N\cdot x(N)(1-(1+{\frac {r}{N}})^{-Nt})}{r}}}$

Applying the limiting expression developed above we may write present value as a purely time dependent function:

${\displaystyle P_{v}(t)=\lim _{N\to \infty }P_{v}(N,t)={\frac {M_{a}}{r}}(1-e^{-rt})}$

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:

${\displaystyle P(t)={\frac {M_{a}}{r}}(1-e^{-r(T-t)})={\frac {P_{0}(1-e^{-r(T-t)})}{1-e^{-rT}}}}$

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 here.

Maps Continuous-repayment mortgage

Comparison with similar physical systems

Define the "reverse time" variable z = T - t. (t = 0, z = T and t = T, z = 0). Then:

${\displaystyle P(z)={\frac {M_{a}}{r}}(1-e^{-rz}).}$

This may be recognized as a solution to the "reverse time" differential equation:

${\displaystyle {\frac {M_{a}}{r}}={\frac {1}{r}}{\frac {dP(z)}{dz}}+P(z).}$

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.

${\displaystyle V_{0}=RC{\frac {dV(t)}{dt}}+V(t).}$

The key characteristics of such equations are explained in detail at 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.

src: www.americanfinancing.net

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_N (divided into N equal payments made at time intervals ?t where ?t = 1/N years), we may write:

${\displaystyle {\begin{aligned}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{aligned}}}$

If N is increased indefinitely so that ?t -> 0, we obtain the continuous time differential equation:

${\displaystyle {\operatorname {d} P(t) \over \operatorname {d} t}=rP(t)-M_{a}}$

Note that for there to be a continually diminishing mortgage balance, the following inequality must hold:

${\displaystyle P_{0}\leqslant {\frac {M_{a}}{r}}}$

P₀ 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:

${\displaystyle P_{t+\Delta t}=P_{t}(1+r\Delta t)-M_{N}\Delta t\;}$

Using the notation P_n to indicate the mortgage balance after n periods, we may apply the recursion relation iteratively to determine P₁ and P₂:

${\displaystyle P_{1}=P_{0}(1+r\Delta t)-M_{N}\Delta t\;}$

${\displaystyle {\begin{aligned}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{aligned}}}$

It can already be seen that the terms containing M_N form a geometric series with common ratio 1 + r? t. This enables us to write a general expression for P_n:

${\displaystyle {\begin{aligned}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{aligned}}}$

Finally noting that r ? t = i the per-period interest rate and ${\displaystyle M_{N}\Delta t=x}$ the per period payment, the expression may be written in conventional form:

${\displaystyle P_{n}=P_{0}(1+i)^{n}-{\dfrac {x[(1+i)^{n}-1]}{i}}}$

If the loan timespan is m periods, then P_m = 0 and we obtain the standard present value formula:

${\displaystyle P_{0}={\dfrac {x[1-(1+i)^{-m}]}{i}}}$

Solving the differential equation

${\displaystyle {\operatorname {d} P(t) \over \operatorname {d} t}=rP(t)-M_{a}}$

One method of solving the equation is to obtain the Laplace transform P(s):

${\displaystyle P(s)={\frac {M_{a}}{s(r-s)}}={\frac {M_{a}}{r}}\times {\frac {(-r)}{s(s-r)}}.}$

Using a table of Laplace transforms and their time domain equivalents, P(t) may be determined:

${\displaystyle P(t)={\frac {M_{a}}{r}}(1-e^{rt}).}$

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:

${\displaystyle {\begin{aligned}&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{aligned}}}$

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_n in the difference equation, the expression for P(t) may be written in the following algebraically equivalent form:

${\displaystyle P(t)=P_{0}e^{rt}-{\frac {M_{a}}{r}}(e^{rt}-1)}$

src: iannellisbakery.com

Calculation of accumulated interest and principal payments

Re-arranging the original differential equation we obtain:

${\displaystyle M_{a}=rP(t)-{\frac {dP(t)}{dt}}}$

Integrating both sides of the equation yields:

${\displaystyle M_{a}.t=\int _{0}^{t}rP(t)\,dt\,-\int _{0}^{t}{\frac {dP(t)}{dt}}\,dt\,}$

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_at. Evaluating the first integral on the right we obtain an expression for I(t), the interest paid:

${\displaystyle I(t)=M_{a}t-{\frac {M_{a}(e^{rt}-1)}{re^{rT}}}}$

Unsurprisingly the second integral evaluates to P₀ - P(t) and therefore:

${\displaystyle I(t)=M_{a}t-P_{0}+P(t)\,}$

The reader may easily verify that this expression is algebraically identical to the one above.

src: www.americanfinancing.net

Loan cost factor

The cost of a loan is simply the annual rate multiplied by loan period:

${\displaystyle C=M_{a}T={\frac {P_{0}rT}{1-e^{-rT}}}}$

Let s = rT. Then we may define loan cost factor C(s) such that C = P₀C(s) i.e.: C(s) is the cost per unit of currency loaned.

${\displaystyle C(s)={\frac {rT}{1-e^{-rT}}}={\frac {s}{1-e^{-s}}}}$

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₀rT (rT >> 0).

By way of example, consider a loan of 1000000 at 10% repaid over 20 years. Then s = 0.1 × 20 = 2.

${\displaystyle C=1000000\times {\frac {2}{1-e^{-2}}}\approx 2.313\times 10^{6}}$

The product rT is an easily obtained but important parameter in determining loan cost according to the equation C=P₀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.

src: iannellisbakery.com

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_e where s_e=r_et and r_e is the equivalent simple interest rate:

${\displaystyle {\frac {s}{1-e^{-s}}}=1+s_{e}}$

It is straightforward to determine s_e 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_e.

In his book Problem Solving with True Basic, 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.

src: gustavoroubert.com

Period of a loan

If a borrower can afford an annual repayment rate M_a, then we can re-arrange the formula for calculating M_a to obtain an expression for the time period T of a given loan P₀:

${\displaystyle {\begin{aligned}&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{aligned}}}$

src: gallerlaw.com

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:

${\displaystyle k={\frac {M_{\min }}{M_{a}}}={\frac {P_{0}r}{M_{a}}}}$

Now we may consider a small re-arrangement of the equation for loan period T:

${\displaystyle T=-{\frac {1}{r}}\ln \left(1-{\frac {P_{0}r}{M_{a}}}\right)}$

${\displaystyle rT=s(k)=-\ln(1-k)\;}$

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).

src: bloximages.newyork1.vip.townnews.com

"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:

${\displaystyle {\frac {P(t)}{P_{0}}}={\frac {1}{2}}={\frac {1-e^{-r(T-t)}}{1-e^{-rT}}}}$

Solving for t we obtain:

${\displaystyle t_{\frac {1}{2}}={\frac {1}{r}}\ln \left({\frac {1+e^{rT}}{2}}\right)}$

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.

src: iannellisbakery.com

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:

${\displaystyle P_{0}={\frac {M_{a}}{r}}(1-e^{-rT})}$

we may write:

${\displaystyle rP_{0}=M_{a}(1-e^{-rT})\,}$

${\displaystyle \Rightarrow rP_{0}-M_{a}+M_{a}e^{-rT}=0}$

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₀, M_a 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:

${\displaystyle f'(r)=0\,}$

${\displaystyle \Rightarrow P_{0}-M_{a}Te^{-rT}=0}$

${\displaystyle \Rightarrow r={\frac {1}{T}}\ln {\frac {M_{a}T}{P_{0}}}}$

Since the function is approximately parabolic between the roots at r = 0 and the sought value, we may estimate the required root as:

${\displaystyle r\approx {\frac {2}{T}}\ln {\frac {M_{a}T}{P_{0}}}}$

Using this as a starting point, increasingly accurate values for the root may be determined by repeated iterations of the Newton-Raphson algorithm:

${\displaystyle r_{1}=r_{0}-{\frac {f(r_{0})}{f'(r_{0})}}.\,\!}$

Some experimentation on Wolfram Alpha reveals that an exact analytical solution employing the Lambert-W or "product log" function can be obtained. Setting s = M_aT/P₀ we obtain:

${\displaystyle r={\frac {1}{T}}\left(W(-se^{-s})+s\right)}$

In the region of interest W(-se^-s) 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 analytical solution using the Lambert W or "productlog" function on Wolfram Alpha.

Newton-Raphson iterations

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:

${\displaystyle P_{v}(t)={\frac {M_{a}}{r}}(1-e^{-rt}).}$

In similar fashion, a future value formula can be determined:

${\displaystyle F_{v}(t)={\frac {M_{a}}{r}}(e^{rt}-1).}$

In this case the annual rate M_a is determined from a specified (future) savings or sinking fund target P_T as follows.

${\displaystyle 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}}.}$

It will be noted that as might be expected:

${\displaystyle F_{v}(t)=P_{v}(t)\times e^{rt}.\,}$

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):

${\displaystyle P(t)=P_{0}e^{rt}-{\frac {M_{a}}{r}}(e^{rt}-1).}$

Example

The following example from a school text book 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:

${\displaystyle x(12)=PMT(1\%,120,500000)=2173.55}$

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:

${\displaystyle M_{a}={\frac {500000\times 12\%}{e^{0.12\cdot 10}-1}}=25860.77}$

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:

${\displaystyle M_{a}=\lim _{N\to \infty }N\cdot x(N)\,}$

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_a 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".

The following table shows how as N (annual compounding frequency) increases, the annual payment approaches the limiting value of M_a, 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.

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

Annual payment rate (mortgage loan):       ${\displaystyle M_{a}=\lim _{N\to \infty }N\cdot x(N)={\frac {P_{0}\cdot r}{1-e^{-rT}}}}$

Annual payment rate (sinking fund):       ${\displaystyle M_{a}=\lim _{N\to \infty }N\cdot x(N)={\frac {P_{T}\cdot r}{e^{rT}-1}}}$

Future value:       ${\displaystyle F_{v}(t)={\frac {M_{a}}{r}}(e^{rt}-1)}$

Present value:       ${\displaystyle P_{v}(t)={\frac {M_{a}}{r}}(1-e^{-rt})}$

Loan balance:       ${\displaystyle P(t)={\frac {M_{a}}{r}}(1-e^{-r(T-t)})}$

Loan period:              ${\displaystyle T=-{\frac {1}{r}}\ln \left(1-{\frac {P_{0}r}{M_{a}}}\right)}$

Half-life of loan:       ${\displaystyle t_{\frac {1}{2}}={\frac {1}{r}}\ln \left({\frac {1+e^{rT}}{2}}\right)}$

Interest rate:           ${\displaystyle r\approx {\frac {2}{T}}\ln {\frac {M_{a}T}{P_{0}}}}$              ${\displaystyle r={\frac {1}{T}}\left(W(-se^{-s})+s\right){\text{ with }}s={\frac {M_{a}t}{P_{0}}}}$

Universal mortgage calculator. Given any three of four variables, this calculates the fourth (unknown) value.

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

References

• Beckwith, R.E. (June 1968). "Continuous Financial Processes". The Journal of Financial and Quantitative Analysis. 3 (2): 113-133. JSTOR 2329786. 
• Georgia Institute of Technology; Hackman, Steve. "Financial Engineering: ISyE 4803A Course Notes" (PDF). Georgia Institute of Technology. Retrieved 2009-04-27. 
• Glencross, M.J. (2007). Mathematics: Grade 12 OBE. Effective Teaching Publishers, Cape Town, RSA. ISBN 978-1-920116-36-1. 
• Munem, M.A.; Foulis D.J. (1986). Algebra and Trigonometry with Applications. Worth Publishers, USA. ISBN 0-87901-281-1. 
• Hahn, Brian D. (1989). Problem Solving with True Basic. 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.

Source of the article : Wikipedia