Numerical Methods in Computational Finance. Daniel J. Duffy

Чтение книги онлайн.

Читать онлайн книгу Numerical Methods in Computational Finance - Daniel J. Duffy страница 34

Numerical Methods in Computational Finance - Daniel J. Duffy

Скачать книгу

solution is given by:

upper U Superscript n Baseline equals upper A Superscript n Baseline upper U 0 plus upper B StartFraction 1 minus upper A Superscript n Baseline Over 1 minus upper A EndFraction comma n greater-than-or-equal-to 0

      where we note that upper A Superscript n Baseline identical-to n Superscript italic t h power of constant upper A and upper A not-equals 1.

      In order to prove this, we need the formula for the sum of a series:

1 plus upper A plus ellipsis plus upper A Superscript n Baseline equals StartFraction 1 minus upper A Superscript n plus 1 Baseline Over 1 minus upper A EndFraction comma upper A not-equals 1 period

      For a readable introduction to difference schemes, we refer the reader to Goldberg (1986).

      Contrary to popular thinking, there is more to ODEs than trying to find analytical solutions for them. Very few ODEs have analytical solutions, and we must resort to ODE solvers.

      The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in essence the same as a small and concrete special case.

      Paul Halmos

      In Chapter 2 we discussed both systems of ODEs and scalar ODEs. The focus was mainly concerned with notation, the structure of ODEs and finite difference schemes to approximate them. We implicitly assumed that the solution of the corresponding initial value problem existed in an otherwise unspecified time interval and that the solution was unique. These assumptions constitute a huge leap of faith. In this chapter we discuss existence and uniqueness results for ODEs and stochastic differential equation (SDEs). We also introduce several important numerical schemes and code in C++ and Python.

      We turn our attention to a more general initial value problem for a non-linear system of ODEs:

      where:

StartLayout 1st Row 1st Column Blank 2nd Column y colon normal double struck upper R right-arrow normal double struck upper R Superscript n Baseline comma upper A element-of normal double struck upper R Superscript n Baseline comma f colon normal double struck upper R times normal double struck upper R Superscript n Baseline right-arrow normal double struck upper R Superscript n Baseline 2nd Row 1st Column Blank 2nd Column and colon 3rd Row 1st Column Blank 2nd Column f left-parenthesis t comma y right-parenthesis equals left-parenthesis f 1 left-parenthesis t comma y right-parenthesis comma ellipsis comma f Subscript n Baseline left-parenthesis t comma y right-parenthesis right-parenthesis Superscript down-tack Baseline where f Subscript j Baseline colon normal double struck upper R times normal double struck upper R Superscript n Baseline right-arrow normal double struck upper R comma j equals 1 comma ellipsis comma n period EndLayout

       Does System (3.1) have a unique solution?

       In which interval does this solution exist?

       What is the asymptotic behaviour of the solution as ?

      To this end, let B be a region of n plus 1 dimensional space, and let f be continuously differentiable with respect to t and with respect to all the components of y at all points of B. We assume that the following inequalities hold:

      and:

Скачать книгу