Quasilinearization numerical scheme for fully nonlinear parabolic problems with applications in models of mathematical finance

Abstract In this paper, on the basis of Newton’s method, we propose a fast quasilinearization numerical scheme, coupled with Rothe’s method, for fully nonlinear parabolic equations. General conditions that provide quadratic, uniform and monotone convergence of the quasilinearization method (QLM) of solving fully nonlinear ordinary differential equations that arise on each time level, are formulated and elaborated. The convergence of QLM and its rate are examined numerically, on a simple test example with an exact solution. The first few iterations already provide extremely accurate and stable numerical results. The second goal is to consider three applications of the proposed schemes in financial mathematics. Namely, numerical results for three nonlinear problems of optimal investment are presented and discussed. The numerical experiments of the last problem are based on the data from statistic information of the Bulgarian National Bank and Bratislava Interbank.

[1]  Martin Schechter.  Periodic solutions of second-order nonautonomous dynamical systems , 2006 .

[2]  R. Bellman,et al.  Quasilinearization and nonlinear boundary-value problems , 1966 .

[3]  Naoyuki Ishimura,et al.  Existence of solutions for the nonlinear partial differential equation arising in the optimal investment problem , 2008 .

[4]  V. Lakshmikantham,et al.  Generalized quasilinearization versus Newton's method , 2005, Appl. Math. Comput..

[5]  Adam M. Oberman,et al.  Pricing early exercise contracts in incomplete markets , 2003, Comput. Manag. Sci..

[6]  A mathematical theory for numerical treatment of nonlinear two-point boundary value problems , 2006 .

[7]  G. Barles,et al.  Convergence of approximation schemes for fully nonlinear second order equations , 1991 .

[8]  Matthias Ehrhardt,et al.  Nonlinear models in option pricing --- An introduction , 2010 .

[9]  I. Blank,et al.  Convergence of Rothe's method for fully nonlinear parabolic equations , 2004 .

[10]  A. Samarskii.  The Theory of Difference Schemes , 2001 .

[11]  Ansgar Jüngel,et al.  Convergence of a High-Order Compact Finite Difference Scheme for a Nonlinear Black-Scholes Equation , 2004 .

[12]  Tanya G. Melton,et al.  Generalized quasilinearization method and higher order of convergence for second-order boundary value problems , 2006 .

[13]  Werner C. Rheinboldt,et al.  On Discretization and Differentiation of Operators with Application to Newton’s Method , 1966 .

[14]  Astronomy,et al.  Quasilinearization approach to nonlinear problems in physics with application to nonlinear ODEs , 2001, physics/0102041.

[15]  V. B. Mandelzweig,et al.  Numerical investigation of quasilinearization method in quantum mechanics , 2001 .

[16]  Tetsuro Yamamoto,et al.  Historical developments in convergence analysis for Newton's and Newton-like methods , 2000 .

[17]  Daniel Sevcovic,et al.  On the risk-adjusted pricing-methodology-based valuation of vanilla options and explanation of the volatility smile , 2005 .

[18]  Daniel Sevcovic,et al.  Weakly Nonlinear Analysis of the Hamilton-Jacobi-Bellman Equation Arising from Pension Savings Management , 2009, 0905.0155.

[19]  Martin Stynes,et al.  Numerical Treatment of Partial Differential Equations , 2007 .

[20]  I. Stakgold.  Green's Functions and Boundary Value Problems , 1979 .

[21]  Matthias Ehrhardt,et al.  On the numerical solution of nonlinear Black-Scholes equations , 2008, Comput. Math. Appl..

[22]  L. Kantorovich,et al.  Functional analysis and applied mathematics , 1963 .

[23]  Enrique A. Navarro,et al.  Numerical solution of linear and nonlinear Black-Scholes option pricing equations , 2008, Comput. Math. Appl..

[24]  J. Pratt.  RISK AVERSION IN THE SMALL AND IN THE LARGE11This research was supported by the National Science Foundation (grant NSF-G24035). Reproduction in whole or in part is permitted for any purpose of the United States Government. , 1964 .